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CXX. The Interaction in the Theory of Beta Decay. 


By D. L. Pursry*, 
Department of Natural Philosophy, Glasgow University}. 


[Received July 5, 1951.] 


ABSTRACT. 


The cross terms occurring in the expression for the theoretical nth 
forbidden beta spectrum are given for an arbitrary mixture of interactions. 
Relations are derived between the six first forbidden vector nuclear 
matrix elements. Comparison with experiment is made for first forbidden 
transitions with allowed shape beta spectra. The results are discussed 
in relation to other recent work in the theory of beta decay. It is 
concluded that the interaction is most probably tensor with little or no 
admixture of other interactions. 


§1. InTRODUCTION. 


Tue Fermi type theory of beta decay postulates an interaction between 
heavy and light particle fields proportional to a linear combination of 
the five Lorentz invariant expressions : 


S=(¥"*BD)(*B4), 
V=(¥*0)(f*4)—(V*ab)(p*ad), 


T= (P*Bo®)(b*Bop) + (P*Ba®) (b*Bag), (1) 
A= (P*o)(b*ad)—(P*y5P)(p* ys), 
P=(P*By5)($*By 5); J 


where ¥Y, &, %& and ¢ are the proton, neutron, electron and neutrino 
fields respectively, and a, o, 8 and y, are the usual Dirac matrices. 
There is no good theoretical reason for preferring any particular linear 
combination of these interactions over any other, and consequently 
the interaction is free to be fitted with experiment. The large amount 
of accurate experimental information which has become available during 
the past few years not only has firmly established this type of theory 
for allowed and certain forbidden transitions, but helps considerably 
towards the experimental determination of the interaction. 

Because of the arguments favouring the Gamow-Teller selection 
rules (summarized by Konopinski 1943), the interaction has long been 
thought to contain at least a proportion of either T or A. The most 
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convincing evidence for this now available, however, is the experimental 
discovery of several beta spectra having the characteristic shape, 
predicted only by T and A, for first forbidden transitions with an angular 
momentum change of two units and parity change. An excellent survey 
of recent experimental work, including the results just mentioned, has 
been given by Wu (1950). 

With a mixture of interactions, cross terms may occur which affect 
the predicted spectrum shape. For allowed transitions, such cross terms 
occur only for mixtures of S with V or of A with T, and have the effect 
of multiplying the allowed spectrum by a factor (1--a/W), where W 
is the beta particle energy in units of mc?. These cross terms have been 
calculated by Fierz (1937), who neglected the Coulomb field of the 
nucleus, and more recently by De Groot and Tolhoek (1950). The latter 
authors give an upper limit of 4 for a, and conclude that mixing of 
S with V or of A with T, if it occurs at all, must be very slight. 

It is clear from the foregoing that to be consistent with recent 
experimental evidence on forbidden transitions the interaction must 
contain either T or A, while evidence on allowed spectra forbids 
interactions containing both T and A, or both S and V. There remain 
unconsidered, however, a fair number of transitions which, although 
forbidden according to the semi-empirical ft classification, nevertheless 
have allowed shape beta spectra. Some of these transitions may be 

_interpreted as first forbidden with an angular momentum change of 
0 or 1 unit and change of parity, and they will be used here to gain 
further information about the beta decay interaction. 


§ 2. CoRRECTION Factors, 


Greuling (1942, equations (11)-(16)) has given the correction factors 
by which the allowed spectrum must be multiplied to give the theoretical 
nth forbidden spectrum for each of the five interactions. For this paper, 
however, the correction factors associated with cross terms are also 
required. To obtain these, the nth forbidden correction factor for an 
arbitrary mixture of interactions was calculated by the method given 
by Konopinski and Uhlenbeck (1941). The correction factors for the 
cross terms are rather cumbersome, and to avoid interrupting the 
argument, they are given in an appendix. These cross terms when used 
in conjunction with Greuling’s formule give the complete nth forbidden 
correction factor for an arbitrary mixture of interactions, except for a 
few terms unimportant except for J=0 to J=0 transitions, and a few 
terms which can be regarded as small corrections to the (n—1)th forbidden 
correction factor. As an example, the nth forbidden correction factor 
for the mixture A+yS is C,,+y?C,,.+-yC,(S, A), where C,4 and C,. are 
the correction factors given by Greuling for the pure axial vector and 
pure scalar, interactions respectively, and C,,(S, A) is the cross term given 
in the appendix. 
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§ 3. Tue PsEuposcaLaR INTERACTION. 


By convention, P in this paper is regarded as having no allowed 
transitions. It is then seen from Greuling’s formule, and from the cross 
terms given in the appendix, that the selection rules for an nth forbidden 
transition according to the interaction P are the same as those for an 
(n—2)th forbidden transition according to the interactions T and A. Hence 
the correction factors will be almost independent of the amount of P in 
the interaction unless this is great enough to overcome the extra 
forbiddenness. If this were the case, one would expect the first forbidden 
pseudoscalar matrix element of P to give apparently superallowed 
transitions ; all the observed superallowed transitions, however, with the 
exception of B', are already well explained by Wigner’s theory of 
supermultiplets (see, for example, Konopinski 1943), and for none of 
them does a parity change seem reasonable. Furthermore, if the third 
forbidden P matrix element were strong enough to affect first forbidden 
transitions with unit angular momentum change, it would certainly 
dominate the much weaker transitions with angular momentum change 
of two units, and this does not appear to be the case. For these reasons, 
it is concluded that, while the amount of P in the interaction may be 
quite large, it is not sufficient to affect the correction factors for first 
forbidden transitions with unit angular momentum change, The 
correction factor for first forbidden transitions with no angular momentum 
change, however, may well be modified or even dominated by the presence 
of P in the interaction. 


§ 4. Nuc“teaR Matrix ELEMENTS. 


In the first forbidden correction factor, terms occur involving twelve 
different matrix elements. Of these, two are the second order tensors 
responsible for the transitions with angular momentum change of two 
units; their contribution to other first forbidden transitions is small, 
and will be ignored. Four of the matrix elements are pseudoscalars, 
and therefore cannot contribute to transitions involving an ‘angular 
momentum change. In this section, relations will be found between the 
remaining six matrix elements, all of which are vectors. 

The notation (f|a|i) will be used for these matrix elements, instead 
of the notation fa used by Konopinski and Uhlenbeck (1941); by this 
is to be understood (for electron emission) 


(fai) )== [art . . , facts 3 a0 (evil ale gee aie (2) 


where Y, and Y; are respectively the nuclear wave functions for the 
jnitial and final states, a is a charge independent vector operator, 
7, and 7, are the usual isotopic operators, the superscripts refer to the 
particles in the nucleus, and the integration includes summation over 
spin and isotopic spin states. As is usual in the theory of beta decay, 


4M2 
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that representation of the Dirac matrices is chosen in which f and 
a« become in non-relativistic approximation —1 and —v/c respectively, 
while the units of action, velocity and mass are chosen to be f, c and m,. 

In the absence of any relativistic model of the nucleus the results of 
this section must be regarded as essentially non-relativistic, although 
relativistic notation will be used for the sake of convenience. It is 
therefore permissible to replace 8 by its non-relativistic approximation 
—l, provided it commutes with the other operators appearing in the 
matrix element, and hence 


(fer |)=—(f|r]t); (f[Boar|t)=—(floar]i). . . - (3) 


Further relations can be derived using the Hartree type of single 
particle model. A sufficient condition that the single particle states 
should be eigenstates of the operator k=f[(o.m)+1], where m is the 
orbital angular momentum, is that the magnitudes / and j of the orbital 
and total angular momentum of the states should both be good quantum 
numbers. This condition is certainly satisfied by the shell model of 
Mayer (1949, 1950) and of Haxel, Jensen and Siiss (1949, 1950). The 
operator k is the well-known angular momentum operator of relativistic 
theory (see, for example, Dirac 1947), and has eigenvalues 


(je). ee 


The following relations, in which [] and {} denote commutators and 
anti-commutators respectively, are easily verified : 


(4) 


[k, r] =—iBoar, 
(5) 
[k?, a]= {k, Ba}, 
[k, [k?, r]] = {k, r}, 
(6) 
[k, [k?, a]]= {k, a}, 


(6) expresses the selection rules k;+-k,=0 or k;—ky=+1 for the matrix 
elements of r and a; from (4) they are seen to be equivalent to 
Aj=0, +1, yes. (5) gives the relations 


(f[Boar|i)=—te(f|r|4); (f[Balt)=—e(flalé), . . . (7) 


where «=k;—k,y. If k;+-k,=0, the second of these relations is true only 
in non-relativistic approximation. These relations, together with (3), 
enable the six vector matrix elements to be reduced to two, namely 
(f|r|é) and (f|a|i). 

The selection rules for (f|o|i) are found similarly to be k,=k, or 
k,+k,;=+1, while, provided k;k;, 


(f[Bo l= +e)floli. . 2 2 1. . (8) 
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Since 8 commutes with co, it may be replaced in non-relativistic 
approximation by —1. Consequently the matrix element for transitions 
with k;+k,= +1 must vanish in non-relativistic approximation, and 
must therefore be, of order (v/c)?; indeed, it is clear that (8) is self- 
consistent for such transitions only if the large components of the wave 
functions do not contribute to the matrix element. Obviously, transitions 
of this type will appear to be forbidden, and they must not be confused 
with genuinely first forbidden transitions, all of which involve a parity 
change. From (4) it is seen that the selection rule k,-+k=+1 is 
equivalent to 4j= +1, 4J/J= +2; for the allowed transitions with 
k,+k,= —1 or k=k,, Al=0. 

To obtain a relation between (f|« |i) and (f|r|i), the nuclear Hamiltonian 
is assumed to be 


Ha — 2 {(@. pi)-+40U(1+7)M, +(1 7M} 


14 A 


to 2 2 {Vole +(e! cH)V jk) —(o! +o). (2°28), (pV ik) 
ah VG, aire al. ee we ba art (9) 


where V,, V,, V, and V, are the ordinary, charge exchange, spin-orbit 
and Coulomb potentials between two nucleons respectively, and 
M,, and M,, are the neutron and proton masses. As before, the use of 
relativistic notation is for convenience only, and this Hamiltonian is 
to be considered valid only in non-relativistic approximation. For the 
sake of simplicity, and in the absence of any evidence to the contrary, 
the spin-orbit coupling is supposed nh to be of charge exchange type. 


5 5 (73 —i7})r/ with H, it is found that 
j=l 


Fle iW Wolf e194 104, Mf 1Be19 
Base gat ol { de P(r) —ir}) (o8(2) —2®)V (Jk) 


g=1k9 
~i(oi +o") a (ri—2*)V (jk) —H— eV (Gk) }M%. (10) 


The last term in this expression will be averaged over the particles k, 
assuming they form a spin saturated spherical core of uniform density 
and radius p equal to the nuclear radius. 

For both the exchange force and spin-orbit coupling contributions 
to (10), the problem is essentially that of averaging (r7/—r")V(7), where 
p= [r/—r* |, over the core formed by the particles k. If V(r) is expanded 
in a series of Legendre polynomials of the cosine of the angle between 
ri and r* and the spherical symmetry of the core is used, this averaging 
becomes elementary, and gives 


By calculating the commutator se 


(ci—r)V(7) )=1(3/16r;p ule V(r)[(p+7;)2—2?][r?—(p—1;)?|rdr, . (11) 


yl 


1198 D. L. Pursey on the 


where 7,= [r?|. Since the transforming particle is presumably in an outer 
shell, 7; may be taken approximately equal to p. For a Yukawa well 
of depth v and range 1/k, the averaging then gives 
riv(3A/2k%p%){(kp)-2—3(kp)-2-+ exp (—2kp)[2+5(kp)-? + 6(kp)- 2+ 3(kp)*]} 
(12 a) 
(Hughes and Le Couteur (1950) have also found this result, using a 
rather different method), while for a square well of depth v and range 
1/k, one obtains 
rip(ZA/16k4p4)[1—3(kp)-*)]. . , . - - (126) 
In the charge exchange contribution to (f|a|¢), the factor 7§ gives on 
averaging a factor (N—Z)/A where N is the neutron number of the initial 
nucleus and Z the proton number of the product nucleus. V, is assumed 
to be of the form [a,-+a,(o/ . 6”)]V(jk) but because of the spin saturation 
of the core, only the term with coefficient a, will contribute ; for the same 
reason, any tensor forces of the usual type will vanish. Using the value 
given by Rosenfeld (1948) for a,, based on the ratio of the triplet and 
singlet s potentials and the requirement of saturation for heavy nuclei, 
the effective V, is found to be —(1/10)V(?s) where V(?s) is the triplet 
s potential. Taking the nuclear radius p to be 7A, and using both 
a Yukawa well of depth 100mc? and range 7, and a square well of depth 
40mc? and range 27, for V(’s), the exchange force contribution to 
(fla|¢) is found to be 


Yukawa well : —15iA-48(N—Z)(1—3A~*8)(f|r|z).  . . (13a) 
Square well : —12iA-48(N—Z)(1—ZA—*8)(f|r|z).  . . (136) 


The two wells give results in fair agreement. This exchange contribution 
is so small in any case that the use of more accurate well parameters is 
unnecessary. 

Hughes and Le Couteur (1950) estimate that for V, a Yukawa well 
of range ry and depth 1 MeV. is consistent both with the probable 
spin-orbit separation in He® and with the separation necessary according 
to the spin-orbit coupling shell model to account for the apparent level 
order in heavy nuclei. (The factor 2 difference’ in the depth given above 
from that given by Hughes and Le Couteur is due to a different definition 
of c; in this paper the spin angular momentum is $6). On averaging, 
the contribution to (f|«|%) from the spin-orbit coupling is found to be 


—3A~19(1—3A~°)(f lo ar|i)=—3icA-¥8(1—3A—28)(f|r|i) (14) 
on using (7). 7 
The Coulomb contribution to (f|a|i) can be estimated by a similar 


averaging procedure, and is found to be u(Za/p)(f|r|¢). Collecting these 
results, one obtains 


(fla |2)=tx(f |r|?) 
where 
«=(W,—W,)—(M,, —M,,)-+(Za/p)—15A-4(N—Z)(1—3A-28) 


—3<A-48(1 — 34-213), oo 
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For most first forbidden transitions, the first two terms in approximately 
cancel, while the exchange force and spin-orbit terms are small compared 
to the Coulomb term. For positron emitters, the second, third and 
fourth terms in x change sign. 


§ 5. COMPARISON WITH EXPERIMENT. 


It was pointed out in the introduction that the interaction cannot 
contain both S and V, or both T and A, while the pseudoscalar interaction 
P was discussed in §3 and found not to contribute to first forbidden 
transitions involving an angular momentum change. To avoid considera- 
tion of P, attention will be confined to transitions in which an angular 
momentum change of one unit for the transforming single particle state 
seems plausible according to the shell model of Mayer and of Haxel, 
Jensen and Stiss. Consequently the only mixtures that need be discussed 
are those of S with T, A with V, V with T and S with A. The first two 
of these have almost identical correction factors, and therefore only the 
first of them will be considered. The correction factors for these mixtures 
and the type of transition under consideration may be found from the 
general formule of Greuling together with the cross terms given in the 
appendix, and may be further simplified using the relations between 
matrix elements found in the last section. In low Z approximation, they 
are 


T+yS : C= {2+ (2/3W)+ (n/18) +5} 
— 2ye{€[(Za/2p)+ (2/3)] + (Wob/9)+ (Za/p)/(6W)} 
ODS CASEI StS Oma iemabite 
l6a 
T+yV: C= (1+y?) {E+ (2€/3W)+ (9/18) +5} +-y(n/6) 
— 2ye {(€2/W) + (2£/3)+-(Wo/9)(2—W,/W)}. . . . (166) 
A+yS : C=(1+y"){(Za)2p)2+ (Za|p)(£/3)+ (n/6)—3(W»—W)(W—1/W)} 
NT SAEED eS leat 
1l6¢ 
~ where €=x— (Za/2p)—W,/3, 
n= W2—1—2W,W-+2W?, OTe tee (1'7) 
¢=2W—W,-—1/W, 


a is the ratio of (f|a|i) to i(f|r]2), and e=k,—k,= +1, depending on the 
transition. A factor |(f|r|¢) |? has been omitted from all of these correction 
factors. For positron emitters, the sign of Z must be changed throughout, 
and the sign of the cross term (proportional to y) changed in the mixture 
T+y8. 

Only transitions with an ft value greater than 10% and whose experi- 
mental Kurie plot is straight to within 5 per cent for a reasonable energy 
range will be regarded as sufficiently forbidden and _ sufficiently 
accurately investigated to warrant detailed comparison with theory. 
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A search in the literature for such transitions has yielded the following : 
Na%2 (Good, Peaslee and Deutsch 1946), Na** (Siegbahn 1946), Sc* 
(Peacock and Wilkinson 1948), V* (Peacock and Deutsch 1946), 
Pr146 (Feldman, Lidofsky, Macklin and Wu 1949), Pm!47 (Agnew 1950, 
Langer, Motz and Price 1950, Lidofsky, Macklin and Wu 1949), 
Hf18! (Chu and Wiedenbeck 1949), Re1** (Beach, Peacock and Wilkinson 
1949, Langer and Price 1949), W187 (Peacock and Wilkinson 1948), 
Au$ (Langer and Price 1949, Saxon 1948, Steffen, Huber and Humbel 
1949). P82, Ga®® and Cs184 are more probably of the 4j=1, d/=2 type 
of transition mentioned in the last section, while Mo% and I1%! are not 
easily interpreted on the shell model picture; these five interactions are 
therefore not considered. ; 

The sign of ¢ is obtained from the shell model interpretation of the 
transition. The level scheme given by Haxel, Jensen and Siiss (1950) 
is used in conjunction with the measured spins of neighbouring odd 
nuclei to suggest the initial state, and, if the transition is between ground 
states, also the final state, of the transforming nucleon. If the transition 
is to an excited level of the product nucleus, a reasonable transition 
involving parity change must be postulated. In this case, for heavy 
nuclei with a large neutron excess the transforming neutron need not 
belong to the outermost occupied neutron shell, and several interpretations 
may be possible. 

For Na*? and Na*4, the shell model interpretation is complicated by 
the anomalous filling of the d;). shell; however, for Na” the most 


probable transition seems to he p3i2—>81;2, for which «= —1, while for 
Na** the transition is probably d;.—fz. with «= +1. For both 
Se** and V48, the only reasonable transition is f,.—d5)., with e= —1. 


Pr'# and Pm!’ undergo ground state to ground state transitions, and 
are almost certainly of the type hg2—>gz. with e= +1. For Hf18}, 
Ret and W1%? the transition may be any of t13/2 >hr1/2) fs)2 >9 722 
P1232, frj2—dsj, all with e= —1, or hyig>G7/2 With e= +1, or 
fu2972 11 which case the transition is not of the required type ; for 
these transitions « was assumed to be —1. For Au?®8, the transition is 
probably f52>d3). with e= +1. | 
The correction factors (16a) and (166) are not very sensitive to the 
value of x unless € is small. This is the case for Na®2, and fig. 1 shows 
theoretical Kurie plots for this transition according to the tensor inter- 
action for different values of x. (15) gives for 2 the value —2-5, and this 
gives an almost straight Kurie plot in fair agreement with experiment. 
The Possibility of error, however, cannot be denied, and in particular, 
if « is +1 instead of —1 as assumed, (15) gives c= —3-9 and a Kurie 
plot far from straight. It seems unlikely, therefore, that any reliance 
can be put on conclusions about mixtures involving T or V based on 
this transition. It is also found that for the positron emitters Na22 and 
v* the Kurie plots for all mixtures of the type A+yS are reasonably 
straight in the experimentally measured region. Consequently, no 
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reliable information about any mixtures of interactions can be obtaine d 
from the available evidence on Na®, and this transition will no longe r 
be considered. For all the other. transitions, the theoretical spectrum 
shapes are not very sensitive to the value of x, and the estimate of 
§ 4 should give fairly reliable results. 

For heavy nuclei with fairly low end points, the spectrum shapes 
predicted by all the pure interactions are indistinguishable from the 
allowed shape. For mixtures, however, if the cross term is negative 
there may be a large cancellation causing the theoretical spectrum to 
deviate strongly from the allowed shape. This is illustrated for 


Fig. 1. 


+Experimental Points 
(After Good et al) 


2:2 (mc*) 


Theoretical Kurie plots for Na?? according to the tensor interaction, with 
different values for the ratio « of —i(f|a|i) to (f|r}{i). The estimate of §4 
gives x= —2°5. 


Pr143 in fig. 2. For all mixtures, if ye=—1 the theoretical Kurie plot 
is similar to that shown for the mixture A+S, while if y increases or 
decreases, the Kurie plot approaches a straight line. The mixtures 
T-+yS and V+yA are remarkable in that the theoretical Kurie plot 
changes from ‘‘ concave down” through an S-shape to “ concave up ”’ 
before becoming straight ; this change is very rapid, and for Pr!®* it 
is complete when y=0-8. Frequently the experimental Kurie plot shows 
an excess of low energy electrons which may be due to experimental 
difficulties such as source thickness effects, or to a second low energy ~ 
transition, or it may be genuine. As a precaution which should be more 
than adequate to prevent false conclusions being drawn, only that part 
of the Kurie plot which is experimentally straight is compared with 
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theory. Theoretical and experimental Kurie plots are regarded as 
consistent with each other if the discrepancy between them is not more i 
than 5 per cent at the low energy end and 10 per cent at the high energy 7} 
end of the straight portion. For Pr‘, only that part of the Kurie plot 1 
above 1:4 mc is used, and the Kurie plots for the mixtures A+3S and 
T-+4S both show the maximum allowed discrepancy. Since the sign of 
the cross term depends on « and therefore on the transition, it clear ; 
that by considering different transitions both lower and upper limits 
can be obtained for the ratio y of the imteraction constants mM any 
mixture. a 

.< For lighter nuclei such as Na**, the picture is very similar, except 
that the predicted spectrum shapes for pure A and S are no longer 


Fig. 2. 
VPWF +Experimental Poiats 
(After Feldman et al) 


ww 
I +2 4 6 -8 20 22 a+ 26 28 @,5 
Theoretical Kurie plots for Pr™* for various mixtures of interactions. 


allowed, and indeed for Na** no mixture whatsoever of S with A gives 
an allowed shape spectrum. Fig. 3 gives theoretical Kurie plots for 
Na** for the interactions A, A—S and T. of which only the last is 
consistent with experiment. It was pointed out in the introduction 
that the interaction must contain either T or A, but not both: the 
spectrum of Na** seems to favour T, but A with a large admixture of ; 
V cannot be ruled out. Of the nine transitions under consideration, ; 
Na** is the only one capable of distinguishing between T and A in this : 
way, and it should be kept in mind that the Na™ result is valid only if © 
the transition is genuinely first forbidden and not a highly unfavoured — 
allowed transition, and if the first forbidden pseudosealar matrix 
elements do not contribute. This will be the case whatever the angular 
momentum change of the nucleus as a whole. provided. as was assumed 


a ee, ee, 


—_—/- 


> 
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above, the single particle transition is of the type d52—>f72; if, however, 
the transition is, for example, d,.—f,., the result will be valid only 
if the nucleus as a whole changes its angular momentum. 

The limits to the amount of mixing found from each of the nine 
transitions considered are given in Table I. The first-named interaction 
in each mixture is assumed to predominate, so that y is numerically less 
than one; since, however, the mixtures StyA and V+ ¥V are almost 
identical with the mixtures A+yS and T+wyV respectively, and since 
the Jimits are given for both T+yS and S+yT, this means no loss of 
generality. The mixtures V+yA and A+yV are almost identical to 
the mixtures T+yS and S+yT respectively. 


Fig. 3. 
As 
pwr 


: + Experimental Points 
(After Siegbahn) 


b5 2 25 3 ts a oe 


Theoretical Kurie plots for Na * for various interactions. No mixture of the 
type A+ is consistent with experiment. 


From Table I. it is seen that very little admixture of S with A, and 
little admixture of V with T is permissible. Mixtures of T and §, or of V and 
A, are possible, but one or other of the two interactions must predominate. 
If the evidence of Na**4 is accepted, then the only permissible interactions 
are T with perhaps a small admixture of either V or S, and a mixture of 
A and V with V predominating. The tensor interaction T is the only 
pure interaction consistent both with the interpretation given for 
Na‘ and with all the other experimental data on beta spectrum shapes. 

The validity of these results, of course, depends on the accuracy of the 
estimates used for the ratios of the nuclear matrix elements, and on the 
correctness of the shell model interpretations of the transitions. An 
error in the relations between the matrix elements will cause a 
proportionate error in some or all of the limits given in Table L. ; such an 
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error, however, would have to be very large before it would alter the 
main features of the table. An error in the shell model interpretation 
may cause the adoption of the wrong sign for «; in this case, the limits 
given will be true for —y instead of for y. 


TABLE I. 
Upper and Lower Limits for the Mixture Parameter y. 


Mixture Na?4 Sct6 ys Pri48 Pm!47 Hf181 Rel8é* Wwis7 Au198 


The first-named interaction is assumed to predominate. The upper and lower 
rows of figures for each mixture are respectively upper and lower limits for y 
while the columns indicate the transitions from which the limits are derived. 
X denotes that no values of y are consistent with experiment. 


* See note added in proof. 


§ 6. CONCLUSION. 

In the introduction, it was pointed out that because of the arguments 
in favour of the Gamow-—Teller selection rules, and because of the discovery 
of spectra having the characteristic first forbidden spin change of two 
shape, the interaction must contain either T or A. From the allowed 
spectrum shapes, it was seen that the interaction cannot contain both 
S and V, or both T and A. Taken in conjunction with these results 
the work of this paper shows that the interaction must either be rien 
bi with possibly a small admixture of S or V, or mainly V with an 
admixture of A. If the evidence of Na24 is not accepted a predominantl 

axial vector interaction’ is also possible. The decays of RaE of (isi 
are almost certainly second forbidden, and their spectrum shapes appear 
to demand an interaction in which either T or V creda 
(Konopinski and Uhlenbeck 1941, Langer and Moffat 1951), though not 
all possible mixtures have been considered. By considering the lifetimes 
of superallowed transitions, and assuming Wigner’s theory of super- 
multiplets, Moszkowski (1951) has concluded that the interaction is 
puhee almost entirely T or almost entirely A. Marshak (1942) and 
Greuling (1942) have calculated theoretical half-lives for various forbidden 
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transitions using the tensor interaction, and found reasonable agreement 
with experiment. It appears therefore that the experimental evidence 
strongly favours the pure tensor. interaction, although none of the above 
evidence precludes an admixture of P. 

Tiomno Wheeler and Rau (1949) and Tiomno and Wheeler (1949) 
have shown that the beta decay and the capture by protons of y-mesons 
can both be explained by quadrilinear interactions with the same 
interaction constant as in nuclear beta decay, and this has led to the idea 
of a universal quadrilinear interaction. For such an interaction, symmetry 
properties with respect to the interchange of fields having similar physical 
properties may have physical significance. The best known interaction 
having such properties is the Critchfield interaction S-A—P (Critchfield 
1943), which, if the fields are quantized so that any two anti-commute, 
is symmetric with respect to the interchange of any pair of fields. Apart 
from the pure interactions, all of which are either symmetrical or anti- 
symmetrical with respect to the interchange of the light particle or the 
heavy particle fields, the only two other interactions having special 
symmetry properties which are consistent with allowed beta spectra 
are S—T+P and 38+T+3P, which are respectively symmetrical and 
antisymmetrical with respect to the interchange of the two neutral or 
the two charged fields. Of these three interactions, only the last is 
consistent with the results of Table I, and it is inconsistent with 
Moszkowski’s results, and probably also with the RaE and Cs!*’ spectra. 
On changing the order of the charged and neutral fields in the interaction 
to that used by Michel (1950) in calculating the theoretical beta spectrum 
from p-meson decay, 38+T+3P becomes pure tensor, which predicts 
a spectrum in complete disagreement with experiment. Consequently 
@ universal interaction, if one exists, can have no special symmetry 
properties, and the ordering of the fields in the interaction must 
therefore be uniquely determined by their physical properties such as 
charge, or sign of magnetic moment of a neutral particle. This need not 
preclude the hypothesis of a Majorana neutrino, provided there are some 
neutral particles which are not of Majorana type; this is, however, 
the case, since, both because of its magnetic moment and because of the 
stability of nuclei, the neutron cannot be a Majorana particle. By 
comparing either with Michel’s results or with the figures given by 
Tiomno, Wheeler and Rau, it is seen that the tensor interaction with the 
ordering of the fields usual in the theory of nuclear beta decay is consistent 
with the beta spectrum from p-meson decay, whether the neutrino is 
a Dirac or a Majorana neutral particle. 
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APPENDIX. 


The correction factors for the cross terms occurring in the expression 
for the theoretical nth forbidden beta spectrum for an arbitrary mixture 
of interactions are given below, using Greuling’s (1942) notation except 
that Q,,(a) is written for Greuling’s Q,(a, r)/(m!). The units of action, 
velocity and mass are chosen to be #, ¢ andm, respectively. The method 
of using these correction factors was explained in §2. 


Cy(8, V=—[Q,(BF)QK(e)+e.0] E (A,,Q,9-9-24D, Pg} 


+1[Q,,(Br)Qi(a)—c.c.] y fA, ,R,g—-)-2-LC, ,P,g2ir--1}, 
v=0 
C,(8, T)=i1Q, (Br) Q4(Ba)—c.c.] F {A,,N,g2- 24.6, Lge} 
v=0 


+7Q,,(Br)Q*(Boar)—c.c.] z {A ,M,q2-")-2 
v=0 


+(Dy.—B,,)b,g%"-}; 


C nl A)=—{Q, (Br) n(GAL)—C. c.] z {A,,Q, te "2 —2C,,,R,g"-"-1 


+(Dy,—By,)P,g?r-}, 


C,,(V, ye [Q,,(o Q? *(Bar) )+e. c.] z= A P g2"-”)-2 


—1Q,, (x) Q; n (BO) )—¢. c.] z= as Anke a oe 24-C,, ak Br picid 7 


—1Q,(1)Qi(Boar)—c.c.] z tn, g-9-24 20, R,gur-na 


+ (D, rage eas! 


1, (a)Qir(Boar)+e.c.] 2 {A Bg'- "24.6, ,Pger-o-1}, 


C,( (V, A)= [Q,(a a)Q*(oar)+c.c.] z {A,,,N,q%"-210, L q2in—)-1} 
y=0 ¥ 


+7[Q,,(r)Q*(oar)—c.¢.]- z {A ,M,g2(n-"-2 
v=(0 


ny 


Be (D,, ee B, Liga. 


C,.(V, P)=—7Q,_ i(aar)Q Cie 1(By5r) )—c.c.] = ac Q,qrn-*-b-2 


m1,» 


—2C,-, ,R,q%-"-1)-1 ; 
n-1,» of +(D ‘ge sists ep cate 
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C,(T, A)= —[Qn4 1(Ps)Q ea o)+¢.¢.] a etree eS 


—[Q, (Boar) Seni aaey z {A,, ,Q,q%-"-2 


— —— 2(n—”) 
oF (0, DNR, “ 1 By } Pg } 


+[Q,(Ba)Q*(oar)-+e.c.] ¥ {A,, ,B,g%"--24 C,, ,P,q2—)-1}, 
v=0 

C,,(T, P)=—1[Q,,(8e)Q*(Bysr)—c.c.] 2 {A,,N,g--2+.6, qu}, 
e 2—() 


C,(A, P)=i[Q,,(6)Q¥(Bysr)—c.c.] Z {A,,R,g2"--246, ,P gr}, 
»v=0 


These last two correction factors represent a mixing of the 
(n+1)th forbidden P spectrum with the (n—1)th forbidden T and 


A spectra respectively. 
The quantities P,, Q, and R, are defined analogously to L,, M, and 


N, by 
Py (27°F 5) 1935S otalp 


Q,= (2p7F y)-4[9? (+-2)—f 2] 2», 
R,= (2p7F 9) “1,9 +f -G+29-O+ oh aia 


By expanding the radial functions g and f in powers of p and neglecting 
the first and higher powers, one obtains 


2! 2/s,(v-+1+s,)\ 1 
=P (se) Capa) Ww 
aor NS wee Za\? 1 (Lat (La 
Q,=(F,/Fo) (ori) {- ae Ae a) int ae & ) 


S,(v-+1) as, |e 
oe Gs 1) ak gi eo ae Ww? 


(4s +3)(Za)? Za \2 
a5 bee y+l1-+s, a | (sa) wh. 


pentrara (Sahn) {5 (+r 


For low Z, s, may be put equal to v-+1, F,=Fo, and (Zx)’=0 to give the 
simpler expressions 


2p | 201 | oy Cae Ne 1 =) (ya 
oo Gears We ae) {5 T \ar3) Wi 


2’ 2 Lo \ ae 
eg care Cys 
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i i i d Tolhoek 
For positron emitters, the sign of Z and, as De Groot an 
(1950) have shown, the signs of C, (8, V), C,(8, 1) Cas ADCs Ee 
C,,(V, P) and C,(T, P) should be changed. 


Note added in proof—Since this paper was written, papers by Metzger 
and Hill (1951) and Steffen (1951) have been received confirming the 
result of Grant and Richmond (1949) that the spectrum of Re? is 
complex. In calculating Table I., this spectrum was assumed to be 
simple: consequently the limits in Table I. derived from the Re1%é 
transition are narrower than experiment now warrants. 
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CXXI. Ultrasonic Propagation in Liquid Helium near the Lambda-Point. 
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ABSTRACT. 


It is suggested that the sharp minimum in velocity and maximum in 
attenuation of ultrasonic waves in liquid helium in the region of the 
X-point is a consequence of fluctuations which lead to the production 
of inclusions of He II in He I above the A-point, and of He I in He II 
below the A-point. Expressions are derived for the velocity and 
attenuation in an inhomogeneous fluid, and are applied, with some 
modifications, to liquid helium. It is found that the shape of the velocity 
curve may be accounted for satisfactorily, and that the observed 
attenuation just above the A-point requires the inclusions of He II to 
consist of about 850 atoms. It is shown that fluctuation theory predicts 
the presence of inclusions of this size in the concentration required to 
explain the experimental results. 


§1. InTRODUCTION. 


EXPERIMENTS on the propagation of ultrasonic waves in liquid helium 
have revealed peculiarities of behaviour at temperatures near the A-point 
(2-18° K.); the velocity drops to a sharp minimum at the A-point, and 
at the same time the attenuation rises rapidly to a peak. The effect which 
culminates at the A-point begins to be apparent not far below 3° K., at 
which temperature the specific heat also displays the first symptoms of 
its rise to a sharp maximum, so that it is probable the two effects are 
intimately related. Naturally the specific heat variation is reflected in 
a corresponding variation of y, the ratio of the principal specific heats, 
which drops to a value near unity at the A-point, and this provides a 
satisfactory phenomenological explanation of the shape of the velocity 
curve. It does not, however, in itself suggest a mechanism for the rapid 
increase in attenuation, for which a number of tentative proposals have 
been made, none of which has been worked out fully enough to carry 
conviction. The present paper is intended as a crude analysis of a model 
of liquid helium above the A-point which provides a qualitative explanation 
of the observed attenuation. A complete theoretical investigation would 
be of far too great a complexity to be practicable at present. 

It has been suggested by Keesom and Keesom (1935) that the increase 
in the specific heat of HeI between 3:0° K. and 2-18° K. is to be 
attributed to thermodynamical fluctuations of sufficient magnitude to 
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SER. 7, VOL. 42, NO. 334.—NOV. 1941 4¥ 


1210 A. B. Pippard on the Ultrasonic Propagation 


lead to the formation of inclusions of He II in the liquid, and, although 
their analysis is not sufficiently refined to give more than a rough estimate 
of the importance of this effect, there can be little doubt of its essential 
correctness. As a starting-point, therefore, we shall assume that in this 
temperature range the liquid is inhomogeneous in the sense that it 
contains small inclusions whose thermodynamical properties, if not 
exactly the same as those of He II in bulk, are at least markedly different 
from those of HeI. Their influence on y may be readily understood by 
considering the reduction which they bring about in the temperature rise 
during adiabatic compression. For the present purpose it is convenient 
to think of the effect of adiabatic compression as equivalent, so far as 
temperature changes are concerned, to the injection of heat at constant 
pressure. The inclusions are now seen to exert a two-fold influence. In 
the first place they reduce the temperature rise by increasing the mean 
specific heat of the liquid, and in the second place they actually reduce 
the equivalent heat input ; for while in order to simulate the temperature 
changes accompanying adiabatic compression heat must be injected 
into He I, it must be removed from He II on account of the negative 
expansion coefficient of the latter. The net result is that the temperature 
rise of the liquid is' less than it would be in the absence of inclusions, and 
correspondingly y and the velocity of sound are reduced. 

As mentioned above, a detailed analysis of the processes taking place 
in a real liquid subject to fluctuations during passage of a sound wave 
would be a very difficult matter to treat rigorously, and we shall therefore 
in the next section consider a model of an inhomogeneous liquid from 
which the major theoretical obstacles have been eliminated. We shall 
then be enabled to discuss qualitatively the extent to which the behaviour 
of this model parallels that of liquid helium, and hence to demonstrate 
that the fluctuation theory is capable of accounting for the general 
features observed experimentally. 


§2. PRopaGaTION OF SOUND IN AN INHOMOGENEOUS MEDIUM. 


Consider a medium consisting of a number of spherical inclusions of 
uniform size embedded at random in a matrix of different properties, and 
initially at the same temperature. In order to determine the propagation 
characteristics for a compressional wave of angular frequency w in this 
medium it is necessary to calculate the adiabatic compressibility. It is 
clear from the previous section that if the temperature rise of the 
inclusions during compression is less than that of the matrix, and if there 
is thermal equilibrium between the two, the compressibility of the 
medium as a whole will be affected by the presence of the inclusions. 
Moreover, if w is so high that thermal equilibrium is not perfect, the 
volume oscillations of the medium will lag behind the pressure oscillations, 


and there will be a continuous conversion of the mechanical energy of the 
wave into heat. ‘ 
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If a fraction a of the medium is in the form of inclusions, whose 
compressibility under the conditions of the experiment is k,, while the 
compressibility of the matrix is k,, the mean compressibility will be given 
by the expression : 


EA eae nat rah tas Lae FL) 
The actual compressibility of the matrix, k,, is related to its adiabatic 
compressibility, k,,, by the expression : 
dT 
Risen icant sigs tates eatin eee) 


where £, is the expansion coefficient, 1/V (0V/@T),, of the matrix, and 
dT,/dp is the rise of temperature of the matrix as a result of compression 
under the conditions of the experiment. Since for adiabatic compression 


of the matrix, 
(=) ania! E 
CD sere 


we may rewrite (2) in the form, 


kaka =Yi— (V4 ee ea ew (8) 
where h,= {dT,/dp}/(0T,/0p),. Similarly for the inclusions 
keo/Kog=Ve—(Ye—1 hg. ee aes ae) 


If h, and h, vary from point to point they must be replaced by their 
average values h, and h,. Thus the determination of k is resolved into 
the problem of determining h, and hy. 

' As explained earlier, we may simulate the effects of adiabatic 
compression by injecting heat; for a given rate of change of pressure 
let the appropriate rate of injection of heat be q, per unit volume of 
matrix and g, per unit volume of inclusions. It is a simple matter to 
show thermodynamically that 


J2/91=P2/Pi=A (say). 


If the heat is injected and withdrawn periodically with frequency w, the 
resulting temperature oscillations of the matrix and inclusions would 
be represented, in the absence of thermal contact, by the equations 


Noes Lay and Vaiss Lo==G5, 


in which T, and T, are the departures from the mean of the temperatures 
of matrix and inclusions respectively, and C, and C, are specific heats 
per unit volume. As is common in the solution of periodic problems, 
d/dt is replaced by tw. 

In order to calculate the effects of thermal transfer it is convenient 
to divide the heat sources into two parts, which may be superposed to 
give the complete solution. Let us first supply heat at a rate q, per unit 
volume of matrix and pq, per unit volume of inclusion, where »=C,/C). 


ANZ 
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This we shall refer to as the standard state, in which the temperature 
oscillation is the same, g,/twC,, for both matrix and inclusions, without 
the necessity of heat exchange between the two. To represent the actual 
state of affairs we must in addition supply heat to the inclusions at a 
rate q,(A—p) per unit volume. We may first consider one inclusion 
only to be supplied with this additional heat, the rest being in the 
standard state, and solve the resulting conduction problem. Afterwards 
the result may be summed over all inclusions, since all the effects are 
additive. 

If the material outside the one inclusion considered be smoothed out 
into a uniform material of specific heat C(=(1—x)C,+2aC,), the 
conduction problem is straightforward. The equation of heat conduction 
in the material surrounding the inclusion takes the form, for spherical 
symmetry, 

PL 2 dle 400 

dr 'rdr « 
where x is the thermal conductivity ; and the appropriate solution, in 
which T tends to zero as r tends to infinity, may be written : 


ce “exp (—fr), where £=,/(iwC/x). 


The integration constant, u, is determined by the boundary condition 
at the surface of the inclusion, where r=a. For simplicity we shall 
suppose the inclusion to be at a uniform temperature at all times, the 
same as that of the matrix immediately surrounding it, that is [T],. 
Then since heat is supplied to the inclusion at a rate $7a%q,(A—) from the 
exterior source, and at a rate 4za?«x[dT/dr], by conduction from the 
matrix, the rate of rise of temperature of the inclusion, iw[T],, is governed 
by the equation 


3 
—* [aT /dr},—ieoCg[T],+9,(A—n)=0. 


From this equation w may be determined, and the amplitude of the 
temperature oscillation of the inclusion is found to be given by the 
expression 
oe a® q1(A—p) 2 
[Tt] = ne 1h or ances ; where t=C,a pl: eo (5) 
Thus of the heat ¢,(A—p) supplied to unit volume of the inclusion, the 
amount which remains within it is given by twC,[T],, 2. e. 


q1(A—p) 
Sear cents: + Soke heh pe 


and the rest, which is conducted to the surrounding material is given by 
91(A—pz)—iwC,[T],, ¢. e. 


twr 


(Ia). are MPN es 
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Now the surrounding material is made up of matrix and inclusions in the 
standard state, and the heat (7) will be divided among them in the 
proportions : 
ie XL 
;_—_——.— to the matri eee 
1+(u—1)a 7 I=) 
Let us now add together the effects of all the inclusions, acting in the 
same way. In unit volume of material there will be a volume a occupied 
by inclusions, so that, combining (7) and (8), we see that the heat supplied 
to the matrix by all the inclusions amounts to 


to the inclusions.* « (8) 


QA—z) _ x12) 
O48) Thal-piwe Te” 


or, since the total volume of matrix is (l—z2), 
q1(A—z) ou bat 
ltal+tiwr 14(u—l1)a 


per unit volume of matrix. To this we must add the heat supplied to 
unit volume of matrix in the standard state, which is just q,, to give for 
the resultant total rate of heat supplied per unit volume to the matrix 


Aa fi- , tS 
Pam e (1-Fal-Fies)f1+(n—l)a]f* 


Similarly we may calculate the resultant total rate of heat supply 
per unit volume to the inclusions 


(1—x)(1-+-al)(u—A) 
(fal fiorI+@—Nals 


Now di lqy=hy and G3 /Wo=ho ; hence from (1), (3) and (4), 


(1-+-aZ) 


. l 
ad, 1+ i! 


ate —2)kyg+thog 


_ &(1—ax)(u—A)(1+al) ; 
Been ic Ge Tie] [he ya(¥1— 1) —kealy2— V/A]; 


con (1) u—A)E- a1 = Dk 
=U Vira Maat pal pion u— le] 


since, by a well-known thermodynamical formula, y—1=/?T/Ck. 

From this expression for the mean compressibility it will be seen that 
at sufficiently high frequencies only the first two terms are important, 
and this clearly corresponds to such rapid oscillations of pressure that 
no heat is exchanged between the inclusions and the matrix. The last 
term expresses the modification brought about by heat exchange at 


(9) 


* This is true only if there is thermal equilibrium between the matrix and 
inclusions, 7. e. at low frequencies, but the error involved at higher frequencies 
should not be great, provided that x is not too large. 
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lower frequencies. When w=0, at and wr vanish, and for frequencies 
low enough for af and w7 to be much smaller than unity, (9) may be 
rewritten in a slightly simpler, approximate form : 


a(1—2x)(u—A)*(y —1)k a . 
, k+(1—a\hyat thyet ——— aye? ree ali) 


Since the velocity of sound, v=1//(kp), where p is the density, the 
velocity and attenuation of a sound wave may be immediately deduced 
from (9) or (10). If v is the limiting velocity at high frequencies, 
corresponding to a compressibility given by the first two terms of (10), 
and Av is written for v)—v, then from (10), 


Av, &(1—a)(u—AP1— Nha 


Ate te Ar eg i ee 11 
% > Dall + — Ne — a) thea] (it) 
and the attenuation per unit distance, 
ae TABU APY = Dos i 


2401+ (u—1)a] 


This expression for « is, of course, valid only when w7<1; as the 
frequency is increased « rises as w? at first, but more slowly as wr 
approaches unity, and finally, when w7>1, it increases only as 1/w so 
that the attenuation per wavelength actually decreases. If k,, and ky, 
are not greatly different, equations (11) and (12) may be combined to 
yield a simple expression for « at frequencies for which w7<l, 


_wtAv  w*C,a* dv 
Sing ee (13) 

Before discussing the relation between the model we have analysed 
and the actual state of affairs in liquid helium, we may take the argument 
one stage further. Up to now we have assumed that the medium contains 
permanent inclusions, but it would more nearly approximate to a 
fluctuating liquid if we considered the inclusions as transient, appearing 
and disappearing in a random manner at all points in the medium. 
As will appear below, this transience of the inclusions does not modify 
equation (11) for the velocity, but it may reduce the attenuation. In order 
to understand the effect let us return to an early stage in the previous 
calculation, and consider one inclusion surrounded by the material in 
the standard state. The inclusion oscillates in temperature with an 
amplitude [T], given by (5), in addition to the oscillation of the standard 
state. Thus the actual amplitude under these conditions is 


q 
tw, 


(l—tw7’), if wr’ <1, where 7’ =7(1—A/p). 


The oscillation of the inclusion is thus of very nearly the same amplitude 
as that of the rest of the material, but with a time-lag of 7’. 
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Now consider a small region of the material, which, if the inclusions 
are transient, will spend a fraction x of the time as part of an inclusion, 
and the rest (l—a), as part of the matrix. So long as it is part of the 
matrix its temperature will stay in phase with that of the standard state, 
but as soon as it becomes part of an inclusion it will tend to lag behind 
the standard state, and if it remains part of an inclusion for a time long 
compared with 7’ it will assume the time-lag +’ corresponding to 
dynamical equilibrium. When, with the disappearance of the inclusion, 
it becomes once more part of the matrix, the time-lag will revert to zero 
with a time constant approximately equal to 7’. Thus if the time of 
persistence of the inclusions is much greater than 7’, the mean time-lag 
of the inclusions will be the same as if they were permanent, and the 
results of the previous calculation will apply. On the other hand, if the 
inclusions persist only for a time shorter than 7’ they will not take up the 
full time-lag before they disappear, their mean time-lag will be less than 7’ 
and the attenuation will be correspondingly reduced. For any given 
time of persistence of the inclusions, their mean time-lag, f7’ (say), where 
f<i, may in principle be calculated, and the previous arguments may be 
carried through exactly as before, with + everywhere replaced by fr. 
Clearly equation (11) for the velocity will be unaltered, since it does not 
involve 7, but (12) and (13), expressing the attenuation, will be modified, 
and a new approximation to the attenuation coefficient may be written : 


_ w*fC,a* Av 
eg ———. 


14 
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We shall return later to a discussion of the value to be assigned to /.. 


§ 3. APPLICATION TO Liguip HELIUM. 


The model which we have analysed in some detail resembles the actual 
state of affairs in liquid helium only superficially, and the analysis has 
been given mainly in order to make clear the nature of the processes 
which may affect the velocity and attenuation in an inhomogeneous 
fluid. The major difference between a fluctuating liquid and the model 
lies in the fact that in the former the inclusions exist essentially as a 
result of temperature fluctuations, so that it is not possible to consider 
the initial state of the material, before the passage of a sound wave, as one 
of uniform temperature. In spite of this difference, however, a close 
analogy may be drawn between the two systems, and to a large extent 
the fluctuations may be disregarded except in so far as they are the 
mechanism responsible for the appearance of inclusions. This may be 
understood by reference to a new model which more closely resembles 
a real fluctuating liquid. Let us suppose that the liquid is initially 
entirely free from spontaneous fluctuations and is uniform in temperature 
and in all other properties ; now let heat be injected into and withdrawn 
from the liquid from point to point in a random manner, in such a way 
that the total internal energy remains constant. The temperature of the 
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liquid will now fluctuate randomly, and the energy flow resulting from 
its thermal conductivity will be similar to that which occurs spontaneously 
in a real liquid. If the mean temperature is not too much higher than 
the A-point the local withdrawals of heat may be sufficient to produce 
transient inclusions of He II. Now let us superimpose on this random 
heat injection an orderly injection which will simulate the effect produced 
by adiabatic compression; that is, we inject heat at a uniform rate 
into every region of the matrix (He I), and at a different, but uniform, 
rate into every inclusion (He II). If the amount of heat injected in this 
way is very much smaller than the random injections its influence on the 
distribution of the inclusions will not be very great, though, as we shall ~ 
see below, it cannot be entirely neglected. Since the equation of heat 
conduction is linear the energy flow due to each cause is independent, 
so that the ordered heat flow is the same as if matrix and inclusions were 
initially at the same temperature. Moreover, the effect of each type of 
energy flow on the volume of the liquid is independent, so that each 
contributes independently to the compressibility. The random fluctua- 
tions, however, do not affect the volume, so that they may be disregarded. 
To this degree of approximation, therefore, the real liquid may be treated 
as exactly analogous to the model which we have studied, and the 
expressions for the velocity and attenuation of a sound wave may be 
taken over directly. . i 
There is, however, one assumption in the preceding argument which is 
not strictly valid, namely that the passage of a sound wave does not 
affect the distribution of the inclusions. In fact the heating of the material 
by compression shortens the lifetime of each inclusion in such a way as 
to decrease appropriately the number of inclusions present at a higher 
temperature. If this effect were negligible it would be possible, as in the 
model, to express the specific heat, C, in the form (l—2)C,+aC, ; but, 
in fact, there will be an important contribution arising from dx/dT. Any 
attempt to calculate the specific heat of liquid helium above the A-point 
would involve a complete statistical description of the assembly, and 
would be at least as difficult as an exact calculation of the specific heat 
arising from an order—-disorder transition. Since, therefore, we cannot 
justify the analogy between the actual liquid and our model so far as the 
specific heat is concerned, we shall modify the argument so as to incorporate 
experimental values of the specific heat. The calculation may be 
simplified without serious error by taking as zero the expansion 
coefficient, f,, of He IT. This assumption is equivalent to neglecting ) in 
comparison with « in the model already analysed, and is justified by the 
fact that the relative change in specific heat at the A-point is much greater 
than the relative change in expansion coefficient. At the same time we 
must for consistency put y,=1 and hence k,=k,,. A short calculation 
now leads to an expression for the real part of the mean compressibility 
at frequencies such that wr<1, which, as may be seen from equation (10), 
may be obtained by assuming thermal equilibrium between matrix and 
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inclusions. The rate of input of heat per unit volume is q,(1—z), since 
9o=Boq,/B,=0, and thus the rate of rise of temperature of the matrix 
is q,(1—a)/C. Since for the pure matrix the corresponding rate would 
be ¢,/C,, we have immediately that h,=C,(1—a)/C. Hence, from 
(1) and (3), 


k[keyg=1+aAk/k,,+(y;—1)(1—2)[1—0(1—2)], 
where Ak=k,,—k,, and b=C,/C. 


If now we write v, for 1/4/(p,k,,), the velocity of a sound wave in an 
extended volume of the pure matrix (He I), and v, for 1/4/(p,k,,), the 
velocity in an extended volume of the material of the inclusions (He IT) 
and further introduce the definitions : 


> 


Av'=v,—2, 
where v is the velocity in the actual material, we have that 


Av'=2Avy+4v,(y,—1)(1—a)1—b(1—a)].  . . . (15) 


It is worth emphasizing that y, in (15) is not the experimentally 
determined value of y in actual Hel, but the value which would be 
observed if there were no inclusions of He II. This “ideal y”’, as we 
shall show later when comparing theory and experiment, can be fairly 
precisely estimated by extrapolation from higher temperatures. 

So far as the attenuation coefficient is concerned, it seems likely that 
the expression (14) derived from the model, for values of wz much less 
than unity, needs little modification when applied to the real liquid ; 
the attenuation arises from the lag of the temperature oscillation of the 
inclusions behind the pressure oscillation, which should not be influenced 
to any great extent by the variations in the mean lifetime of the inclusions 
brought about by the temperature oscillation, provided that the 
amplitude of the temperature oscillation is small. We shall therefore 
use (14) in comparing theory and experiment, choosing some plausible 
value for C,, the specific heat of the material which composes the 
inclusions. It is also necessary to ascribe a value to f, which can only 
be done by guessing the mean life-time of the inclusions. It will be 
remembered that f=1 if this life-time is much greater than 7’ and 
decreases to zero as the life-time becomes much smaller than 7’. Now 
since 7’ is the relaxation time for heat transfer between the matrix and 
the inclusions it is reasonable to suppose that the life-time of the 
inclusions will be of the order of 7’, under which conditions f will be 
approximately 4. Probably at temperatures well above the A-point, 
at which inclusions are only rarely formed as a result of an unusually 
large fluctuation, the life-time of such inclusions will be rather shorter 
than 7’ and f will be correspondingly smaller; conversely, near the 
A-point the inclusions will be readily formed, and may persist for times 
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longer than 7’. In the absence of any detailed theory of the fluctuations 
the best that can be done is to take } as a compromise value of f with 
reasonable assurance that except very near and very far from the A-point 
it will not be wrong by more than a factor 2. We shall therefore write 
for the attenuation coefficient : 

_ w*C,a* Av 

is 9K«v?. 
where, as in the model, v, is the limiting velocity at high frequencies, 
corresponding to a compressibility (1—2x)k,,+%k gq, and Av=v)—v. 


(16) 


b 


§ 4. COMPARISON OF THEORY AND EXPERIMENT. 


In order to compare the predictions of (15) and (16) with experiment 
it is necessary to choose values for the physical quantities entering into 
these expressions. The two which particularly need discussion are C, 
and y, the specific heat and ratio of principal specific heats of 
“ideal Hel” containing no inclusions of HeII. These must be 
estimated by extrapolation from higher temperatures, and, since the 
properties of helium in the range from 2-18° to nearly 3° K. are to varying 
extents influenced by the presence of inclusions, the required extrapolation 
is considerable. Fortunately, however, a few plausible assumptions 
enable the extrapolation to be carried out without too much uncertainty. 
The data which will be used are to be found most conveniently exhibited 
in Keesom’s (1942) book, and the references in the following paragraphs 
are to diagrams therein. 

We consider first the ‘‘ ideal specific heat ” of helium I, C,. It is clear 
from diagram 4.17 that between 4° and 3° K. the specific heat falls 
almost linearly, that is, much more slowly than might be expected for 
a normal liquid, and this trend is exhibited in an even more striking 
fashion by the entropy diagram 4.40. From the latter diagram it is 
clear that the entropy of solid He and HelIlI at 0°K. takes a value 

about —0-85 (the entropy of liquid HeI at the normal boiling point is 
taken as zero), and it is very nearly to this value that the entropy curve 
of HeI above 3°K. extrapolates. It is likely, therefore, that the 
A-transition in helium is in this respect analogous to that in super- 
conductors rather than to order—disorder transitions; that is to say, even 
if the phase-transition did not occur there would be still no excess 
entropy at 0° K.* In parenthesis it may be remarked that if liquid He, 
is similar to He, I, then there need be no cause, on thermodynamical 
grounds, for surprise that no transition occurs. All that is necessary to 
eliminate the entropy by the time the liquid has cooled to 0° K. is that the 
specific heat shall vary roughly as T, as apparently occurs in He I. To 
return to the immediate problem, we shall assume that the entropy 
curve 4.40 can be extrapolated almost linearly to a value —0-85 at 0° K., 
and hence derive an almost linear ideal specific heat curve for He I. 
EE OF Oi A eg a en 


* Note added in proof.—This behavi ; 
and Mendelssohn (194 4 ehaviour has already been noted by Daunt 
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In fig. 1 are exhibited the experimental curve of C and the assumed 
curve for C,, from which the ratio b(=C,/C) may be determined as a 
function of temperature. In this diagram C and C, are referred to unit 
mass rather than unit volume, and are measured at the vapour pressure 
of the liquid. . 

The ideal ratio of the principal specific heats, y,, may be derived from 
the thermodynamic relation : 


atl ee Ah ( Op Op 
eC NCL) NOL): 
where C,, is the specific heat per unit mass at constant pressure. 
We have already considered the extrapolated curve for C,, which is 


measured at the vapour pressure, and a small correction for the variation 
of vapour pressure with temperature must be applied to convert it into 


Fig. 1. 


Specific heat : cals./deg. gm. 
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Temperature; °K, 


Specific heat of HeI: (a) Experimental curve. 
(6) Assumed curve for “ ideal Hel”. 


a curve for C,. The density, p, changes so slowly with temperature that 
there is no need to correct the measured values, but the two differential 
coefficients in the expression are affected by the approach to the A-point 
and require correction. Keesom’s curves 4.38 for (0p/0T), at various 
densities give a clear indication of how the extrapolation should be made, 
since they lie close together and run parallel above 2-7° K., only diverging 
as each approaches its discontinuity at the A-point. Now the A-point 
for helium of density 0-18 gm./cm.? lies 0-3° K. lower than that of helium 
under its vapour pressure, so that at 2-2° K. it will hardly be affected 
by the approaching discontinuity. The ideal curve for (dp/0T), of 
helium under its vapour pressure may therefore be constructed by 
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continuing the curve above 3° K. parallel to the curve corresponding 
to a density of 0-18 gm./em.3. No similar method is available for (dp/0T),, 
which is exhibited in 4.39, but the extrapolation presents no difficulties. 
Values of (p/0T), were plotted against T for temperatures above 3-0° K. 
and the curve was extrapolated by eye to 2-1° K., bearing in mind that 
it must tend to zero, according to Nernst’s theorem, as T tends to zero. 
From these extrapolated curves a curve for (y,;—1) was constructed, and 
the result is shown in fig. 2. 

We are now in a position to see whether the expression (15) provides a 
satisfactory explanation of the variation of the velocity with temperature, 
of which the experimental curve (Atkins and Chase 1951) is shown in 


Fig. 2. 
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Ratio of principal specific heats of “ ideal HeI”. 


fig. 3 (a). One might tackle this problem by guessing a reasonable form 
of the curve for v, (which should join smoothly on to the experimental 
curve at about 3° K.); the difference between this curve and the 
experimental curve would give values of Av’, which could then be 
compared with the predictions of equation (15). We prefer, however 
to work in the reverse order, estimating dv’ from (15) and hence 
constructing a curve for v,. The aim is now to see whether a curve may 
be thus obtained which behaves in the expected way, becoming steadily 
flatter as the temperature is lowered, and tending +6 a zero gradient 
at 0° K. im accordance with Nernst’s theorem. Since the major 
contribution to Av’ arises from b rather than 2, especially at the higher 
temperatures, we may as a first approximation to Jv’ put x=0 in 
equation (15), obtaining the simpler expression 


Av’ =43v,(y,—1)(1—b). 


% 
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Values of Av’ calculated in this way for different temperatures and added 
to the experimental values of v give a first approximation to the curve 
of v, versus T, which is shown as curve (6) in fig. 3. It is encouraging 
that this curve joins smoothly to the experimental curve above 3° K. 
Below 2-5° K., however, it takes a less satisfactory shape, and it is in 
this region that non-zero values of x must be assumed in order to smooth 
_ out the curve. It is now necessary to use values of Av, in (15), and 
therefore to estimate v,, the velocity of sound in the material of the 
inclusions. Since it is not possible without a detailed theory to determine 
the extent to which the inclusions resemble small regions of He IT, there 
must obviously be considerable latitude in the choice of v,. For the 
present comparison a value of 217-5 m./sec. has been chosen, corresponding 
approximately to the velocity just below the A-point, but other values 
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Velocity of sound in helium ; 
(a) Experimental curve. 
(b) First approximation to curve for “ideal Hel”. 
(c) Final curve for “ ideal HeI’. 


are equally likely, and the calculations which follow could readily be 
modified to suit them. It is now a matter only of finding such values 
of x that a smooth curve for v, may be drawn, using the complete 
expression (15) for 4v’. One possible variation of x with T is shown in 
fig. 4, and the corresponding curve for v, as curve (¢) in Ligon Lhe 
values which it is necessary to ascribe to x appear to be fairly reasonable, 
but, since the velocity is not highly sensitive to variations in x and the 
parameters occurring in (15) are not very accurately known, it would be 
unwise to regard the curve for x as having more than qualitative validity. 
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However, we may conclude from this analysis that the hypothesis of a 
non-uniform composition of helium above the A-point is capable of 
accounting for the observed velocity curve, since the form of curve 3 (c): 
is very much what one would expect for a uniform liquid, such as 
‘ideal He I’ is assumed to be. : 
Finally we must see whether equation (16) predicts an attenuation of 
the right order of magnitude. Here the comparison between theory 
and experiment is even more uncertain than with the velocity, since 
depends critically on the size of the inclusions, about which nothing is 
known. We shall therefore calculate the value of a which will lead to the 


Fig. 4. 
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right value of « at a given temperature, and see whether inclusions of 
this radius are likely, according to simple fluctuation theory, to be present 
at that temperature in the concentration, x, found above. It will be 
necessary to choose a value for C,, and here again there is very little 
guidance from theory which will enable a wise choice to be made. 
According to the measurements of Keesom and Keesom (1935) the 
specific heat of He IT in bulk rises, just below the A-point, at least to 
6 cals./gm. deg., but the rise is very rapid, and it is unlikely that the 
inclusions on the average take so high a value. We shall here assume Cy 
to be 1-5 cals./gm. deg., or 0-22 cals./em.3 deg. According to the recent 
data of Atkins and Chase the attenuation constant «, at 2:3°K. anda 
frequency of 14 Mce./s. is 0-27/em., of which approximately 0-10/em. is 
explicable in terms of viscous damping and thermal conduction in 
“ideal He I” (Pellam and Squire 1947). The difference, 0-17/em. we 
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suppose to result from the presence of the inclusions. Ruane 
Y= 224-3 m./sec., Av=4-5 m/sec. and K=4:4x.10-5 cals./deg. em. gec.* 
in equation (16), we find that the observed attenuation may be 
accounted for if the mean radius of the inclusions is about 2:1 10-7 em., 
corresponding to a cluster of about 850 atoms. 

In order to determine whether inclusions of radius 2-1 10-7 cm. 
might be present in such numbers as to occupy 9 per cent of the total 
- volume, as required by our analysis of the velocity curve, we may carry 
out a simplé calculation which, though based on quite inexact 
assumptions, should give a qualitatively correct answer. Suppose the 
liquid at 2-3° K. to be divided into regions whose volumes, V, are equal 
to the estimated volume of an inclusion, 3-9 x 10-2° cm3. _ If these region's, 
having the properties of ideal He I, are regarded as weakly interacting 
with one another we may immediately calculate the proportion which 
through fluctuations have temperatures less than 2-18°K. This 
proportion, by standard fluctuation theory, is given by }{1— erf (I.5T)] 
where 5T=0-12°K. and [=4/(C,V/2kT?) deg.-1, and substitution of 
numerical values leads to a value 11 per cent. The agreement with the 
desired value 9 per cent is probably better than the assumptions warrant, 
especially in view of the uncertainties in the values to be ascribed to 
various parameters occurring in the theory. We may, however, conclude 
that the suggested mechanism of the extra attenuation near the A-point 
is capable in principle of explaining the magnitude of the effect. 

As the A-point is approached the attenuation rises very rapidly, and 
probably increases by a factor greater than 20. This cannot be explained 
on the assumption that the inclusions stay constant in volume, but, on 
the other hand, if they increase in size nearer the A-point, as it is very 
plausible to assume, there is no difficulty in explaining the rapid rise, 
since « depends on a?. Below the A-point the attenuation drops again 
sharply, presumably as the liquid becomes more and more completeiy 
homogeneous, but we shall not attempt to calculate the attenuation when 
the liquid is mainly He II, with inclusions of He I, since in the absence 
of any adequate microscopic theory of He IT the foundations of such a 
calculation would be even less firmly secure than those of the present work. 
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SUMMARY. 


The static field of an infinitesimal loop of dislocation lying in its glide 
plane is calculated. The field of any plane loop lying in its glide plane 
may be obtained by integration of the result. The field produced by the 
sudden creation of an infinitesimal loop is calculated as a function of 
time. The field of an arbitrarily moving plane loop may be obtained by 
integration. The method is illustrated by simple examples. 


§1. INTRODUCTION. 


TuHE process of plastic deformation in a crystal may be represented 
(cf. Cottrell 1949) by the following deformation of an isotropic elastic 
continuum. A cut is made over a finite area of a plane, the glide plane, 
and the opposite faces of the cut are moved relative to one another a 
distance b in a direction, the glide direction, which lies in the glide plane. 
The elastic deformation is thus a particular type of dislocation, as discussed 
by Volterra (1907). Burgers (1939) has shown that the static displacements . 
produced by such a dislocation can be expressed by integrals along the 
line bounding the area which has slipped. We shall here express the same 
general result as a double integral over the area which has slipped, and 
extend the method to give the instantaneous displacements produced by 
a moving dislocation. In §2 we derive the static displacements produced 
by an infinitesimal loop of dislocation in an infinite body ; -§3 digresses 
to give the static displacements produced by such a loop at the centre 
of a finite sphere. In §4 the method of synthesis is illustrated by the 
screw dislocation, in which the cut surface is a half plane bounded by a 
line parallel to b, and the edge dislocation, in which the cut is made 
over a half plane bounded by a line perpendicular to b. The displacements 
produced if an infinitesimal loop of dislocation is suddenly created are 
calculated in §5, and, as an example of the method, the displacements 
produced by a rigid motion of a screw dislocation are calculated in § 6. 


§2. THe Frevp oF an InFrntresmmaL Loop oF DistocatTion Lytna IN 
Irs GLIDE PLANE. 
Consider an infinitesimal loop of dislocation with Burgers vector { 
b=(0, 0, —b) lying in the az plane close to the origin. Let the area of 
+ Communicated by the Author. 


d The vector b used in this paper is identical with the vector f of Burgers 
(1939). The current convention (cf. Frank 1951) is to take b= —f. 
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the loop be oA, while a typical element of its perimeter is dé. The stress 
field of this dislocation in an infinite body may be found from the formula 
of Burgers (1939), which gives at the point x a displacement u of the form 


ae 4k 
ere. Bi eA DLAC Vis ral ae va L's |. aE (1) 
Geb pee ata eke ee. ek (2) 

b - 
uta2al SF “) 
and 

A+ ; d 

b= tg | bate) 2. Ns a ee 


Here Q is the solid angle subtended by the loop at the point x, A and 
are Lamé’s elastic constants, and r=|x—§]. 
In this case, 


47u*=(0, 0, —b8Ay/r3), | 
47u**=(0, —bdAz/r3, 0) 
and 4p (A+p)bd A ye/(A+ 2u)r3, 
; A+p 3xyz 
lead Sy eae ee 
eading to Bau, See ee \ (5) 
A+p (3y22 z z 
Sal ae —_ = 
3 re ae re =) r? 
Ae (SY YN CY 
nd oem) Pa Raw ae 
im hae =e ( f° Z)+ $ | 
where 
B=:—b8A/4z. (6) 


This field may be derived from the “ first type of simple solutions ”’ 
of Boussinesq, given by Love (1927) at §§131, 132. In Love’s notation 
it is 

(Uz, Uy Uz) = 30 (Us, Va; We)/0z+ 50(Us, U3; ws)/dY. GP tees, © (7) 

If w;, is the component in the direction of x, of the displacement at the 
point x, produced by a force in the direction x; at the origin, this 
displacement field is of the form 

Ou,;; . Ou; 
—s ud pant , ‘7 
Unis (ei + mt) Geely eel 1 teers oS) 

The forces applied at the origin have zero resultant and zero moment. 

The corresponding diagonal terms with i=k are “ double forces without 
moment ’’ (Love). The strength of an infinitesimal loop of dislocation 
is thus a tensor, in which one suffix represents the normal to the plane of 
the loop and one the Burgers vector. The non-diagonal terms represent 
elementary acts of slip, while the diagonal terms represent non- 
conservative motions (Nabarro 1951). 
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§3. A Loop aT THE CENTRE OF A SPHERE. 


If the dislocation exists in a finite body with free surfaces, a comple- 
mentary field must be added which neutralizes the traction across the 
boundary. For a sphere of radius R centred on the loop, such a field is 


given by 
yB Mu, = 48(A+p)(2A-+ Tp )ey2, 
yB-M,,=48(A+p)(2A-+ Te) y?2— 24(A+ pw) (5A Tp) rez 


4.5(27A2-+ 56Au-+ 282) R22 . (9) 
and = yB-lw,=48(A+p)(2A-+ Tu)y2?— 24(A+p)(5A+ Th)ry 
+5(27A2+ 56Ay+ 287) Ry, 
where 
y= 2(A+-2u)(19A-- 14p) Ro) er 


These displacements may be derived by a superposition of two solutions 
of the types given by Love at §172, taking w.=yz and ¢,.=yz. 
On the surface r=R, the total displacement derived from (5) and (9) 
is given by 
yB MU, =210(A+p)(A+ 2p )ayz, 
yB-%u,=210(A-+p) (A+ 2n)y22+ 15A(A-+ 2) R22, pr eh 
and yB-1u,=210(A+p)(A+ 2u)yz?+ 15A(A+ 2p) Ry. | 
‘In spherical polar coordinates with Ox as axis, these become 
yB-lu,=15(8A+7)(A+ 2) R8 sin? 6 sin 24; 
yB-lu,=15X(A-+ 2n)R sin 6 cos 6 sin 246 oh kn cart 
and *  yB-luw,=15X(A+ 2u) R38 sin 6 cos 2¢. 
Suppose a homogeneous stress p;,; is exerted on the sphere before the 
dislocation is made. The work done by the displacements (12) against 
this stress is obtained by integrating w,p,.,-- upp ro UgPr, Over the surface 


of the sphere. The only components which do not vanish in the integration 


over ¢ are those arising from p,, and p,,, which are 


Prr=Pyz Sin? 8 sin 2¢, | 


* Pro=Py, 8in 0 cos 6 sin 26 (13) 
| 
and Pro=Pyz Sin 8 cos 2¢4. | 
The work done W is given by 
yBW=8(A+2y)(19A+ 14y)7R5p,, 
14 
or W=—bdAp,,,. oe 


Synthesis of Elastic Dislocation Fields 1227 


§4. SYNTHESIS OF THE EDGE AND THE Screw. 


Returning now to the infinite body, we may take the loop to be 
rectangular, with 5d=dx6z. By replacing z by z—¢ and integrating (5) 
over all ¢, we obtain the displacements produced by a pair of screw 
dislocations of opposite sign lying close to the z axis, and separated by da. 
These are 

uU,=0, Uy=9, Uz= —byda/2n(w?+y?), 2 . . (15) 


Replacing x by x—é and integrating again from ¢=— co to £=0, we 
obtain the displacements produced by a single screw dislocation along 
the z axis, which, apart from a constant, are 

Ub ga= OF Uy =, Uz=—(b/2r) tan-ty/e. . . . (16) 


Alternatively, we may first integrate over all € to obtain a pair of edge 
dislocations close to the x axis, with displacements 


U.==0, 

2 Sa ee ) 

vA +2) PF2P 2r(AF2n) FPA? b . . (17) 
bdz(A+p) yz? bdz y 


Ne aA 2) (YAP Br 2u) PE 


To obtain an edge dislocation lying along the x axis, we integrate from 


f=—aotol=0. The resulting displacements are, apart from a constant, 
U,=0, 
SY eS 0) Cee ee ares 
Uy P(A Du) eet arpa) oY TF | as) 
b(A+h) YZ 


b y 
~ tan-12, 
FALSE) page te 2 


$5. THE SUDDEN CREATION OF AN INFINITESIMAL Loop. 


The movement of a dislocation in its glide plane may be regarded as 
the result of the continuous addition of infinitesimal loops to its boundary. 
If the displacement produced by the sudden creation of such a loop is 
known as a function of time, the field of a dislocation moving arbitrarily 
in its glide plane with subsonic speeds may be obtained by superposition. 
The known solutions for the field of a moving dislocation have each been 
obtained by some special artifice: the present method gives an explicit 
solution in the general case, although, as will be seen from §6, ‘its 
application in simple cases is rather laborious. ; 

The displacement arising from the sudden creation of a loop of 
dislocation may be derived from that arising from the sudden application 
of a single force by the same process of differentiation and combination 
(7) as can be used to derive (5) from the known expressions for the static 


402 


{228 F. RB. N. Nabarro on the 


displacements produced by a single force. Love (1904) has given the 
displacements arising from a single force suddenly applied and then 
continuously maintained. If the force is applied instantaneously, the 
results involve discontinuous functions, and the differentiation of these 
is not straightforward. It is more convenient to start from the formule 
of Stokes (1849) giving the displacement produced by a force which is 
an arbitrary function of the time, and to carry out the differentiations 
before passing to the special case of a force suddenly applied. A typical 
formula is given by Love (1927, §212). The displacement produced at 
time ¢ at the point x, by a force in the x, direction applied at the origin 
and proportional to y(t) has a component w,; in the direction of x; given by 
2,—-1 rie 
Uiy= a | t'y(t—t’) dt’ +1 Fax,[a-?x(t—r/a)—€ x (t—1/€)] 
+38.¢7c*y(t—4/¢), | ee 
where a and c are the speeds of dilatational waves and of shear waves. 
The displacements corresponding to the formation of a loop of dislocation 
in the orientation already considered are given, in analogy with (7), by 


ria 


U j= Oy; | OZ OU OY. "ea en ce 
Substituting (19) in (20) leads to the expression 
o3r-1 r/e 
U;=2 Ox ,0yOz iB t'y(t—t’) dt’ 
=~ Se (12-2 x) [a-*y(t—r/a)—c-*y(t—r 
ie dr X r/c)| 


1 d(yz)[, d 
- iors | 2a-tx(t—r]a)—30-8x(t—rle) +r Tetylt—e) | : . (21) 
If the force is applied suddenly at time t=0 and then maintained, 
x(t)=0 for i<0, x(t)=F for t>0. The expression (21) then takes different 
forms in the three regions t<<r/a, r/a<t<r/c, t>r/e. It may be written 


t<+r/a: F-1u,=0. © ey oi es Ae Re en 
tlast<rje: F-t4,= ae {= _ 15t—27fa-85(1—r/a)—e-*5(t—r[)]} 
+3 Ge [MH Bore 80—r)| Sain, they 

and 
t>r/¢: Fu,=— our (5 - =) _ a 2 : Bee eae uit tb 


In (23), terms have been omitted which occur only at the end points 
of the range, and are there finite, so that they contribute nothing to the 
integral over a distribution of loops. If the density of the medium is 
we have i ey 

Ge=(A4- 20) pps ne Dag 
and (24) agrees with (5) if 


F=—@B=beAj4r, . . . . . . . (25) 
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§6. Rigtp Motion or a Screw DisLocaTion IN A STRAIGHT LINE. 
Consider a straight screw dislocation lying parallel to the z axis, and 
moving so that at time 7 it intersects the plane z—0 at the point (€(z), 0, 0), 
where |£’(7)|<c. The motion of the dislocation is equivalent to the: 
creation of rows of infinitesimal loops of the type already discussed. 
The displacement corresponding to the creation of a single loop spreads 
out in a wave of finite thickness, bounded by the waves of dilatation and 
of shear. The displacement at (x, y, 0, t) corresponding to the creation of 
a single loop at (€, 0, ¢, 7) begins as soon as t—r=r/a, where 
ese y eyes ette i Oe Ve the (26) 
The component w, arising from the point ¢=0 does not vanish, and arrives 
when t—7=p/a, where 
(el Se Ee ae, Inte aa ee eee nme 76 
According to (15), the resultant displacement obtained by integrating 
over all ¢ is without dilatation in the final state. It does not follow from 
this alone that the disturbance spreads from the line of the dislocation 
with the speed c of a shear wave. For example, (Love 1927, § 213, (iv.)) 
the displacement around a centre of compression, which in the static 
case (ibid. §132, (a’)) is a pure shear, is propagated entirely with the 
velocity a of longitudinal waves. However, in the present case the 
resultant displacement is at all times parallel to the z direction, and 
independent of z. The dilatation always vanishes, and the disturbance 
spreads as a shear wave, being first observed when t—t=p/c. Disturbances 
from different parts of the dislocation line propagated with speeds lying 
between a and ¢ arrive before this time, but continuously cancel one 
another. . 
- In synthesizing the displacement, there are three situations to consider. 
t—7<p/a: the longitudinal wave from the nearest element of the 
. dislocation has not yet reached (x, y, 0) ; 
plaxt—r<p/c: the longitudinal wave from the nearest element has 
arrived, but the shear wave has not ; 
i—7t>p/c: the displacement produced by the nearest element has 
reached its final value. 
The displacements u, and u,, are odd in Z, and vanish on integration. 
. In the case t—7<p/a, the expression t—7—r/a is negative for every 
element of the dislocation. - It follows from (22) that wu, vanishes. 
In the case p/ja<t—r<p/c, the expression t—7—7/a is negative when 
(2>Z and zero or positive when (?<Z7, where 
eee ZO (bet) aus ass et, erase (29) 
The expression —7—7/c is always negative. The contribution of elements 
with ¢2>Z? is zero, and the contribution of the remaining elements is 
given by (23) and (25) as 


dou Za ((3 150 BOS] 
eo ( Go) tae a 


beysé | zy 
= 2¢%r-¥a-*3(t—7—r/a)| Dea -(29) 
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where the integral is taken over the closed range. The terms in 6(t—7—7/c) 
have been omitted, since they vanish in this range. The terms not 
involving 8(t—r—r/a) are readily integrated, giving 4Z _/a>(t—7)®. The 
contribution of the § function from the points ¢?=Z? exactly cancels this 
at all times in the range. There is thus no disturbance until the shear 
Wave alrives. 

In the case t—r>p/c, there is again no contribution from the points 
with 2>Z?. When Z?>¢2>Z?, the elements contribute according to 
(23), and when £2<Z? they contribute according to (24), where 

t+ Ze c8t—7\8.0) OF ee 
The contribution from the second region is given by 
4ru, ‘ZEN See LOLs eects 12¢2 2 
bcp 2, 1G a or lta) le) a8 Tes 
— 20?r-4[a-8 (t— r—r/a)—c 38(t—7—1€)] 


—r-te-85(¢—7—rfo)} Th hat aan 2 SS 


taken over the closed range. The terms not involving 6 functions yield 
4, 4Z,, 1. Pome 
The 6 functions yield 
a 4Z, 


and the total contribution from this region is 


bo 


UP LE 1 i| 2 
beyse  c%K(t—7)3 (= a, 3) it nea hee (32) 


The contribution from the region ¢? <Z? is given by 


47ru, fe 362/41 1 1 
botySE 5 vu) cen eb (5 e x) a saat gs 


223 22. 
=— pcr(t—7) a @c3(t—7)3 : ot hye cae OG Roel 
Adding (32) and (33) gives the total displacement in the region t—7>p/c 
as 


_ by®é c(t—T) 
Ast increases, this tends to the static value given by (15). 
In this particular problem the result could be obtained much more 
easily by noting that the displacement is ultimately equal to (15), and 
travels as a pure shear wave. This suggests a trial solution of the form 


w= 4s (4). 1. Gs aie ten ses 


(34) 
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The solution must obey the equation 


(Ch 0 Wane : 
Cee ee arc. (36) 
which is satisfied if ; 
eae) falth)=<(aue—e 1) f (ab) vee eile (ST) 
where 
w= p/ct. sgl ere YO coke 4) coalesce 3S) 


The solution of (37) subject to the boundary conditions at w=0 and 
w= © leads at once to (34). 
The total displacement is 


by (7 E'(r)c(t— 7) dr 
eae = 2c | o Pct—T)°— py” (39) 
where 7, is given by 
[g—f(79) Py et—t)P (40) 


For the particular case of steady motion, é(r)=v7, and (39) may be 
integrated by introducing the variable 7, where 


COS 7=COsS 7)/COS «, 


sin = p/c(t—7) 


and sin «=vy/c| (x—vt)?+y?]*. 
This leads to the known result 
b (1—v?/c?) ty 
es ee -1 A 4] 
We Qa — x—vt ( ) : 


_ I am indebted to Dr. J. D. Eshelby for suggesting ge in 
the presentation of this paper. 
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SUMMARY. 


The scattering constant has been determined by measurements of the 
multiple scattering along proton tracks ending in the emulsion. A value 
of 27-8-0:5 has been obtained (for 62~0-03 and an average cell-size of 
72 p). 

A simple method involving the combined measurements of multiple 
scattering and residual range has been used to identify the particles. 
A parameter P= 2-37 log &,-+ 1-37 log R was determined for each track 
(R being the range in microns and &,, the mean angle of multiple 
scattering—per 100 .—measured along the faster half of the trajectory) ; 
P is a linear function of the logarithm of the mass of the particle. 
A scale of values of P has been calibrated experimentally using tracks 
of protons from artificial sources and those of o-mesons. Using this 
method tracks of 83 protons have been identified amongst a group of 
136 particles of charge |e| ejected from cosmic ray stars. The mean 
value of P for the 83 protons was found to be 4:275-+0-013. The resolution 
_ between z-mesons, protons, deuterons and tritons was found to be 
satisfactory. An appendix gives the coupling matrix calculated for the 
sagitta method, taking into account the loss of energy. 


$1. INTRODUCTION. 


THE scattering constant as defined in a previous paper (Gottstein et al. 
1951) can be determined by measuring the multiple scattering along 
tracks which end in the emulsion, and which are produced by particles 
of known mass. ; 

The main advantage of this method is that the value of the initial 
energy of the particle can be determined reasonably well by measuring 
its residual range. The inaccuracy inherent in such a range measurement 
is of the order of 2 per cent—this being the probable error due to 
straggling (Rotblat 1950). The calculations are more complicated 
because, in contrast to the normal case in high energy tracks, the energy 
here can no longer be considered constant in that part of the track used 
for scattering measurements. The loss of energy along the track by 
1onization must be taken into account. 
ee 


* Communicated by Professor C. F. Powell, F.R.S. 
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In a previous paper, Goldschmidt-Clermont et al. (1948) have used a 
statistical variable «; defined as the product of each angle a, by the 
energy T; corresponding to the centre of the cell. The mean Je| is a 
parameter which depends only on the mass of the particle and the 
scattering constant. In principle, this method enables one to carry on 
the measurements to the very end of the track—decreasing the cell size 
as the end is approached, and so increasing the statistical weight of the 
result. But such a procedure would make it difficult to determine the 
scattering constant, as energy is lost rapidly towards the end of the track 
and it is hard to eliminate effectively the Rutherford single scattering 
from the measurements. 

However, if the scattering be measured in a part of the trajectory 
where the change of energy is not too large, then a constant cell length 
could be maintained for all the measurements, thereby introducing 
a simplification into the determination of the parameter e. 

In the present investigation, we have arbitrarily limited the scattering 
measurements to the faster half of the track. In our calculations, the 


1 
average 7a le ;I';| has been replaced’ by the statistically equivalent 


1 1 1/2 
quantity (52 | xl) > é aT; s) . The identical mathematical expecta- 


tions and the slightly differing variances of these two parameters are 
indicated in the Appendix. 

For the present investigation, we have employed 83 tracks of slow 
protons which ended in the emulsion. The experimental procedure 
adopted for identifying the particles producing these tracks is described 
below. 

§ 2. EXPERIMENTAL PROCEDURE. 


1. Selection of Tracks for Measurement. 


One hundred and fifty-four tracks ending in the emulsion, of length 
greater than 500, originating from cosmic ray stars, have been 
systematically investigated. Tracks due to particles of charge greater 
than one (i. e. He, Li, Be, . . . nuclei) were identified by visual examina- 
tion of the track for delta rays (Sorensen 1949, 1951) and by gap 
measurements (Hodgson 1950). Highteen such tracks were found amongst 
the 154 tracks investigated. The remaining 136 tracks judged to be 
caused by particles of charge one, were used for the range and scattering 
"measurements. 

. 2. Measurement of Range. 


The tracks were divided into segments of 48 » and smaller. Horizontal 
and vertical projections were measured in each segment. In calculating 
the vertical projection in the unprocessed emulsion a shrinkage factor of 
2-5 was adopted. From repeated measurements on the same track the 
total probable error in these determinations has been found to be less 
than 1 per cent. Tracks inclined at angles greater than 15° to the surface 
of the emulsion were not used. 
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3. Measurement of Multiple Scattering. 

The first half of the range determined above was used to measure the 
multiple scattering. This distance was divided into “n” cells, each of 
size “s”. ‘‘s” will be referred to as the “ primary cell size”. The 
tracks actually measured varied in length from 500, to about 4000p, so 
the numbers “n” and “s” have also varied. ‘“s”’ was so chosen that 
“nn? was at least equal to 30, the values of “s” being limited to 
sequence 8, 16, 24, .. . 48 (8 being the smallest unit on our eye-piece 
scale). In one eye-piece we used a scale of 48, divided into segments of 
8 each, placed parallel to the displacement of the stage. The track 
was moved in steps of sy along this scale. In the other eye-piece and 
perpendicular to this displacement, we used the scattering scale as 
described in previous papers. The position of the track was read on this 
scale at each step. 

4. Calculation of the Mean Angle «. 


Second differences were computed in cell-sizes which were integral 
multiples of the primary cell-size. This enabled us to choose the optimum 
cell length, defined as that for which the d has a value greater than four 
times the total noise in the cell. The noise of the microscope M4074 was 
determined for different cell sizes by use of a primary track 2 cm. in 
length, which had a finally determined mean angle of scattering of 
0-0007°/100,. when measured in 4000, cells. The results of these 
measurements have been described in a previous paper (Gottstein et al. 
1951). Noise elimination as described in previous contributions (Menon 
et al. 1951) cannot be used in this case owing to the simplification we 
have introduced. A signal to noise ratio of at least four must be 
maintained to keep the influence of noise on the final result to within 
a few per cent. The absolute values of the second difference thus 
obtained were transformed into an angle & expressed in degrees per 100 x. 


$3. IDENTIFICATION OF Mass. 


Using this angle and the value of the range already measured, we have 
determined a parameter P=1-37 log R+2-37 log & (R in p and & in 
degrees/100 1). This parameter is a linear function of the logarithm of 
the mass of the particle. It has been computed from the range energy 
relation expressed in the form of a power law R=E!*3xM and from the 
definition of % (see Appendix). 

Before P can be employed as a parameter to identify masses of the - 
particles, a scale of values of P must be calibrated for particles of known 
mass. This calibration has been carried out in the following manner. 

We chose twenty-four tracks of o-mesons which ended in the emulsion and 
which were longer than 5004. They have all been attributed to 7--megons. 
The possibility that these o-mesons are due to --mesons is extremely 
small, for these reasons: the tracks were chosen from plates exposed 
in balloon flights at 65,000 ft. and 75 per cent of the o-meson tracks 
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chosen had produced stars with more than two prongs. It is known that 
o-stars produced by »~-mesons are mostly 1-pronged (George and Evans 
1951). The probability of their producing stars with two or more prongs 
is very small. It is therefore reasonable to assume that almost all the 
o-meson tracks are in fact due to 7-mesons. 

By measuring the range and the scattering in the first half of the 
trajectory, one may compute a value of P for each of these tracks, and 
plot them in a histogram as in fig. 1(a). These values are distributed 


Fig. 1. 


Histogram showing the values of the parameter P=1-37 log,) R-+-2-37 logy 
&(& in °/100u, R=total range in microns), for : 
A—24 tracks of o-mesons. 
B—28 tracks of protons from the N1 (d, p) N?° reaction. 
C—30 tracks identified previously by grain density versus multiple scattering. 
D—136 tracks due to particles of charge 1 ejected from cosmic ray stars. 


The ordinates give the number of tracks per interval of 0-05 of P. 

The arrows, P, D, T, refer to the calculated values of the parameter P for 
protons, deuterons and tritons, assuming P;—5-1 and mass of proton/mass 
of 7-meson=6-7. 


with mean 5:1. Adopting P=5-1, for 7-mesons, and a value 6-7 for the 
ratio of the mass of the proton to that of the 7-meson, values of P for 
different particles may be calculated as:  P(proton)=4.27, 
P(deuteron)= 3-98, P(triton)=3:80. 


\ 
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The values of P(proton) were calibrated experimentally using tracks 
of protons (1200 »>range> 1000 which came to the end of their range 
in the emulsion) from the reaction N'4(d, p)N¥. The Q value of this 
reaction is 8-165 MeV. ‘These were obtained in plates exposed to 
8 MeV. deuterons in the Liverpool cyclotron. Fig. 1 (B) shows the P values 
for 28 of these tracks. 

To study the overlap of the values P(proton), P(deuteron), P(triton), 
in a mixed group of particles, the calibration has been extended to include 
all particles of charge 1, using tracks ending in emulsion of length 
+3 mm. originating from cosmic ray stars. The particles producing these 
tracks had been identified in this laboratory, by measurements of the 
grain density and multiple scattering by Camerini, Fowler, Lock and 
Muirhead, who have very kindly given us their results. The P values 
obtained for this group of 30 tracks are shown as a histogram in fig. 1 (Cc). 

We then determined the P values of 136 tracks (previously identified 
as being due to particles of charge 1). These were from cosmic ray stars 
in Ilford G5 plates, 400 thick, exposed in a high altitude balloon flight. 
The P values obtained have been plotted as a histogram in fig. 1 (D). 
The peak for P(protons) is in agreement with those obtained in fig. 1 (B) 
and (c). It is also consistent with the value of 4-27 obtained using 
Pr=5-1 and m,/m,—6-7. Corson (1951) has made a plot of the 
parameter p=(a(x).R°°%/,/x),,, where x and R (the cell-size and 
range respectively) are expressed in units of 1004. He obtains a value 
of p~5-9 for 76 proton tracks. Our results similarly plotted, give 
p~56. The agreement may in fact be even better than these figures 
show, since Corson remarks that since he uses a cell of 20 for a range of 
500 his values are perhaps a little large, owing to noise introduced into 
the measurements. 

To determine the scattering constant, we have taken the group lying 
between 4:05 <P <4-60 as protons. A certain number of protons which 
exhibit low multiple scattering would have been lost in the deuteron and 
triton peaks. On the other hand, the proton group will contain a certain 
number of deuterons and tritons as impurity. The magnitude of such 
contamination is less than 10 per cent, and its effect on the determination 
of the scattering constant will be very small. 


§4. DETERMINATION OF THE SCATTERING CONSTANT. 


Using the identified group of slow protons, a plot has been made of 
log « against log R. This is shown in fig. 2. The best line fitted to these 
points gives a slope of —0-554+0-011 (the error being the probable 
error). This value is in good agreement with that obtained (0-578) for 
the range-energy relation by more accurate methods. However, as the 
linearity of the regression line is not certain in the region considered, 
the scattering constant has been finally determined as shown in Table I. 

The points on the log «-log R plot (fig. 2) have been divided into groups 
using equal intervals of 0-1 along the log R axis. The coordinates of the 
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centres of gravity of these groups are given in columns 1 and 2 of Table I. 
Column 3 gives the effective range (=-72 R—see Appendix) and 
column 4 the corresponding energy as obtained from the range energy 
relation (Rotblat 1951, Catala and Gibson 1951, Bradner et al. 1950). 


Fig. 2. 


27 29 31 32 3 3-4 35 3: 37 


=!) 
N 


= 
o 


LOG R 


Plot of logy) & oon versus log,) R for tracks having values of the parameter P 
in the interval 4-:05<P<4-60. The best line fitted to these points hasa 
slope—0-554-40-011. 


TABLE I. 


log R | log & | Reg=0-72 R| T(Re) in MeV. | T(Rei) x & 
2-758 | 0-215 8-50 
2-857 | 0-166 9:75 
2-951 | 0-093 11-00 
3-053 | 0-056 


3-152 | 1-965 
3-235 | 1-944 
3-333 | 1-902 
3-412 | 1-845 
3-579 | 1-765 
3-665 | 1-658 
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The product %T(Reg) (which is equal to half the scattering constant in 
the units described previously) is given in column 5—the mean value of 
which is 13-9-.0-25. The final value of the scattering constant as thus 
obtained is 27-8-L0-5 (the errors given being the standard deviation). 
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APPENDIX. 


APPROXIMATION INVOLVED IN CALCULATING |<|j00,- 


The n-dimensional distribution function of successive angles measured 
on the track between the points of residual ranges R, and R,, is 
characterized by a square matrix (A,,). Its elements A;, are equal to 
the quantities «,«,—this being the expectation of the product of the 
ith angle by the kth. 

These elements may be calculated using the range energy relation 
T=K,R” in the case of the coordinate method, following Moliére’s 
20 oT 
l+o / 
(o being the kinetic energy in terms of the rest mass) we obtain 

(see equation, p. 1239). 
In the present experimental investigation none of these terms exceeds 
0-01. For the following argument they may be neglected. 

The experimental procedure as described consists in calculating the 
arithmetic mean of “n” successive angles. This quantity being taken 
as proportional to the r.m.s. value, we may set 


nx 1 l $ { 
lal=521eh~/ 2a). os he me eee Eee (2) 


By definition, the expectation of the last term in (2) is given by 


JGzeaet)= seve (3) fhe. ab 9). @ 


If we replace, in the square root, the Y by the integral we get finally 
ae ee Bn dK a 
le =oK- V8 els ere aE Pu ate oath) 


LK 
®em= 3K Res”, SViviaee Ps Gc tees mae) 


theory (1951). Making the non-relativistic approximation T 


Introducing 


| 0 00.7 4K, /T 
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(2) AU ='@ AjaaTyo0dser oe 23 pue ’@ jo sonyea ouL 
ao tae) ( )te >> pur (peu g ATOAT 


‘uoryRjer AS1oua oSue oy Jo queysuoo =" 57 
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The variance given in (8a) would be greatly increased in comparison 
with that in (8 b) if ioe portions of the track were used. 
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By U. Camerini*, J. H. Daviss, P. H. Fowter, C. FRANZINETTI T, 
H. Murrueap ft, W. O. Lock, D. H. Perkins and G. YEKUTIELI §, 


H. H. Wills Physical Laboratory, University of Bristol ||. 
[Received July 30, 1951.] 
[Plates XLIT.—XLIV.] 


SUMMARY. 


An analysis has been made of the secondary particles ejected from 
nuclear disintegrations observed in electron sensitive emulsions exposed to 
the cosmic radiation at 68,000 ft. Scattering and grain density measure- 
ments have been carried out on the tracks of 2000 particles associated 
with these stars. In addition, the grain density and angular distribution 
of 3070 shower particles and 1508 “ grey” tracks have been measured. 

For 200 stars, the energy of the primary particle which produced the 
disintegration was measured. A detailed analysis was made of such 
events. Single fast 7-mesons of kinetic energy less than 1 BeV. are found 
to interact strongly with nuclear matter. 

An estimate of the frequency of occurrence of neutral mesons was made 
from a consideration of the energy balance in stars of low multiplicity, n,. 


§1. INTRODUCTION. 


In Parts I. to V. of the present series of papers a description has been 
given of the disintegrations produced in photographic emulsions by the 
high energy protons, neutrons and «-particles of the cosmic radiation. 
Parts I. and IT. contained a phenomenological description of the principal 
features of the disintegrations observed at an altitude of 11,000 ft., and 
above 50,000 ft.; in Parts III. and IV. it was shown that a large 
proportion of the shower particles were 7-mesons, which interact strongly 
with nuclei, and their distribution in energy was determined ; and in 
Part V. evidence was given that neutral 7-mesons are emitted, together 
with the charged z-particles, in the “ showers”’, that they have a life- 
time less than 10-13 secs., and that their decay into two quanta gives rise 
to most of the soft component of the cosmic radiation. 

In this and the final paper which follows, Part VII., an attempt has been 
made to elucidate the nature of the physical processes which occur within 
the nucleus under bombardment by energetic protons, processes which 
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result in the production of the “ showers’ of mesons, and nuclear 
fragments. This first paper contains a description of observations on 
15,000 nuclear disintegrations recorded in plates ‘‘ exposed ”’ at 68,000 ft. 
In about 200 stars, the observed track of the primary particle which 
produced the disintegration was long enough to allow the energy to be 
determined. For such events a much more detailed analysis has been 
possible. 

Part VII. contains an analysis of the new experimental material. In 
approaching the problem, it is assumed that the disintegrations can be 
described in terms of a succession of two-body interactions occurring in 
the nucleus (nucleon-nucleon and meson-nuclean) which can lead, if the 
interacting particles are sufficiently energetic, to the creation of a-particles. 
Recent observations, in this and other laboratories (O’Ceallaigh 1951, 
Fowler e¢ al. 1951) strongly suggest that in the nuclear interaction of 
protons and neutrons of energy > 10 BeV. the production of more massive 
mesons may become important. But at the rather lower energies, 
<10 BeV., with which we have been mainly concerned, any such process, 
and those involving the creation of pairs of nucleons, if they occur, can 
hardly play an important role. 


§ 2. NOMENCLATURE. 


In describing the observations, it will be convenient to adopt the 
nomenclatures given in Part I. A closer analysis has been found desirable, 
however, in the case of the more heavily ionizing particles. The tracks 
produced by the particles emerging from nuclear disintegrations are 
therefore divided into three classes :-— 

(i.) Tracks of “ shower ”’ particles. These are tracks in which the 
grain-density, g, is less than 16 grains per 50u ; in the conditions of the 
experiment, this value corresponds to a particle with the electronic charge 
and a specific ionization less than 1-4 times the minimum value, 7 min. 
The great majority of such tracks are produced either by 7-mesons of 
energy greater than 80 MeV., or by protons of energy greater than 
500 MeV. The total number of such tracks associated with a given “ star” 
is denoted by n,. 

(ii.) “Grey” tracks. These are tracks with a grain density in the 
interval from 16 to 80 grains per 50 4, and their number for any star is 
represented by N,. They can be produced by protons with energy in 
the interval from 25 to 500 MeV. Only a very small proportion of the 
protons emitted during the “ evaporation ”’ of highly excited nuclei have 
energies as great as 25 MeV. Most of the “ grey ” tracks are therefore 
due to protons, recoiling from fast nucleons and produced at an early 
stage of the disintegration. Some, however, are due to deuterons and 
tritons, and about 18 per cent to 7-mesons of energy less than 80 MeV. 

(iii.) “ Black” tracks. These tracks have a grain density greater than 
80 grains per 50 «4, and their number is denoted by N,. They are due to 
protons of energy less than 25 MeV. : to deuterons and tritons of energy 
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less than 50 MeV. and 75 MeV., respectively ; and to a-particles less than 
800 MeV. aaa 

The sum N,+N, is represented by N,. A “star” is characterized 
in terms of N,+-n, together. with a suffix to denote the nature of the 
particle producing the disintegration ; viz., n, p, 7, «—for neuton, proton, 
m-meson or x-particle respectively. Stars with n,>2 are referred to as 
“ showers ”’. 

§ 3. SCOPE OF THE OBSERVATIONS. 

In an examination of 87 ¢.c. of emulsion, 15,300 nuclear disintegrations 
or “ stars ” were observed. Measurements of scattering and grain density 
were made on 2000 of the tracks associated with these stars, the only 
criterion for selection being that the track, in the emulsion, was longer 
than 3 mm. 

Measurements of grain density and of the angular distribution of the 
ejected particles were also made on the tracks of 3070 shower particles, 
and 1508 grey tracks including those which were shorter than 3mm. In 
addition, all stars in which the track of the primary particle was lorig 
enough to allow the energy to be determined, were analysed in detail. 


§ 4. COMPOSITION OF STARS. 


The distribution of the 15,300 nuclear disintegrations according to the 
number of shower particles, n,, and the number of heavily ionizing particles, 
N,, is shown in Table I. The results given are those recorded by the 
observers ; but a more detailed study of individual stars, selected at 
random, indicated that the errors in the observations of both N, and 
n, were small. 

It will be convenient to describe first the distributions in the values of 
the energy, mass and angular distribution of the different types of 
secondary particles; and where possible, the variation of these quantities 
with the star parameters N,, and »,. 

4.1. “ Shower Particles’. 

Of the 2000 tracks for which the scattering parameter « was 
’ determined, 470 were shower particles. These were sub-divided into three 
groups—according to whether they could be identified as mesons, or 
protons ; or whether their nature was uncertain. We denote particles of 
the three types by the symbols z,, p, and z, p, respectively. 

z,. Mesons with kinetic energy in the interval from 80 to 1100 MeV. 
The particles in this group constitute about 60 per cent of the 
total number of shower particles. 

p, Protons of energy between 500 and 800 MeV. (20 per cent of 

' the total). x 

a, p. Particles with a mean angle of scattering, «, less than 
0-0226 degrees per 100 x. For such particles, the grain density 
in the tracks, both of mesons and protons, is equal to the 
minimum value, and the two types of particle cannot be dis- 
tinguished in the conditions of the experiment. (20 per cent of 
the total.) 

4P2 


TABLE I. 


ANALYSIS OF 15,300 STARS OBSERVED IN PLATES FLOWN AT 68,000 FT. 
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The distribution in energy of the 7, particles for stars of different types 
indicates that there is no significant dependence on either N, or n,. In 
other words, the energy distribution of the mesons is nearly the same for 
disintegrations in which a heavy nucleus of silver and bromine is completely 
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disrupted [as for those in which only a few heavily ionizing particles 
emerge, many of which are due to the disruption of light elements such as 
carbon or oxygen; and it is the same for “ showers” in which many - 
mesons are created as for those in which there are only a few. This feature 
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is displayed in Table IT., which shows values of the average kinetic energy 
of the particles for stars with different values of N,, and n,. It may be 
emphasized that although we are delaing here with the low-energy end of 
the energy spectrum of the ejected mesons, it contains more than 60 per 
cent of all these particles produced in the disintegratons. 


TABLE II. 


Average kinetic energy of the 7-mesons of energy E, where 
1100 MeV >E>80 MeV for stars of different types. 


Ns {| 2-4 D—900 
E 385 +40 340 +35 395 +40 
N; 3-8 9-14 15-00 
E 385 +45 360 +40 360 +30 


A similar analysis of the distributions in energy of the protons, p,. 
and unidentified particles, (7, p), is subject to large statistical errors, but 
the results are given in Tables ITI. a and III. 6. These show the propor- 
tions of protons and unidentified particles, (7, p), among the shower 
particles of all types, for different kinds of stars. The most important 
fact established by the results is that the proportion of these energetic 
particles (7, p) is greater for large than for small showers. If we consider 
all the stars in which a given number of shower particles, n,, is produced 
the proportion of fast shower particles (7, p) is greater for the small than 
for large values of N,. 


TaB_eE III. a. 
Number of z-p particles per shower particle. 
Nz, 3-8 Nz; 9-2 N; 3-00 
Ns 1 0-15+0-:04 0-:09+0-05 0-13+0-03 
Ns 2-4 0:22+0:05 0-22+0-:05 0-22+0-03 
Ns 5-00 0-41+0-:09 0-14+0-04 . 0-23-+0-04 
Ns 1-0 0-23+0:03 0-16+0-:02 0-20-40-02 
TaBue III. 6. 
Number of proton shower particles per shower particle. 
N; 3-8 Nj, 9-2 Ni 3-00 
Ns 1 0-21+40-04 0-18-+0-06  0-20+0-04 
Ns 2-4 0-16+0:04 0-13+0-:04 0-14-40-04 
Ns 5-0 0:05+0-03 0:12+0-03  0-10+0-03 
Ng l—o 0-16+0:03 0:13-0:02 0-14+0-02 


We shall see later that the first of these facts is a simple consequence of 
the greater average energy of the primary particles of stars with greater 
multiplicity ; and that the second is due to secondary interactions and 
loss of energy of the shower particles in escaping from the nucleus in which 


they are produced. 
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In order to calculate the average energy required to produce the 
x, p particles, it is necessary to make assumptions about the proportions 
of z-mesons and protons among them. This has been done in the following 
way : the energy spectrum of the -mesons in the cosmic radiation, most 
of which are produced by the decay in flight of z-particles, has been 
determined by Sands (1950). From this spectrum, the energy distribution 
of the parent z-particles can be computed ; and thence, from the low 
energy distribution of 7-mesons observed in the present experiments, the 
total number to be expected among the (7, p) particles can be calculated. 
The result is about two-thirds of all the (7, p) particles, and the rest are 
therefore attributed to protons. Assuming that this proportion is constant 
for stars of all types, the average energy of a (7, p) particle is found to 
be 1930 MeV. 

Using this value, and taking into account the rest-energy of the 
z-particles, the average energy of the shower particles can be found for 
stars of different types. The results are shown in Table IV. ; they reflect 
the different proportions of fast shower particles (7, p) in different kinds 
of stars. 


TABLE LV. 


Average total energy of shower particles (7,+p,+7, p) as a 
function of N,, and n,. 


Ng il 2-4 5-00 

E 750+40 MeV 850+40 MeV 910 +45 MeV 
N; 3-8 9-14 15-0 

E 855-+40 MeV 775 +55 MeV 730+45 MeV 


4.2. Angular Distributions. 


The angular distribution of the shower particles with respect to the 
line of motion of the charged primary particles has been determined by 
observations on more than 1000 “‘stars”. The results have been 
separated into a number of classes according to the characteristics of the 
stars and are shown in figs. | and 2. Most of the shower particles are 
emitted in directions making angles, 0, less than 40° with that of the 
primary particle, but some are ejected backwards, @>90°. The angular 
distribution of the backward-moving particles is nearly isotropic. 

Fig. 1 shows that the “width” at half-maximum of the angular 
distributions decreases—but only very slowly—as n,, the number of 
shower particles in a star increases ; fig. 3 shows the variation with 
multiplicity, n,, of the median angle of the distribution, n, for stars of 
a given type. The observations can be easily interpreted in terms of the 
increase of the average energy of the primary particles required to produce 
stars of increasing multiplicity. This involves an increase in the velocity 
of the C-system of the interacting nucleons, and therefore an increasing 
tendency for the mesons to be thrown forward, in the L-system, in 
directions near that of the primary particle. | 
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The ‘ half-widths”’ of the angular distributions increase with the 
value of N,, (fig. 3), an effect which can be interpreted in terms of inter- 
actions of the shower particles with the parent nucleus. Such interactions 
will contribute, first, to excitation energy of the nucleus, and thus to the 


Pipe 


No. of tracks. 


60 80 100 120 140 160 


20 


Angle to primary. 

‘The angular distributions of shower particles in the laboratory system for 
different ranges of multiplicity. The ordinate gives the number of tracks 
per angular interval. 
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Angle to primary. 
The angular distribution of shower particles for stars of N, <6 and N;, >7. 


completeness of the disintegration of the nucleus, 7. ¢. to an increase in 
N,,, and, secondly, to a scattering of the emerging mesons, and so to the 
breadth of the angular distributions. 

The fraction of the shower particles projected at angles greater than 
90° to the direction of the primary, is denoted by B. Observations on 
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this quantity also have an important bearing on the question of secondary 
interactions of the shower particles in escaping from the parent nucleus, 
and it will be convenient to refer to it as the “ backward fraction ”’. 
As shown in fig. 4, it decreases only slowly with increase in the value of 
n,, the multiplicity of the star. There is, however, a strong dependence 
on N,. Thus, taking all stars with shower particles, n,>1, the value 
of B is equal to 0-113--0-016 for N,<7; and to 0-1710-012 forN, >7. 
Almost all stars for which N, >7 must be due to the disintegration of the 
heavy elements in the emulsion, silver and bromine; and also some, 
at least, of those for which N, <7. It follows that the ratio of the backward 
fraction due to the disintegration of the heavier elements, By, to that for 
the lighter elements, B,, must be greater than 17-1/11:3, 4. e- 
By/By > 1-5 £0-24. 


Fig. 3. 
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Multiplicity n,. 
The variation of the median angle of the showers with the multiplicity, for 


stars of N;,>7 and N,<7. The curve shows the variation expected on the 


seat ies of multiple production in a single nucleon-nucleon collision 
(K=1-0). 


4.3. “ Grey” Tracks. 

The grey tracks are due almost entirely to protons, deuterons and 
tritons, and the results of a study of the distribution in energy and the 
relative proportions of the different types of particles have been given in 
Part IV. These quantities appear to be constant for stars of different 
types, with the exception of type O,, and fig. 5 shows their average 
number per star for different values of n,. 

From a knowledge of the number and the average kinetic energy of the 
particles producing the “grey” tracks the average energy required to 
produce the heavily ionizing particles which appear as mh grey’ and 
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. black ” tracks has been calculated. Following the procedure described 
in Part I. (Brown et al. 1949) the total energy Ey, of these particles, 
including neutrons of the same speeds, may be calculated. The results 
may be expressed by the realtion :— 


Ey =/(N,)=155 N,—100 MeV. (N,>1). 


This formula applies for stars of all types with the important exception 
of type O,.. The formula given in Part I., viz : Ky ,=37N,+4Nj, was 
obtained using a selection of stars which were predominantly type O,, for 


low values of N,. 


B 
—_Nuy7 => 1540.25 
By <6 


Fig. 4. 
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The variation of backward fraction with multiplicity. Experimental points ©. 
Curves marked “‘iso”’ calculated variation on the assumption of pure multiple 
production with isotropic angular distribution of the mesons in the C-system 
of the colliding nucleons for (1) and (2) generations. “ cos‘ @’” variation 
calculated on the assumption of an angular distribution in the C-system of 
form cos! 6’sin 6’d6’. The full line curves indicate the fraction projected 

« backward on the assumption that each meson is elastically scattered by a 
nucleon in its exit through the nucleus, either by Compton-type scattering, 
or by scattering isotropic in the C-system of meson and nucleon, with 
cross-section independent of energy. 


The angular distributions of the particles which produce the grey 
tracks, divided into three groups according to the energies of the particles, 
are shown in fig. 6. The observed tendency for the fast protons to be 
emitted at only small angles to the direction of motion of the primary 
particle is in accord with the results of Crussard (1950), Harding (1951): 


and Osborne (1951). 
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4.4, “ Black” Tracks. 

The “black” tracks are emitted almost isotropically, and their 
observed energy “ spectrum ” is consistent with the view that they arise 
as a result of the evaporation of the nucleons from a highly excited 
nucleus (Harding, Lattimore and Perkins 1949). 


§5. Stars PRoDUCED BY PRIMARY PARTICLES OF Known ENERGY. 


In analysing the experimental material, a track has been identified 
as due to a singly-charged primary particle producing a star if it satisfies 
the following conditions :— 

(1) Its grain-density must be less than 16 grains per 50 y. 

(2) The track must be in the upper hemisphere. 


Fig. 5. 


<— Ny. 


<— N, protons 25-500 MeV. 


<— Protons 80-500 MeV. 
Protons >800 MeV. 
Mesons >1100 MeY. 


<— Protons 500-800 MeV. 


Average number per star. 


| 2 5 10 20 
Ns 


hae a sc} . Aare . 
Variation with multiplicity of the number of protons and mesons of different 
energies and of the number of heavily ionizing particles Nj. : 


(3) For stars of higher multiplicity, n,>2, the track of the particle 

; assumed to be the primary must be not widely inclined to the 
average direction of the secondary shower particles. 

(4) The energy of the assumed primary particle, if measurable by the 
scattering method, must be sufficiently great to have produced the 

_ Star and the shower particles. 
Of 2000 tracks on which scattering measurements and grain counts 
have been made, 200 satisfied these conditions. 
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It is well known that mesons produced in nuclear interactions may be 
emitted backwards, the direction of motion being nearly opposite to that 
of the parent particle. It is therefore possible to make errors of 
identification and to attribute a track to a primary particle when it was, 
in fact, due to an ejected meson—the true primary particle being a 
neutron. Such errors are most serious in the case of stars of low 
multiplicity. A star produced by a fast neutron in which a single fast 
meson is ejected vertically upwards, type 1,, may be wrongly classed 
as of type 0,. The proportion of such errors can be estimated in the 
following ways :— 


Number of tracks. 


Angle in degrees. 
The angular distribution of the “‘ grey ”’ protons as a function of energy. 


The fraction of secondary shower particles from stars of type 1,, which 
have directions of motion in the upper hemisphere, is found to be 0-15. If 
it is assumed that the angular distribution of the secondary fast. particles 
produced in stars of type 1, is similar to that in type 1,, then it follows 
that only 10 per cent of the “ primary ”’ particles of stars of type 0, are, 
in fact, secondary particles. The probability of wrongly identifying 
a charged primary particle in events from which several shower particles 
are emitted is small, because the average energy of “ backward ”’ moving 
shower particles is low so that the scattering measurements prove them 
to be insufficiently energetic to have produced the star. 
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A plot of the values of g and & for all the track of “ primary ” particles 
identified by these criteria is shown in fig. 7. The distribution in energy 
of the particles—deduced from the results shown in fig. 7, after correction 
for geometrical factors and with the assumption that all the fastest 
particles (7, p) are protons, is represented in fig. 8. 

It is possible to determine directly the distribution in energy of the 
7m-mesons among the primary particles in the interval <10? MeV. 
Extrapolating the results to greater energies, it was found that only about 
15 per cent of the fast primaries (7, p ; E<1-2 BeV.) are due to 7-mesons. 
A second estimate of this ratio was obtained by the following method :— 

The energy spectrum of the secondary (7, p) particles is known, and the 
cross-section for their interaction with nuclei is close to the geometrical 
value. It is therefore possible to estimate how many of the high energy 


Fig. 8. 
Integral spectrum of “‘ primary ” protons and z-particles. 


| pe eee ee 


100 


Intensity in arbitrary units. 


Energy in BeV. 
The integral distribution of the “ primary ”’ protons and z, p particles. 


(x, p) particles, produced in nuclear disintegrations in the plate assembly, 
will interact with another nucleus in the local matter. Such calculations 
show that only about 10 per cent of the fast primaries (7, p) have been 
generated in the plate assembly and, because of the short life-time of the 
7-mesons and the very weak interaction of u-mesons, most of the others 
must be protons. Further, some of the fast (7, p) primaries locally 
generated are protons, thus only about 6 per cent of the fast primary 
particles of energy >10° MeV. can be attributed to 7-mesons. 

Similar calculations allow estimates to be made of the proportion 
of primary particles of lower energy which it should be possible to identify 
as m-mesons. The expected number thus obtained, for the particular 
conditions of the experiment, is 40, a number in good agreement with that 


observed, viz., 43. 
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5.1. Variations of n,, N, and N, with the Energy of the Primary Particle. 


The energy of sufficient primary particles has been ‘measured to allow 
the characteristics of the stars to be studied as a function of the primary 
energy, and fig. 9 shows the variation of n,, N, and N, with E,. 


* 1000 
Kinetic energy MeV. 


iA . . AY — _ . . . 

Variation of N;, Ny, and 7, with the kinetic energy of the primary particle ; 
x stars produced by 7-mesons ; @ proton initiated stars; O stars initiated 
by “a, p”’ particles. 

Fig. 10. 
Egrey Versus primary energy. 


Energy MeV. 


Primary energy in BeV. 


The average energy of a “ 
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Fig. 9 also includes the results of observations on showers produced by 
x-particles of great energy. For these events the mean multiplicity per 
nucleon has been plotted against the mean energy per nucleon. The 
results are in good accord with those for primary protons, with which 
they have been incorporated. Fig. 10 shows the variation with primary 
energy, E,,. of the average energy of the particles producing the grey tracks. 

The distributions of the individual values of n,, for a number of intervals 
in the value of E,,, are shown in fig. 11. There is a wide spread in the 
values of n, for a given value of E,,, so that the multiplicity of a shower 
gives only a rough indication of the energy of the primary particle which 
produced it. We have been unable, except rarely, in very favourable 
cases, to measure the energy of primary particles above 20 BeV., and very 
-few above 10 BeV. In analysing showers of high multiplicity we have, 
therefore, sometimes taken the observed value of n, as an indication of 


E,, but the results are then subject to large errors. 


5.2. Secondary Particles from Stars with Primaries of Known Energy. 


Whenever possible, detailed measurements have been made on the 
scattering, grain-density and range of the secondary particles arising in 
the disintegrations produced by primary particles of known energy. 
It will be convenient to consider separately the results for different 
classes of primary particles ; for protons, 7-mesons and 7, p particles :— 
(a) Primary protons with energy from 0-2 to 0-9 BeV. 

Among thirty-four stars produced by protons in this range of energy, 
only one secondary meson was identified. Allowing for geometrical 
factors, the estimated frequency of occurrence of the charged secondary 
mesons was found to be 0-07 -+0-07 per star. 

(b) Primary 7-mesons with energy between 0-2 and 1-1 BeV. 

In contrast with the above result for primary protons, seven secondary 
mesons have been identified in a study of thirty-seven stars produced by 
a-mesons, a result which, allowing for geometrical factors, corresponds 
to 0-46 0-23 secondary mesons per disintegration. Details of the 
observations on the seven secondary particles and the associated 
disintegrations are given in Table V. It will be seen that the mesons 
commonly emerge from these nuclear encounters with small kinetic 
energy. A similar observation has been made for 90 MeV. mesons by 
Bernardini, Booth, Lederman and Tinlot (1950, 1951). 

Pl. XLII. shows an event which is remarkable in that two secondary 
mesons are generated by the incident meson. 

(c) Fast primary particles: ratio of the cross-sections for elastic and 
inelastic scattering. 

When a fast particle interacts with a nucleus, some of its energy may 
or may not be expended in the production of mesons. Following 
Heisenberg (1949), we may refer to the two types of collisions as 
“ inelastic’ and “elastic”’ respectively. By determining the average 
multiplicity 7, for different values of E,,, and by making allowance both 
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Fig. 11. 


a KE=90—330 MeV. 


a KE=330—740 MeV. 


Number of events. 


a KE=740—1100 MeV. 


(a) 


p KE=240—535 MeV. 


p KE=535—810 MeV. 


Number of events. 


(7, p) KE=810—1240 MeV. 


(z, p) KE=1240—1910 MeV. a 
5 
5 
® 
3 
(7, p) KE=1910—4300 MeV. ie 
2 
= 
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(z, p) KE=4300—9400 MeV. 
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Ns 
(0) 


Distribution of multiplicity, ns, in stars produced by (a) z-mesons, (b) protons 
and “‘ a, p”’ particles, as a function of energy. 


for the shower particles which are protons and for the neutral mesons 
which are not recorded, we can estimate the proportion of elastic and 
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inelastic collisions made by primary particles of different energy. In 
making the analysis, it has been assumed that the mesons have a 
probability of 2/3 of being charged (Carlson et al. 1950). 


TABLE V. 
: Kinetic energy of Kineti f Y 
St ous gy netic energy oO Y% energy 
en ype incident meson (MeV.) emergent meson (MeV.) ; loss 
20, 79+8 67+13 15 
7+1, 640-+100 40+8 94 
6+1, 755 +150 450-+90 43 
5-+0, 825-180 17+3 98 
6+0, 900-190 29-+6 97 
3+1, 1040 +210 120+20 88 
375+70 

3+2,, f y 

+2, 1050-4200 365 4.70 70 

Fig. 12. 


Proportion F of 0, stars versus primary energy. 
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Proportion of ‘0p stars” plotted as a function of y,, the total energy of the 
primary proton in terms of its rest energy. The curve represents the 
values predicted by Fermi’s theory. (See text.) 


Fig. 12 shows the proportion of primary particles which do not produce 
shower particles, for different values of the primary energy E,. For 
E,,=1-6 BeV. (total energy y,=2-7), half the primary particles do not 
produce shower particles, and the curve shows the rapid fall of the cross- 


section for “ elastic” scattering as the energy of the primary particle 


SER. 7, VOL. 42, NO. 334.—NOV. 1951 4Q 


1258 U. Camerini et al. on the Nuclear Transmutations 


increases. The theoretical minimum value of E,, for the production of 
a 7-meson of energy 80 MeV., is 360 MeV. The results shown in fig. 12 
indicate that for values of E,,<1600 MeV. most of the collisions are not 
completely inelastic. 

The full curve in fig. 12 represents the values deduced from Fermi’s 
theory (Fermi 1950, 1951) (see also Part VII.). These predicted values 
take into account production of mesons of all energies, including those of 
less than 80 MeV. which are too slow to appear among the “shower 
particles”. Such particles are few, however, so that their inclusion does 
not seriously invalidate a comparison between the theoretical curve and 
the experimental results. The agreement between them is seen to be 
satisfactory. 


5.3. The Production of Neutral Mesons ; Missing Energy. 


It is now established that neutral as well as charged 7-mesons are 
produced in energetic nuclear encounters. An estimate of their frequency 
of occurrence has been made by comparing the total energy represented 
by the secondary particles emitted from stars of a given class with the 
mean energy of the primary particles producing them. It has thus been 
shown that there is some missing energy which it is reasonable to 
attribute to the creation of neutral mesons. 

The analysis was restricted to stars, with n,<4, produced by primary 
particles of velocity Bc, where 8>0-7. The upper limit to n, was applied 
in order to avoid the inclusion of large numbers of primary particles of 
energy greater than 10 BeV. If such particles were included, large errors 
would have been introduced in this analysis. The limit to the value of B 
was introduced because primary particles of lower velocity cannot be 
identified with certainty. The results obtained are summarized in 
Tabel VI. 

Column (1) shows the observed frequency of the different kinds of tracks. 
Values of the average kinetic energies of the particles in group (a) have 
been taken from the results of Harding et al. (1949) ; the other values in 
column (2) have been determined from the present observations. For 
the groups P, D, T, the energies quoted in columns (2), (3), (4), (5), are 
for each proton, whether ejected individually or as part of a deuterium or 
tritium nucleus. It has then been assumed that the total number of 
neutrons of energy greater than 25 MeV. is equal to 1-25 times the number 
of these protons. Column (6) gives the total number of nucleons, or 
uesons, associated with one ejected charged patticle (i. e. one track) of the 
type under consideration. 

The results show that the mean energy represented by the secondary 
nucleons, E,, is 1458+60 MeV.; and by the ejected mesons, 
K,, x(630+40) MeV., where x is the ratio of the sum of neutral and 
charged mesons to the charged mesons alone. The total energy associated 
with the “ average ”’ star is then E,=E,, + xE,. The average energy of 
the primary particles of this group, directly observed, was found to be 
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E,= 2700 +150 MeV. It follows that y=E,, — Equcteon/E,= 2:0 +0-3, a value 
which corresponds to a ratio of neutral to charged mesons of 1:0-+0:3. 
This result does not depend significantly on the assumed proportions 
of mesons and protons among the (7, p) particles, but only on the correction 
factor for the “ grey ” and “ shower ”’ particles. 

Carlson, Hooper and King (1950) have estimated the ratio of neutral 
to charged mesons produced in nuclear disintegrations, by observing the 
materialization of the y-rays into which the neutral mesons decay. They 
restricted their observations to showers with n> 3, average value 4-5, and 
found a value N,./N,+ =0-45+0-10. Gregory and Tinlot (1951), 
using a cloud chamber, obtained a value of ~0:3. 


TaBLeE VI. 
Type of No. per Kinetic Binding Total Energy pease 
crack oe energy energy energy per star a es . 
per track per track pertrack (uncorrected) "CUT 
particles 
Pep. F, 3:6 10 8 18 65 2:75 
os. 18 10 3 13 23 1 
4 P,D,T, 2°48 131 8 139 344 2°25 
BL. 0-29 47 140 187 41 x 
é ies 0-18 625 8 633 114 2-25 
Oa 0-60 395 140 535 321 x 
1 ay 0-063 1582 8 1590 100 2-25 
e ? 7 0-127 1960 149 2100 266 x 


All energies in MeV. 


The discrepancy between the present result for stars with n,=1-1 and 
that obtained for showers with n,=4:5, may be attributed, at least in 
part, to the fact that by selecting small values of n,, we also favour the 
selection of showers in which the ratio of neutral to charged mesons is 
large, and vice versa. In nuclear encounters which lead to the production 
of a given number of mesons, it may be assumed that the division between 
charged and neutral particles in any individual case will be subject to 
statistical fluctuations. If we therefore impose a restriction n,>4, we 
shall include in the sample stars with four charged and no neutral mesons, 
four charged and one neutral, etc. Because of the rapid fall in the number 
of stars as n, increases, such observations will tend to give a lower value 
of the ratio N,,/N,. 

ACKNOWLEDGMENTS. 

We have pleasure in thanking Professor C. F. Powell, F.R.S., for his 
constant interest and encouragement ; Miss M. Stott, Miss M. Merritt, 
Mrs. M. L. Andrews, Mrs. M. J. Harris and Mr. M. J. Harris for assistance 


with the measurements ; Miss J. Burnell and Dr. H. K. Heitler for the 


diagrams ; Miss P. Dyer and Mr. A. R. Gattiker for the photographs ; 
and the team of microscope observers of this laboratory for their careful 


scanning of the plates. 
4Q2 


Total 
energy 
per 
star 
179 
ee 

TUS 
41x 
256 
321yx 
225 


266 x 


1260 On the Nuclear Transmutations Produced by Cosmic Ray Particles 


This work has been carried out as part of a programme of research 
supported by the Department of Scientific and Industrial Research, and 
two of us (J.H.D. and W.O.L.) are indebted to this body for 


maintenance grants. 


REFERENCES. 


BERNADINI, G., Bootu, E. T., LEDERMAN, L., and Trntort, J. H., 1950, Bombay 
Conference on Fundamental Particles, 1951. Phys. Rev. 82, 105. 

Brown, R. H., CaMERINI, U., Fowusr, P. H., Herrter, H., Kine, D. T., and 
PowkE LL, C. F., 1949, Phil. Mag., 40, 862, (I). 

CAMERINI, U., Coor, T., Daviss, J. H., Fowuer, P. H., Lock, W. O., MurRHEAD, 
H., Tost, N., 1949, Phil. Mag., 40, 1073 (II). 

CAMERINI, U., Fowutr, P. H., Lock, W. O., and MurruHap, H.. 1950, Phil. 
Mag., 41, 413 (IV). 

Carson, A. G., Hooper, J. E., and Kine, D. T., 1950, Phil. Mag., 41, 701 (V). 

CRUSSARD, 1950, Bombay Conference on Fundamental Particles. 

Fermi, K., 1950, Prog. Theor. Phys., 5, 570 ; 1951, Phys. Rev., 81, 863. 

Fow er, P. H., 1950, Phil. Mag., 44, 169. (III). 

Fow er, P. H., Menon, M. E. K., Powetu, C. F., and Rocuat, O., 1951, in 
course of publication. 

GREGORY, B. P., and Tintot, J. H., 1951, Phys. Rev., 81, 675. 

Haropine, J. B., 1951, Phil. Mag., 42, 63. 

Harpine, J. B., Lattmmore, 8., and Perrys, D. H., 1949, Proc. Roy. Soc. A, 
196, 325. 

HEISENBERG, W., 1949, Z. Phys., 126, 569. 

O’CEALLAIGH, C., 1951, in course of publication. 

OsBorne, L.8., 1951, Phys., Rev., 81, 239. 

Sanps, M., 1950, Phys. Rev., 77, 180. 


U. CAMERINI et al. 


Phil. Mag. Ser. 7, Vol. 42, Pl. XLII. 


ea a 
s _ ¥ 
‘g °e * 
a * 
?. -# 
@ 2 ts es 
° 
3 oan -§ 
i, yaar i” 
af, a, . 
. a 2 woe 
* a = 
. ° . 
SH goers 7 
. ee ae . 
" . ° ° 
es : 4 cat 
. e a oa! 
. “ . 
é F * 
Ae ens 32 & 
. 
¢ - 
» 
nf ee * o o 
Fe 1 
. 
.. : a eee 4 
° * eae 
A EN, 
“Oe ot & 
ares * 
ry oe 
*/~ . 
+ ¢: 
4 ; 
° 
v.- 
#5y 
, 
Ne Wes 
Si 
sot yea ce 
Coes As 
» A ‘ 
a ee So 
. oe 
. : Y s 
* Meee) 
> %@ 
. 
4 % Dikeis 3 
* © * 
e 
¢ 7 
rs : 
. 
ee 
. 


Vp 
A star of type 2+1p in which a meson of 120 MeV. kinetic energy is produced, the 
primary proton having an energy of 600 MeV. 
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Disintegration initiated by a 7-meson of kinetic energy 1085 MeV. in which two 


a-mesons of energy 365 and 375 MeV. are produced, together with a proton 
of 270 MeV. 
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CXXV. Nuclear Transmutations Produced by Cosmic Ray Particles of 
Great Energy.—Part VII. Interpretation of the Experimental Results. 


By U. Camerini, J. H. Daviss, C. Franzinertr, W. O. Loox, 
D. H. Perxrys and G. YExurtte.t, 
H. H. Wills Physical Laboratory, University of Bristol *. 


[Received July 30, 1951.] 


SUMMARY. 


The observations presented in the preceding paper are analysed and 
compared with the current theories of meson production. In making the 
analyses, account is taken of the interaction with the parent nucleus both 
of the recoil nucleons and the created mesons. 

It is found that the observations do not provide decisive evidence in 
favour of either plural or multiple meson production. The experimental 
results can be explained in terms of either the multiple theory with 
secondary generation of mesons, or the plural theory assuming many-body 
collisions. 

Some other conclusions arising from the work are given at the end of 
this paper. 


§1. INTRODUCTION TO THE ANALYSIS. 


THE simplest conditions for studying the generation of mesons in 
nucleon-nucleon collisions would be provided by experiments in which 
fast nucleons collided with protons in a hydrogen target, but such 
observations have hitherto been impracticable. Although such collisions 
must occur in a photographic emulsion, they are rare, and we are therefore 
forced to attempt to infer the nature of the complex physical processes 
which take place in the nucleus from the results of the collisions of the 
nucleons with the mixture of nuclei found in the emulsion. In addition 
to hydrogen, which is the most common element, this mixture contains a. 
variety of nuclei which can conveniently be divided into two groups ; 
the heavier, mostly silver and bromine ; and the lighter, carbon, nitrogen 
and oxygen. 

In approaching this problem, the basic assumption is made that the 
disintegrations can be described in terms of encounters between individual 
pairs of nucleons. Unfortunately, the experimental data available do 
not permit a definite proof of the correctness of this assumption. If the 
physical processes taking place cannot be described in terms of such 
collisions, the interpretation of the experimental results becomes very 
difficult, for no theoretical guidance, regarding many-body problems is 


yet available. 


2 Communicated by Professor C. F. Powell, F.R.S. 
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If, however, the disintegrations can always be described in terms of a 
number of elementary collisions, nucleon-nucleon and meson—nucleon, 
the problem is reduced to an enquiry into their characteristics, and we 
have to determine the following quantities : 


(i.) The cross-sections for radiative and elastic scattering, and their 
variation with the energy of the incident particle. 


(ii.) The probability of radiating as mesons a fraction, K, of the energy 
of the two colliding particles which is available in the C-system, and the 
dependence of K on the primary energy. 


(iii.) The distributions in ‘energy, and in angle of emission, of all the 
secondary particles involved; for radiative collisions, the number of 
mesons created ; and the dependence of all these distributions on K, and 
on the primary energy. 


The analysis proceeds as follows : 


The experimental results are first compared with predictions based on 
the pure plural theory of meson production (Heitler—-Janossy 1949, 1950) 
in which it is assumed that the mesons are produced individually in 
interactions between pairs of nucleons. Secondly, a similar comparison 
is made for the pure multiple theories (Heisenberg 1949, Fermi 1950, 1951) 
in which it is assumed that many mesons may be produced in a single 
collision. An attempt is then made to reduce the observed discrepancies 
by considering both for the plural and the multiple theories the secondary 
generation of mesons by recoiling nucleons and, for the multiple theory 
only, the effect of meson—nucleon interactions. 


§ 2. PurzE PLuRAL THEORY. 


The plural theory, as originally stated by Heitler and Janossy (1949) 
was based upon the following assumptions : (a) in each nucleon-nucleon 
collision only one meson is produced ; (b) the mean free path for collision 
of a nucleon in the nucleus is of the order of the inter-nucleon distance ; 
(c) the nucelus is completely transparent to the created mesons. 


(a) Comparison of observations with the Heitler—Janossy Theory. 


In the pure plural theory, the total number of mesons, charged and 
neutral is given by the number of encounters, N,, made by the primary 
particle in traversing the nucleus, for it is assumed that one meson is 
created at each collision ; it follows that 


N= 28 eS be U 


loach ot Live Saat a oem eet) 


In this formula y, is the total energy of the primary particle in terms 
of its rest-mass energy and y, is the “ cut-off ” energy in the same units, 
below which the primary cannot generate mesons. o, assumed to be a 
constant, is the fraction of the kinetic energy lost by a primary particle 
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in each collision ; Heitler takes c=0-25. The recoiling nucleon takes a 
fraction « of the primary energy, the meson (s—«). If we consider only 
those mesons for which y,>1-7 (shower particles), then the cut-off is 
given by 

M, 


Meo 


(c—a)R(y,—l=17; R= 
The variation of n, with E,,, taking n,—2N, is then determined by o and « 
only. Heitler assumes further that «=c—«=} and y, then equals 3-1. 
The value of N, given by relation (1), multiplied by a factor 2 is plotted 
in fig. 1, and shows good agreement with the experimental results for 
values of n,<6. 


Bion l 


MULTIPLICITY-PRIMARY ENERGY RELATION 
HEITLER THEORY ——— o singly charged primaries 


If | 1 
6 GENERATIONS e@ alpha primaries 


Variation of average multiplicity with primary energy ; comparison of the 
experimental results with the predictions of the plural theory. yp=total 
energy of the primary proton, in terms ofitsrest energy. For the «-particles 
the multiplicity per nucleon is plotted against energy per nucleon. 


Messel (1951) suggests a value of (s—«) ~ 1/3, and that the remaining 
energy is divided symmetrically between the two nucleons. The multi- 
plicity using these parameters is somewhat higher than the values shown 
morig. Ls 

Heitler and Janossy suggest that at high primary energy the recoil 

-nucleons can themselves produce mesons, thus accounting for showers 
of very high multiplicity. The contribution of the recoils is shown in 
curve (2), fig. 1. 

According to the above picture of plural production, the ratio of mesons 
to recoil nucleons, generated in meson producing collisions will be less than 
or equal to unity. We cannot of course measure directly the number 
S.aq of recoil nucleons originating exclusively from meson-producing 
collisions. However, if we consider all nucleons emerging with energies 
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above 80 MeV. we shall certainly obtain an upper limit for Sas 
Table I. and fig. 2 we have displayed as a function of n, the ratio of 
charged plus neutral mesons {=3/2 (no. of identified mesons-+ 0°66 X no. 
of (7, p) particles)}, to the number of nucleons above 80 MeV. energy 
{—2x(no. of protons between 80 and 800 MeV.+0:33xno. of (7, p) 
particles)}. The data are taken from fig. 5 of Part VI. 


TABLE I. 
No. of mesons Nj. of nucleons .. (nucleons 
Ng including : 80 MeV Ratio ( ) 
neutrals 2 : ee 
l 1-65 3-30 0-50 
9 9-4 4-60 0-52 
10 13-0 10-6 1-23 
20 26:8 14-2 1-9 
Fig. 2. 
30 
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ene at high multiplicity (n,>10), the ratio exceeds unity, and thus 
ere appears to be more than one meson produced per radiative collision 
This conclusion would not be affected by choosing a somewhat lower 
limit (é. e. <80 MeV.) for the energy of nucleons recoiling from radiative 
collisions. The results in Table TI. depends also on the proportions of 
mesons and protons in the (7, p) particles (see Part VI.). However, any 


variation in these proportions cannot i 
affect the ratio of m 
by more than 30 per cent. a i 
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(b) Modifications of the Plural Theory. 


In an attempt to meet these difficulties, Heitler and Janossy (1950) 
and Heitler and Therreaux (Heitler 1950, 1951), have modified the plural 
theory in the following manner. The primary nucleons and the successive 
recoil nucleons are assumed to be confined within a very narrow cone. 
It is then visualized that at a certain stage in the disintegration, n moving 
nucleons may all be interacting with the same nucleon “ at rest ” giving 
n further mesons, but only one more recoil nucleon. After m stages of 
this process, there will be m(m-+1)/2 mesons, with only (m-+1) recoil 
nucleons. With this assumption Therreaux finds that the percentage 
of recoil protons amongst the shower particles decreases with increasing n,, 
and that the number of “ grey” recoil protons is almost constant for 
3<n,< 20. 

The degradation of energy as such a cascade develops will, however, 
lead to an increased angular spread of the emitted particles. If so, the 
true picture lies between the two extremes; that in which the ratio of 
mesons to recoil nucleons is equal to m/2, and that of the original theory, 
where there are equal numbers of mesons and recoil nucleons. 

A ratio of mesons to recoil nucleons lying somewhere between these 
extreme values is certainly compatible with our observations. On the 
assumption of a contribution from such many-body collisions therefore, 
the predictions of the plural theory are in accord with the experimental 
data. 


§3. Purge MuutreLte THEORY. 


The multiple theory of meson production, in its simplest form, assumes : 
(a) that the mesons are created in a single nucleon-nucleon collision; and 
(6) that the nucleus is transparent to both the recoiling nucleons and the 
mesons. Heisenberg (1950) has calculated the variation of N, with y,. 
He assumes that the mesons are produced with an energy distribution 
in the C-system of the colliding nucleons which is independent of y,, 
except that it extends up to a maximum value which varies with the 
total available energy in the system. This assumed spectrum is of the 
form f(y,)dy,=dyi/y2, where y, is the total energy of a meson in the 
C-system in terms of its rest-mass. If the theory gives an adequate 
picture of meson production, the form of the spectrum can be determined 
from the present experimental observations, provided certain conditions 
are satisfied. The method is as follows : 

In the case of 23 showers, it was possible to determine the energy both 
of the primary particle and of one or more secondary shower particles 
by observations on the scattering parameter «. Assuming that all the 
primary particles are protons and that all the shower particles of a 
particular event are mesons produced in a single nucleon—nucleon collision, 
a Lorentz transformation enables the energy and direction of motion of 
the mesons in the C-system to be determined. It is also assumed, 
tentatively, that any loss of energy and change of direction of a meson 
in its escape from the nucleus can be neglected. 
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In making the transformation, it is convenient to measure energy in 
terms of the rest energy of the particle, and velocity as a ratio to that of 


light. 
Let y,=total energy of primary particle in the L-system ; 
y,=total energy of primary particle in the C-system. 


It follows that : 
Yp=VU%p+1)/21. 
The total kinetic energy of the two nucleons in the C-system in terms 
of the meson rest-mass is, before collision, 


2R(y,—1) 
and the amount radiated as mesons is 
e=2R(y,—1)K, 


where K=degree of inelasticity, and R, the ratio of the masses of the 
proton and meson, equals 6-7. 
The maximum energy which a single meson can take is either 
2 
«’=6, OF pe ae srlad 4 58 Rouen a 
Yp 
whichever is the smaller. 

Now y;, and 0’, the energy of the meson, and its angle of emission with 
respect to the line of motion of the primary, in the C-system, are related 
to the corresponding quantities y, and @ in the L-system by the transfor- 
mations :— 


Ve=Ve¥p— V(¥2—-1). V(y2—1).c0s 8, . . . « (3) 
sin 6 / 

tan 0’ = ———~—____,, .......... (4 

Yp(COs 6—B,,/B.) ( ) 

bark 0k co" ns ee a rd 


yp(cos 4’+-B,,/By) ” 


Hach secondary is given a geometrical weighting factor g (@), which is 
introduced to take account of the relative probability of observing long 
tracks for different values of #6. The analysis is confined to those stars of 
which the energy of the primary particle has been determined, and this 
can only be done for long tracks nearly parallel to the surface of the 
emulsion, Shower particles ejected at angles near 6=0° and @=180° 
with respect to the line of motion of the primary particle will therefore 
have a much greater chance of remaining in the emulsion than those with 
§=90°. The observations therefore lead to a series of weighted values 


of 6’ and Ys, the angle of emission and the energy, respectively, of the 
meson in the C-system. | 
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The angular distribution of all such mesons is given in Table IT. 


Taste II. 
0’ (in degrees) 0—30 30—60 60—90 
Number Observed 34-23 70 +30 99+35 
Isotropic distribution 27 74 102 


Within the large statistical errors, it is consistent with an isotropic 
distribution. The ratio of backward to forward-moving mesons is 


0-95 +-0-44, a result consistent with the value expected from symmetry 
considerations. 


Fig. 3. 
2 3e 
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n 
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e 
The energy spectrum of the mesons in the C-system derived on the assumption of 
a single nucleon-nucleon collision. ; is the total energy of the meson in 
terms of the rest energy. 
@ K=0-5; OK=1.0. The curve represents the spectrum p*dp/y*. 


Collecting together all values of yj, we obtain a distribution in 
energy f(y,)dy,;. This however cannot be compared directly with 
Heisenberg’s spectrum which corresponds to infinite available energy in 
the C-system of the interacting nucleons. To make the necessary 
corrections to the observations, we divide by the factor p(y,) representing 
the observed fraction of showers having primaries such that y, is greater 
than the value given either by y,=2RK(y,,—1) 


or by Riyp —1)/% 
whichever is smaller. 

The corrected distribution f(y;) dy, is shown in fig. 3, for K=1 and 
for K=0-5, The peak at y,~2 is presumably a consequence of the 
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conservation of momentum in the C-system. The distribution can be 
represented analytically by the equation 
prdpl ys, 

where p=\/(72—1), which falls off at high energies as 1/y,?._ This is in 
satisfactory agreement with the spectrum chosen by Heisenberg. 

It may be mentioned that, even if the pure multiple theory is 
inadequate, valuable information may be obtained by making such 
transformations in favourable cases. 


3.1. Yield Curve on Pure Multiple Theory. 
The average number of mesons produced in the collision is 


cages € 
N= nS 
sR Ve) 
where 


(dm f vstrrarn | | flveidye ©. +. 8) 


The number of mesons classed as shower particles (energy in L-system 
Y,>> 1:7) is at: $e. 
N,= xN,, 

where the factor y is nearly unity, and is easily calculated. In reaching 
this result, an isotropic angular distribution of the mesons in the C-system 
has been assumed. If one-third of all the mesons produced are neutral, 
the average number of mesons which appear as shower-particles is 
n,=(2/3)N,; and , is therefore determined as a function of e=2RK(y,—1). 
The results for three values of K, 1-0, 0-7 and 0-3—are shown in fig. 4. 
It will be seen that with low values of y,, the results indicate a small 
mean value of K, contrary to the Heisenberg assumption, K=1. For 
such values of y,, some of the collisions will be completely elastic—see 
fig. 12, Part VI.—and the remainder correspond to a finite value of K, 
which must, as shown by energy considerations, be near unity. If only a 
single meson-producing collision occurs in the nucleus, a mean value of 
KX, near to unity, must be assumed to obtain agreement with the 
experimental results for large values of y,. 

The yield curve according to Fermi’s theory is included in fig. 4. 
At high energies, the predicted average multiplicity appears to be too 


small, though presumably the agreement could be improved by a suitable 
choice of his fundamental parameter Q. 


Fluctuation of n, in stars produced by protons of given energy. 


The results represented in fig. 11 of Part VI. which show the large 
fluctuations in the values of n, for stars produced by primary protons 
of a given energy may be due to :— 


(4) variations in the kinetic energy of individual mesons in the 
C-system ; 
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(b) variations in the degree of inelasticity, K, of the collisions ; and 
(c) statistical fluctuations in the ratio of charged to neutral mesons. 


Assuming on the average, that 2/3 of the mesons are charged, the 
effect of (c) is readily calculated. The spread then predicted in Fermi’s 
theory is found to be in good agreement with experiment. 

The effect of (4) was computed numerically on the Heisenberg theory, 
for different values of K. Assuming a spread in K, and mean values of K 
as deduced from the yield curve (fig. 4), there is again found to be 
satisfactory agreement with the observations. 


Fig. 4. 
MULTIPLICITY— PRIMARY ENERGY RELATION 


HEISENBERG THEORY 
K= tO 


1) 
— K=07 t GENERATION 
=O 2) ONS 


20 35 


¥P 

Variation of average multiplicity with primary energy ; comparison of the 
experimental results with the predictions of the multiple theory. 
Yp= VL(vp +1)/2]=energy of nucleon in C-system. 


3.2. Angular Distribution of mesons in the Laboratory System. 

The. question may now be asked “ How far does the assumption of a 
single meson-producing collision, without secondary interactions of the 
mesons in their escape from the nucleus, explain the observed angular 
distribution ?”’. It has been found by calculation that the form of the 
observed angular distribution of the mesons in the L-system cannot be 
uniquely explained in terms of simple assumptions regarding their 
energy and angular distribution in the C-system. On the other hand, 
the median angle y, and the backward fraction of the emitted mesons B, 
as observed in the L-system, do not depend so critically on the assumed 
features of the distributions in the C-system. Thus assuming only that 
equal numbers of mesons are projected backwards and forwards in the 
C-system, the medain angle is given by equation (5) 


Nana d DA (ieee, eR ro hea 


or since f7,/B,~1, tan 7=1/y,,.. The value of tan 7 given by this equation 
is an upper limit because a small proportion of mesons, emitted backwards 
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with low velocity in the C-system, will appear in directions 6<7 in the 
L-system. 

The variation of 7 with »,, as calculated from equations (6) and (7), 
is shown in fig. 3 of Part VI. The observed median angle decreases 
much more slowly with increasing n, than the simple theory predicts, 
and this is evidence against a pure multiple process. The difference in 
the observed value of B for N,,<7 and N,,>7 also indicates that secondary 
interactions play a role. 

If the angular distribution of mesons in the C-system is of the 
form F(6')d6’=cos"6'd6’, the fraction of mesons with velocity £; which 
travel backwards in the L-system is given by the equation 


B=5|1- (BE) | Be>Bo- | heey 


If shower particles only are considered (y,>1-7), the above formula 
is slightly modified. : 

Using the spectrum f(y{)dy, we can therefore find the value of B as a 
function of y’, or y, by numerical integration. Again using the y/, versus 7, 
relation, we get a curve relating B with average multiplicity 7,. 

It is found that the actual slope and peak of the curve does not depend 
very strongly on the value of n ; the absolute magnitude however, varies. 
almost linearly with (n+1) as can be seen by binomial expansion of (8). 
The curves (fig. 4, Part VI.), were calculated for K=1 and n=0 and n=4 
respectively. The slope of these curves is quite different from the 
experimental B versus n, relation. The correct slope can be secured by 
assuming secondary nucleon-nucleon collisions (see 4.1) and in order to 
account for the absolute magnitude of B we have to postulate an 
anisotropic component in the C-system angular distribution. Alternatively, 
it may be assumed that elastic scattering of the mesons by nucleons 
can occur, a possibility which is considered in more detail in § 5. 


§ 4. MutrreLe THEORY WITH SECONDARY PROCESSES. 
4.1. Secondary nucleon-nucleon collisions. 


The agreement between the experimental observations and the 
predictions of the pure multiple theory may be considerably improved 
by assuming that a second generation of mesons is created. Such mesons 
may be produced in collisions involving nucleons recoiling from the 
primary particle. In considering such effects, we first neglect any 
secondary interactions of the mesons in escaping from the nucleus. 

If such secondary generation occurs, a simple Lorentz transformation 
is no longer valid, for the mesons are created in several collisions involving 
nucleons of different energy. If, however, we make a statistical study | 
of a large number of showers, produced by primary particles of similar 
energy, it is sufficient to assume an effective value of Vor Von Pde 
mean value of y;, and y/o, the values of y;, for the first and second 
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generation, suitably weighted according to the number of mesons arising 
-in each. If it is assumed that the nucleus is so large that all meson- 
producing collisions which are energetically possible actually occur, then 
it may be shown that y/*/y/ is nearly independent of the value of K: 
and that it lies within the limits 0-8 and 1. 

The above considerations indicate that the average values for y, for the 
secondary particles (y;),,, as deduced from the spectrum of § 3 in which 
it is assumed that they are produced in a single collision, are at most, 
10 per cent too large. The errors are thus not greater than those due to 
the poor statistical weight of the observations at present available. It is 
therefore reasonable to estimate the proportion of mesons due to 
secondary generation, using the value of (y;),,,, and the results are included 
in fig. 4. The calculations have been made for three values of K=0-3, 
0-7 and 1-0. 


4.2. Angular distributions including effects due to secondary 
generation. 


As previously mentioned, the assumption both of secondary generation, 
and of anisotropy in the angular distribution of the mesons in the 
C-system, is needed to explain the observations on the backward 
fraction, B. The calculations have been made assuming K=1, and 
neglecting the spread in multiplicity n, for a given primary energy y,,. 
It is reasonable to suppose however, that these simplifications will not 
seriously affect the general trend of the theoretical curves. For a given 
value of y,, the average meson energy in the C-system will be greater 
than (1) a, if n,<n,; see equation (6). The backward fraction will therefore 
also be greater. The reverse is true if n,>n,. Thus the general slope of 
the theoretical curve should remain almost unaffected. 

When K<1, calculation of the backward fraction is very complicated, 
for the nucleons from the first collision recoil at an angle to the direction 
of motion of the primary particle. This must lead to a broadening of 
the angular distribution of which it is difficult to estimate the magnitude. 


4.3 Variation of Multiplicity n,, with star size N,,. 


/ 


If the secondary generation of mesons plays an important role at high 
primary energies, it will be more likely to occur in heavy than in light 
nuclei. Stars with N,>8 must almost always be due to the disintegration 
of silver or bromine nuclei, whilst stars of smaller N,, may be due either to 
the disintegration of nuclei of the light elements carbon, nitrogen or 
oxygen, or to those of silver and bromine. If therefore we consider stars 
produced by primary particles of the same energy, secondary generation, 
more likely in heavy nuclei, would contribute to making the average 
multiplicity of stars with N,,>8 greater than those with N,<8. On the 
other hand, any reabsorption of mesons in the parent nucleus will tend to 
reduce the number of mesons and increase the number of recoiling nucleons. 
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In view of these considerations, we have plotted the average multiplicity 
n,, as a function of primary energy, Ep, for the two classes of stars ; see 
fig. 5. Included in the figure are the results of Salant et al. (1950) for 
primary particles of energy between 1 and 8 BeV.., and for those of energy 
greater than 8 BeV. (average energy 14 BeV.). At high energies, there is 
a significant difference between the average multiplicities of the two types 
of stars. It is reasonable to infer that secondary processes frequently 
take place in heavy nuclei bombarded by protons of great energy, and that 
re-absorption of the mesons is relatively unimportant. 


Fig. 5. 
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Observed variation of average multiplicity with kinetic energy of the primary 
particle in BeV, for stars of different prong number N,. 


It may be remarked that the results of Salant et al. relate to stars of 
N,=2—8 and N,S9, while our figures are for stars of N,—4—8, and 


N r>9. This may explain the difference between the two sets of observa- 
tions at 2:4 BeV. 


§5. Mrson—Nucteon INTERACTIONS. 


Hitherto any interaction of the mesons with the parent nucleus has been 
neglected, and an attempt has been made to explain the experimental 
results in terms of nucleon-nucleon collisions alone. The results given in 
5-2 (6) of Part VI., show, however, that meson have a strong interaction 
with nucleons, and we now consider the effects of such interactions. 
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5.1 Elastic Scattering of Mesons. 


In order to estimate the effect of the elastic scattering of mesons during 
their escape from the parent nucleus, we have considered two possibilities. 
In the first, the scattering process is regarded as analogous to the Compton 

effect, the relative probability of different angles of scattering being 
similar to that in electro-magnetic case. On being scattered, a meson may 
then be deflected through a large angle and may lose a large fraction of 
it energy. On the second assumption, the scattering is regarded as 
‘“ isotropic ”’ ; the cross-section for the interaction is assumed to be constant 
irrespective of the energy of the meson, and after the encounter, all direc- 
tions of motion of the meson, in the C-system of the interacting particles, 
are equally probable. 

In calculating the effect of such scattering, we have made the simplifying 
assumption, that the original angular distribution of the mesons can be 
represented by a 5-function ; 7.e. all mesons are initially projected in the 
line of motion of the primary nucleon. Assuming each meson to be once 
scattered in escaping from the nucleus, the fraction B’ of the mesons which 
emerge in the backward hemisphere can be calculated for a given value of 
the primary energy y,. The calculation depends upon the known dis- 
tribution in the values of y,, and the energy “spectrum ”’ f(y) dy;. These 
quantities and the relation between n, and y, allow the two curves in 
fig. 4 of Part VI. to be computed. The actual magnitude of B’ will be 
greater than that calculated since the initial angular distribution of the 
mesons is certainly quite broad. A slight error is introduced when the 
curves are taken as plots of B’ against n,, instead of n,(y,), because of 
the spread of n, for a given primary energy. This effect has already been 
considered in 4.2. 

It is reasonable to suppose that any backward fraction unaccounted for 

_by the Lorentz transformation is due to large-angle scattering, and this 
view is supported by the fact that the angular distribution of the mesons in 
the backward hemisphere is nearly isotropic. 

The backward fraction due to single scattering calculated by the above 
method is certainly an underestimate. Comparison of the observed and 
calculated values of B therefore allows us to obtain an upper limit (for 
n=1, equation (8)) for the probability of large-angle scattering, and the 
value thus obtained is 0-3. For the nuclei considered, this value indicates 
that the cross-section for the scattering process is of the order of 1/10 the 
geometrical cross-section per nucleon. 

In the scattering process, the meson projects a nucleon, and this particle, 
if a proton, produces a “grey” track. The above value of the cross-section 
is consistent with the observed slow increase of N, with n,. The ratio 
Q=B(N,>8)/B(N,<8)=1:540-25. If Q’ represents the ratio (B 
heavy nuclei)/(B light nuclei), then from the effect arising purely from the 
Lorentz transformation, one expects Q’~1-0; from scattering alone a 
value of Q’=1-8. Since the observed value does not represent complete 
separation of light and heavy nuclei, the observed value of Q’, 1-5, is in 
reasonable agreement with that predicted from scattering alone. 
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5.2 Reabsorption of Mesons in the Nucleus. 


The evidence presented in section 4.3 of Part VI. precludes the possibility 
that any substantial fraction of fast mesons created in the primary act, or 
acts, can be re-absorbed and annihilated by nucleons, in the parent nucleus. | 
The average total energy of such a meson is ~860 MeV., and the disappear- 
ance of this energy would result in a considerable increase in N ss Thus 
if strong absorption takes place N,, should increase rapidly with increasing 
n,, which is contrary to the observations. 


5.3. Generation of Mesons in Meson—Nucleon Collisions. 


In traversing nuclear matter, a very energetic meson may produce 
further mesons in two ways: (i.) First, it may be elastically scattered 
by a nucleon; and this nucleon may recoil with sufficient energy to 
produce one or more mesons in collision with another nucleon. 
(ii.) Secondly, in the scattering process, the meson may generate further 
mesons directly. 

A meson primary originates, of course, as a shower particle and, to be 
identified, it must have an energy less than 1100 MeV. Such a particle 
can have only a small probability of producing further mesons by 
process (i.) since it cannot transfer more than 880 MeV. in collision with a 
nucleon, and the average value is only ~ 500 MeV. We should therefore 
expect to observe only the elastic scattering of such particles, and the 
result of their collision with a nucleus should commonly be a “ star” 
without shower particles. 

As described in Part VI. (see Table V. and Pl. XX XVI.) one star has been 
observed in which an incident meson of 1050-200 MeV. energy produces 
two mesons of energy 375+70 MeV. and 365+70 MeV. and a recoil 
proton of 270-++50 MeV. On process (i.) it follows from the conservation 
laws that one of the mesons should be ejected backwards, which is con- 
trary to observation. Thus the event can only be explained in terms of 
process (1l.). 

By analogy with the nucleon-nucleon collisions which result in the 
creation of mesons, we may refer to the reproductive process (ii.) as 
“inelastic scattering”. Heitler estimates that for mesons in the range 
of energy under consideration, the ratio of inelastic to elastic scattering 
is ~10 per cent. This value suggests that among the 37 examples of 
stars produced by z-mesons which we have observed, there should be one 
in which the primary particle leads to the emission of two identifiable 
mesons ; and we have seen that one such event has been found. 

Although meson production occurs rarely in the nuclear interactions 
of 7-mesons of energy less than 1 BeV., it is reasonable to suppose that it 
becomes important at higher energies. This view is supported by the 
experiments of Rosser and Swift (1951) who have exposed photographic 
emulsions, at an altitude of 2860 m., above and below 30cm. of lead. 
They found that the ratio of charged primaries to neutral primaries for 
disintegrations with n,>2 was 0:86+0-21 above the lead, and 1-76-+0-48 
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below it. For the following reasons they suggest that the difference is 
due to mesons created in the lead :— 

The average energy of a primary nucleon which produces a penetrating 
shower of n,>2 is of the order of 4 BeV. If the observed showers are 
attributed to fast nucleons recoiling from the 7-mesons, then the incident 
7-mesons must have an energy greater than this value, i.e. >4 BeV. 
. The number of 7-mesons of such great energy is not sufficient, from this 
point of view, to account for the magnitude of the effect observed. On the 
other hand, the results are readily explained if it is assumed that 7-mesons 
of energy greater than 1 BeV. can produce showers of z-mesons directly, 
the cross-section for the process being approximately equal to the 
geometrical value. 


§ 6. PRopuctTion or Mrsons By ENERGETIC «-PARTICLES. 


It is reasonable to suppose that, in producing mesons by collisions with 
a nucleus, a relativistic «-particle will act as a beam of four nucleons. 
Included in fig. 1 are the results obtained in studying showers produced 
by «-particles of known energy. For such showers, the multiplicity 
per nucleon has been plotted against the energy per nucleon, and the 
agreement between these results and those for proton primaries gives 
support for the above view. 


§ 7. CONCLUSIONS. 


The conclusions which follow from the present experiments may be 
briefly summarized as follows :— 

(i.) At low energies (<1 BeV.) most of the nucleon-nucleon collisions 
are completely elastic. The “ inelastic’’ scattering cross-section does 
not become comparable with that for “ elastic ” scattering until an energy 
of about 2 BeV. is reached. 

(ii.) For a given value of the primary energy, a wide range in the 
number of charged mesons emitted is observed in different disintegrations. 
The multiplicity, n,, therefore gives an unreliable indication of the 
primary energy. 

(iii.) Assuming pure “‘ multiple ”’ production, the energy and angular 
distributions of mesons in the C-system of the two colliding nucleons 
can be calculated. The angular distribution thus found is approximately 
isotropic; the energy spectrum can be represented by p°dp/y* 
(p=momentum, y=energy). 

(iv.) The simple theory is of multiple production, with no secondary 
radiative collisions, is inadequate. To obtain agreement with experiment, 
it must be assumed that mesons are also produced in secondary nucleon— 
nucleon and meson—nucleon colilisions. 

(v.) The modified plural theory, which considers meson production 
in many-body collisions by both the primary nucleon and the recoil 
nucleons, also gives satisfactory agreement with the experimental data. 


4R2 
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(vi.) The observations therefore do not provide a decisive test between 
the two theories, and there appears to be nothing to choose between them 
regarding economy of hypothesis. The multiple picture introduces a 
strong-coupling theory, whereas on the plural picture it seems necessary 
to postulate many-body collisions. 3m 

(vii.) The number of neutral mesons ejected from stars of multiplicity 
between 0 and 4 is found to be 1-:0-L0-3 times the number of charged | 
mesons. ze 

(viii.) Re-absorption of fast mesons by the parent nucleus in which they 
are created is very improbable, although such a process may explain 
the very large stars (N,,) of small multiplicity, n,. 

The cross-section for elastic scattering of these mesons by the nucleons 
(Compton-type scattering) is ~10 per cent of the geometrical value. 

(ix.) On the contrary, single fast 7-mesons of energy less than 1 BeV. 
interact strongly with nuclear matter ; in ~50 per cent of the cases they 
emerge from an encounter with a nucleus having lost a large fraction of 
their energy. It is not possible, from the present observations, to decide 
whether in the remaining cases, the mesons are absorbed, or whether 
they suffer charge exchange, so that the outgoing meson is undetected. 
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ABSTRACT. 


Two cloud chambers have been used at a depth underground equivalent 
to 26 m. of water in a study of the penetrating cosmic ray particles. 
Pairs of associated penetrating particles have been observed originating 
in lead plates in the chamber and also entering it from above. It is 
considered that the process is one in which a p-meson produces a 
penetrating secondary of mean energy ~10° eV., the penetrating 
secondary being either a u-meson or a z-meson. The cross-section for 
the process is estimated to be about 5x 10-29 cm.? per nucleon for lead 
and of the same order of magnitude for sandstone. 


§1. INTRODUCTION. 


THE discovery of associated pairs of penetrating particles underground 
was first reported by Braddick and Hensby (1939) as a result of cloud 
chamber experiments at a depth of 60 metres water equivalent. It is 
probable that these events are not identical with penetrating showers 
as observed at sea level, for it is known that the primaries responsible 
for penetrating showers are so attenuated as to be negligible at this 
depth. . 

Pairs of penetrating particles have also been observed under thick 
absorbers at sea level by Shutt (1946) who interpreted them as pairs of 
p-mesons. 

The present work was undertaken to verify and elucidate these results. 
Cloud chambers have been operated in a cave at 26 m.w.e. depth, using 
simple counter control and a number of events involving associated 
penetrating particles (A.P.P.) were photographed stereoscopically. 
While the work was in progress underground phenomena involving 
A.P.P. have been reported by George. and Evans (1950), using 
photographic emulsions, and by Greisen e¢ al. (1951), using counter 
_ arrangements. The relation of these results to ours will be discussed. 

It seems likely that the events described are secondary to p-mesons 
and that the interaction of «-mesons with matter is greater than was 
thought. 


Na iether aa a et A 
* Communicated by Professor P. M.S. Blackett, F.R.S. 
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§2. EXPERIMENTAL ARRANGEMENT. 


The apparatus was installed in a cave at Stockport under 2600 gm./cm.? 
of rock, which was mainly red sandstone of mean density 2 gm./c.c. The 
vertical cosmic-ray intensity was about } of that at sea level. Two 
cloud chambers were used at different times, they were installed in the. 
same position with their centres 170 cm. below the roof. The whole 
cave was maintained at a temperature of 18° C. 

In the first series of experiments the apparatus consisted of a cloud 
chamber of normal design 30 cm. in diameter and 13 cm. deep,, with 
a lead plate 2 cm. thick placed across the centre. This plate was used 
to distinguish between electrons and penetrating particles, electrons 
usually being recognized by the production of electron showers on passing 
through the plate. 

In order to investigate further the properties of A.P.P. a multiple plate 
chamber was constructed. The chamber was 40 cm. x 39-5 cm. x 36 cm. 
in size and could contain up to thirteen 1 cm. lead plates. The chamber 
was illuminated through the back and the expansion was produced by 
the symmetrical lateral movement of two rubber diaphragms on opposite 
sides of the chamber. 

The selection system usually consisted of a simple 3-fold counter 
telescope. Only one penetrating particle passing through the central 
part of the chamber was therefore required to cause an expansion. 


§3. EXPERIMENTAL RESULTS. 


3.1. Cylindrical cloud chamber. Series I. 


Two thousand two hundred photographs were taken with the arrange- 
ment shown in fig. la. The counting rate of the 3-fold counter system 
was 8142 per hour, being } of that at sea level. The total time 
during which the chamber was available to record incoming particles 
(sensitive time) was 27-5 hours. 

Thirteen photographs of this series each showed two penetrating 
particles which were considered associated (the criteria for acceptance 
as associated will be stated later). 


3.2. Cylindrical cloud chamber with 3 cm. lead above. Series II. 


The second series of photographs was taken with the arrangement 
shown in fig. 1b. The 3 cm. lead absorber above the chamber was used 
to increase the probability of degrading any electronic component. 
The position of the counters was not exactly the same as before and the 
triple coincidence rate was 50-+-2 per hour. 

Two thousand six hundred photographs were taken, the total sensitive 
time being 52 hours. Twelve apparent A.P.P. were observed, four of 
which had originated in the lead absorber above the chamber. 
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3.3. Cylindrical cloud chamber with 3 cm. lead above and a split 
counter tray below. Series III. 

To increase the efficiency of detecting A.P.P. relative to random 

associations of particles, a split counter tray was placed below the 

chamber as shown in fig. 1c. This arrangement was sensitive to electron 
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Experimental arrangements. 


showers in addition to A.P.P. The counting rate was 8-0-9 per hour 
at sea level and 5+0-4 per hour underground. The total number of 
photographs taken was 2400 and the total sensitive time was 480 hours. 
Thirty photographs showed apparent A.P.P., of which twelve originated 
in the 3 cm. lead absorber above the chamber. 
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3.4. Multiple plate cloud chamber containing thirteen 1 em. lead 
plates. Series IV. 

For all the observations with the multiple plate chamber a simple 
threefold counter telescope was used. The arrangement is as shown in 
fig. 1d and the counting rate was 12-5 per hour. Ina sensitive time of 
eighteen hours, 215 photographs were taken, each showing at least one 
particle penetrating all the plates. 

One pair of apparent A.P.P. was observed to originate in one of the 
lead plates. One small penetrating shower was observed and one pair 
of A.P.P. originating in the sandstone above the chamber was found. 


3.5. Multiple-plate cloud chamber containing six 1 em. lead plates. 
Series V. 

In order to increase the visibility of the tracks and to make it possible 
to distinguish possible ‘‘ heavy’? tracks from contamination «-particles, 
alternate plates were removed. 

One thousand two hundred and seventy-six photographs, showing at 
least one particle penetrating all the plates, were taken with this 
arrangement. Seven examples of pairs of A.P.P. originating above the 
chamber and traversing the lead plates were observed. Two cases of 
pairs of A.P.P. originating in the lead plates in the chamber and one 
penetrating shower produced in a lead plate in the chamber were found. 
Two photographs showed three penetrating particles entering and trav- 
ersing the chamber. For ease of comparison an analysis of all the 
penetrating events observed is given in the table below. 


TaBe I. 
No. of apparent | No. of penetrating 


Total pee pairs produced showers produced 
time in hrs. |§ |———__—____ 


Inlead Externally; In lead Externally 


1s 


1 


§4. Discussion or Resuurs on Parrs or ASsocraTED 
PENETRATING PARTICLES. 
4.1. Random associations. 

It is necessary to consider the possibility that these pairs of A.P.P. are 
due to random associations. When the pairs are produced in the lead 
inside the chamber the probability of random association is negligible. 
When the pairs are produced above the chamber we impose the 
following criteria :— | 

(a) The pairs must be visually contemporary in the chamber and 

parallel or diverging downwards. 
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(6) On reprojection the pairs must appear coplanar to within about 2°. 

The number of random associations expected on criterion (a) is 
approximately 1 in 1000 for the cylindrical chamber and approximately 
20 in 1000 for. the multiple plate chamber. Criterion (b) reduces these 
figures by a factor of about 20. Many of the pairs rejected due to 
criterion (b) are probably genuine pairs of associated penetrating 
particles, for which the lack of coplanarity is due to the effect of 
scattering in the lead or sandstone. 


4.2. Knock-on electrons. 


Mesons traversing matter are known to give energy to electrons by 
collision (knock-on process), and it is necessary to consider whether any 
of the events appearing as cases of pair production are really events of 
this type. In the arrangements IV. and V., where the particles are 
observed to penetrate several plates without the production of electron 
showers, the probability of their being electrons is negligible. 

In the arrangements I., If. and III. the particles are observed 
penetrating only one plate in the chamber and the probability of their 
being electrons must be calculated. 

Consider the arrangement used in series IJ. and III. A meson 
accompanied by a knock-on electron produced in the 3 cm. lead above 
the chamber where both particles penetrate the 2 cm. lead plate inside 
the chamber, would probably satisfy the criteria for acceptance as a pair. 
The theoretical energy distribution of electrons knocked-on by the meson 
spectrum was obtained, and for each energy the probability of the 
electron emerging from the upper plate and passing without multiplication 
through the lower plate was calculated. For the final electron distribution, 
the scattering distribution of fig. 2 was calculated and compared with the 
observed scattering. 

It is seen that very few of the pair particles which show small 
deflections are likely to be knock-on electrons, although the few particles 
which are scattered through large angles could be explained in this way. 

A further indication that the pair events in general are not knock-on 
electrons is given by the fact that there is no correlation between the 
angle between the particles of a pair and the angle of scattering of 
the more scattered particle. In a knock-on process, a small pair 
angle would be correlated with high electron energy and small scattering 
angle. We may conclude that most of the pairs of associated penetrating 
particles, observed in series IT. and IIT. and originating in the lead above 
the chamber, are genuine. 

In series I. the same argument may be applied to pairs which originate 
in the sandstone roof, and we conclude that the pairs are genuine. 

In series II. and III. any electrons which originate in the roof would 
have to pass through two lead plates to qualify as components of pairs, 
and we consider that all the pairs enumerated are genuine. The possibility 
that electrons secondary to. meson bremsstrahlung might simulate 
A.P.P. has been considered and found to be negligible. 
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From the above considerations, it is clear that associated penetrating 
particles do exist. It is now necessary to consider whether the production 
of these A.P.P. can be explained by known processes. 


4,3. Nature of primaries. 

The nucleon component which is mainly responsible for the production 
of penetrating showers at sea level is highly absorbed with an absorption 
length ~100 gm./em.2. Thus the probability that a nucleon arrives at 
the underground station is ~e~?*. This is far too small to account for 
the observed number of pairs. 

The possibility that these events are 7-meson induced is excluded if 
it is assumed that the absorption length for 7-mesons is the same as that 
for nucleons. However, the behaviour of high energy 7-mesons 
(energy ~10® eV.) is not at all clear. It is possible to eliminate 
a-mesons from sea level as the primaries responsible for pairs observed 


Fig. 2. 
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No. of particles per 2° interval. 


° 2 4 6 8 Te) 12 \4 16 18 20 
Scattering angle in degrees. 
Differential scattering distribution for the more scattered member of a pair, 
originating in the lead absorber above the chamber, in Series IT. and IIT. 


underground on the following considerations. Let us assume as the 
most unfavourable case that all mesons with E>10" eV. at sea level 
are m-mesons since 7-mesons of this energy produced near the top of the 
atmosphere will not have had time to decay before they reach sea level 

Such 7-mesons would be 1/4000 of the total cosmic ray flux at ms 
level. In addition there are at sea level a small number of 7-mesons 
produced locally ; these have however limited range and decay to 
form p-mesons. In order to account for the observed rate of production 
of pairs of A.P.P. in terms of production by z-mesons, all the high energy 
m-mesons must reach the underground station and then must all interact 
in 3.cm. Pb. This is obviously impossible. 
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The above considerations eliminate nucleons and m-mesons as the 
source of pairs. Unless we assume that these events are induced by 
unknown particles the only particles which remain in sufficient numbers 
are p.-mesons. 

4.4. Nature of the pair particles. 


In these experiments the main indication of the nature of the pair 


particles is given by their scattering on traversing the lead plates in the 
cloud chambers. 


Fig. 3. 
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Comparison of the integral scattering distributions. 


In series I., II. and III. accurate scattering measurements were 
possible to within 0-1° and the scattering of each member of the pair 
in the 2 cm. lead block was measured using a low power microscope 
fitted with a goniometer eyepiece. The pair particles were divided 
into two groups, containing respectively the less and more scattered 
members (L.S. and M.S.) of each pair. The integral scattering distribu- 
tions, together with that for single particles, are shown in fig. 3. It is 
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seen that the less scattered members can be correlated with the single 
particles, i.e. the normal p-mesons beam present at the station. On 
the other hand, the more scattered members are obviously not the 
normal p-mesons. In order to determine whether the difference in the 
distributions is due to the method of analysis the scattering angles for 
single penetrating particles were measured, these were taken in pairs 
and the distributions for the more and less scattered members of these 
artificial pairs plotted. It can be seen that the difference in the 
distributions for the genuine pairs cannot be explained in this way. 
It is concluded that the more scattered member or penetrating secondary, 
as we will call it, is probably a »-meson of mean energy ~10° eV., although 
the possibility of its being a 7-meson cannot be ruled out. 

Observations in series IV. and V. show that as a rule the penetrating 
secondary has energy greater than 10° eV., only one case being observed 
where it has energy less than this. This result is not consistent with that 
just given. 


Tas_e I]. 


Cross-section for production of pairs of A.P.P. in lead. 


Metres of Pb.| Observed | Corrected | Interaction | Cross-section 
Series traversed |No. of pairs} No. of pairs| length in in 
by s accepted | produced | metres of Pb.| em.?/nucleon 


‘me 78 3 3 26 5:7 x 10-29 
TIL. 1150 9 30 38 4-0 < 10-29 


hs 89 3 3 30 5-0 x 10-29 


4-5. Cross-section for the production of pairs in lead. 


Assuming that the primary responsible for the production of a pair 
is a »-meson, and that the interaction is between the u-meson and a 
nucleon, it is possible to calculate a cross-section for the ‘process. 

In the results of series III. corrections have been applied for loss of 
efficiency due to (a) the split tray arrangement, and (b) the narrow angle 
pairs lost due to the finite size of the counters. 


Angular distribution of pair particles. 


A histogram of the angular distribution of pair particles originating 
in the lead absorber in series IT. and III. is shown in fig. 4. A correction 


has been applied for narrow angle pairs lost in the split tray experiments 
as mentioned previously. 


4.6. Pairs entering from outside the chamber. 


Quantitative results for these events are difficult to evaluate owing 
to the effects of scattering in the sandstone. A pair of A.P.P. originating 
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in the sandstone and being scattered through an appreciable angle before 
entering the cloud chamber will be rejected under the criteria for selection 
mentioned in § 4,1. 

The angular distribution is very similar to that for pairs originating 
in the lead, indicating that the events are of a similar nature. 
Cross-section for production of pairs of A.P.P. in sandstone. 

This determination demands a knowledge of the angular distribution 
and energy of the pair particles since these quantities determine how thick 
a layer of sandstone is effective. From what has been said previously 
concerning the nature of the pair particles it seems probable that the 
penetrating secondary is a meson (u or 7) of mean energy ~10° eV. 


Fig. 4. 
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Corrected angular distribution of pair particles, originating in the lead absorber 
above the chamber, in Series II. and ITI. 


Assuming a unique energy of formation of 10° eV. for the penetrating 
secondary an estimate of the cross-section for sandstone can be made ; 
this neglects the effects of scattering previously mentioned and so the 
value obtained for the cross-section is a lower limit. 

An estimate of the absolute lower limit to the value of the cross-section 
for sandstone production can be made by assuming that the penetrating 
secondary has infinite energy (instead of 10° eV.), in this case the whole 
thickness of sandstone above the equipment contributes to the flux of 
pairs crossing the cloud chambers. This calculation has been made 
assuming the Gaussian distribution of pair angles approximating to 


the histogram in fig. 4. 
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The cross-sections calculated are 1/3 of those given in Table IT. Thus 
the cross-section for sandstone production of pairs is in any case greater 
than 5x 10-39 cm.2 per nucleon. 

Apparent pairs produced in ar. 

Three cases in each of series I. and V. have been found in which pairs of 
A.P.P. appear to have originated in the air above the cloud chambers. 
The pair angles and probable angles of scattering in the lead plates 
inside the cloud chambers are such that the cases cannot be explained 
in terms of scattering in the roofs of the chambers. The probability 
that these cases are examples of the random association of two particles 
is low, but this possibility cannot be entirely ruled out. Further 
investigations must be carried out before this phenomenon can be 
established. 

Tasue III. 


Cross-section for production of pairs of A.P.P.’s in sandstone assuming 
a unique energy of formation of 10° eV. 


No. of single 
Series mesons travers- 
ing chamber 


No. of pairs Cross-section 
observed in em.?/nucleon 


2200 Be 10% 


2600 3-5 x 10-29 
24000 In x10 


1500 5 3 10-29 


§5. PENETRATING SHOWERS. 


In series IV. and V. two cases have been observed of penetrating 
showers of low multiplicity produced in the chamber. These showers 
differ from those observed at sea level in that the penetrating particles 
are emitted at large angles to the direction of the incident particle. 

Two cases have been found where three A.P.P. enter the cloud 
chamber, and in each case one of the penetrating particles suffers two 
large angle deflections in the lead plate. It is probable that these deflected 
particles are 7-mesons from a penetrating shower produced in the 
sandstone. Several cases of three or more parallel A.P.P. incident 
obliquely in the cloud chamber and penetrating a few lead plates have 


been observed, these are probably remnants of the penetrating particles 
present in extensive air showers. 


§6. COMPARISON WITH OTHER OBSERVATIONS. 


Braddick and Hensby (1939) using an arrangement similar to that for 
series I-III. at 60 mwe. found five cases of pairs of A.P.P. out of 
2000 photographs taken with a counter selection system requiring one 
penetrating particle. This rate is similar to that found in this 
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investigation. George and Trent (1951), using a cloud chamber containing 
two 3 cm. lead plates, have found two cases of pairs of A.P.P., out of 
1500 photographs, originating in the upper lead plate and traversing 
the lower one without deflection. Both these experiments were carried 
out at a depth of 60 mwe. These experiments show about the same 
cross-section as ours, and this result is in provisional agreement with 
the hypothesis that the pairs are produced by p-mesons. However, 
the statistical weight of these experiments is as yet very low. Experiments 
by Shutt (1946) at sea level using a cloud chamber under 88 cm. lead 
showed the existence of pairs of A.P.P. The cross-section calculated 
from Shutt’s results is ~10-®° em.2/nucleon for pair production in lead 
and agrees with that found in the present investigations within the 
statistical errors involved. George and Jason (1947), using a penetrating 
shower counter set requiring at least three A.P.P. at a depth of 60 mwe., 
found no coincidences in 45 days continuous recording. George and 
Trent (1949), however, using a counter set of lower selectivity, found events 
which could be explained in terms of the production of penetrating 
groups of mesons of low multiplicity, by mesons. A cross-section of 
5x 10-9 cm.?/nucleon has been calculated for the production of these 
showers in lead, this value compares well with that found in this 
investigation for the production of pairs. George and Evans (1950) have 
found stars and penetrating showers produced by u-mesons in photographic 
emulsions exposed underground. The cross-sections found for these 
events are ~10-?° cm.?/nucleon for stars and ~4-10-%° cm.?/nucleon 
for penetrating showers. The shower particles observed by George and 
Evans are 7-mesons and the most probable multiplicity is about three : 
a multiplicity of shower particles as high as eight has been observed. 
The two cases in the present investigation in which three penetrating 
particles are observed are probably events of this type. The majority 
of the pairs of A.P.P. observed in our experiments cannot be explained 
in terms of penetrating showers containing two penetrating particles, 
since the total cross-section for shower production calculated by George 
and Evans is only 4-10-®° cm.? per nucleon and showers containing only 
two particles form a small proportion of the showers observed. It is 
likely that the production of a pair of A.P.P. is an event of a more 
elementary type than the production of a penetrating shower. If heavy 
particles are associated with these events it seems certain that they 
would have been observed in photographic plate experiments, and we 
are led to infer that the interaction observed is one in which a »-meson 
produces only a single penetrating secondary. 


CoNCLUSIONS. 
(1) p-mesons produce penetrating secondaries on passing through 
matter, one penetrating secondary being produced in each event. 
(2) The cross-section for this process is ~5 x 10-29 em.?/nucleon for 
lead and of the same order of magnitude (certainly greater than 
5x 10-39 em.?/nucleon) for sandstone. 
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(3) The penetrating secondary is probably a «-meson of mean energy 
~109 eV., but there is a possibility that it is a 7-meson. 
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SUMMARY. 


It is shown that a method of direct quantization of field equations 
{completely avoiding the canonical formalism) is applicable to a wide 
class of fields with a non-local interaction. In §1 the extended source 
theory of Peierls and McManus is quantized and discussed. The 
S-matrix is computed by the method of C. N. Yang. A general rule for 
an immediate transcription of the traditional S-matrix elements into the 
new form is given. In §2 another type of a non-local interaction is 
‘discussed and it is shown that all the usual covergence difficulties may 
be easily avoided if we give up the strict localizability of charge, energy, 
-etc., and if we abandon the postulate of gauge invariance of the second 
kind. This is acceptable since it is possible to guarantee the integral 
conservation laws without the continuity equations and to secure a 
‘vanishing photon rest mass by other devices, not necessarily by the 
formal gauge invariance. In the limit A+0 where A is a fundamental 
length, the discussed field theory goes over into the traditional quantum 
electrodynamics. 


§1. QUANTIZATION OF THE PEIERLS AND McManus FIELD. 
Ler us discuss first the field described by Peierls and McManus (1948) 
with the following Lagrangian 
a Oe Lie ree Som aS a btage rah ani’ Sew L) 


where L'® is the usual Lagrangian for the electromagnetic and spinor 
fields without interaction while the interaction part L’ is 


tilt fe eA Te) ie tater oR ys ey 218 (2) 

where 
NOR SOD DG eee eral 2 ny oath in! eS) 

while 
Geen bie — wala), ay cen See = (4) 


where F(x) is a function of the argument 2,x,. In order to get a 
correspondence with the tranditional theory F(x) should be a delta-like 


function 
F(x) = 


Gn) — | di exp (wet,,) i UcA)s ie (0) 
ena peed PE Se a pd ab FO Ta a A A A Al i hel 
* Communicated by Professor R. E. Peierls, F.R.S. 
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where the Fourier transform f(k2, A) is a function of k,k, and of a 
parameter A with the dimension of length. We may assume, €. g. 


f (2, AN=1/1+ (A2Kz)? | or f(E, Aj=exp {—(A2k)* > =) = (8) 
which become em in the limit A>0. We notice the following properties : 
f d4aF (x)= 1): f(0,A)=1_ > ee 


whence we get the identity 
AOR A), Sk et aa ee 


where the index zero means an unperturbed field obeying the Laplace 
equation 


[JAO20"! Pel Se ee oe 
The interaction part of the Lagrangian may also be written as 
Li}  dhc{ aR x), cA (2) ee 


The Lagrange equations are 


PAs — ino oe he ee 
ne ave mebs= (1e[2) {yhs, APY eee 

a lt yi pt ad err 
Be, —nup=— (te/2){AB, py}. . 2. | (11) 
The above equations are invariant under the usual gauge transformation 

0A ‘ 
A,>A, — oe p> exp (—ieA), mR Big Ee 
where 

(A=0. ce ca 8 eee 


To prove the invariance we use the fact that A® is identical with A. 
Equations (11) may be replaced by the equivalent integral equations 
expressed by means of retarded and advanced potentials (compare 
C. N. Yang (1950) and G. Kallén (1950). 


A,,(2)=AR(a)+ J d4a’D**(a— a Fl (a ca a ee eek eee 


wh x )=yir(a) + (e/2t) [ d4n’S™**(w—a') fy"yp(a’), AR(a’)}, .  (13") 
(x)= P""(w) + (€/2t) J dé’ {AR(a’), p(a')y"JS"t(a’'—2), . (13"”) 
where D™=D—1D; S"=S—18, = ee 
A,(%)=A0(a)+ J d4x’D™¥(a— —x') 92 (2’), oy ee Ronn Lees 


yb (ae) yr ae (¢/2t) J d4a’S*(a—a'){y'p(a'), AR(@’)},  (15’”) 

(ae) =y(a) + (6/24) f dx’ {AR a’), ple’ jy Sea), . 15") 
where 7 D“’—D-+-4D, SY =8448" ees 
We may write equations (13’) and (15’) also in the form 

AR(a t)= Ak in Sehr: a)+ f dq Dk R ret (ady (a— x") 5 (&'). 
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By adding together the equations (13) and (15) we get 
Al(x)=AO* (a) + fd’ DPB (x—2' p(x’), . . . 6. . (IT) 
h(x) =f (x)-+ (e/22) f d4x'S(x—ar') fy"b(x'), AP (easier CL Ta) 
lw) =P (a) + (e/28) f dxf’), P(e’ yy"S(a’—a), . (17"”") 


where A(?), (9), % mean half the sum of the ingoing and outgoing 
waves, €. 9. 
ee AS gs ie ie mee 1S) 
By subtracting the equations (13) and (15) we get 
eee eee de DEP (ean (x), 9 ot se (19%) 
yp (a2) — p(x) = (e/2t) J dAx’S(x—2')fy"f(a’), A,(e')}, . . (19) 
p(w) —yp"'(a)—= —(e/24) J d4x' {A,(x’), f(a)" Sle’ —a). 
(19) 
The traditional Hamiltonian formalism is not adequate for the 
quantization of fields with a non-local interaction. On the other hand 
it is possible to quantize directly the fundamental integral equations (13) 
and (15) or the equivalent (17) and (19) by assuming the usual 


(interaction-free) commutation relations for the ingoing and outgoing 
waves. The same relations hold for A”, 4 and 4! : 


[AO(c), 00), Gate ta) gt Fb + tt 04a oe (20) 
Dn Tie (ee) el ee) ho") 


{yp (a), ys) (ar’) } = {op (av), yf) (ac") } 
HAM), Px =[AMe), Fe")I=0. (207) 


The commutation relations between the perturbed waves follow 
automatically from the field equations (13) or (15) or (17) which may be 
solved successively in a power expansion of the coupling parameter e. 
For the quantization of the field equations it is essential to symmetrize 
the expression for the product of % and A! since these quantities do not. 
commute with each other ‘in case of a non-local interaction. We should 
also use the symmetrized current 


Ose 2a le at eee eo (21) 


Although we do not know the exact form of the commutation relations. 
between the perturbed wave functions (as the exact solution of the field 
equations is not known), we may surely find some general properties of 
the commutators or anticommutators, e. g., the anticommutator between 
two perturbed spinor functions must be a covariant function of the 
argument x,,—2%,, : 

{ys,(x), pg/X") 1 Ung -0 ere ee te) 


yee 
\ 
482 
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which follows from the invariance of the formalism under space and time 
translations. Hence, for the perturbed wave functions taken at the 
same space-time point x, we have the same commutations relations as 
for the unperturbed ones 
U5(0)==S,5(0); 2 i ee 
since we may translate the common point x, to the infinite past (or future) 
where we have to do with the ingoing (or outgoing) waves only. 
This last property enables the proof of the validity of the continuity 
equation for the perturbed current j,. By multiplying (11’’) and (11’”) 


to the right or left by is and ys we get 
ay . — 
By, Je (UA aal ta), pa(a)} AP (e120; ae at fs (24) 


since the anticommutator between two field functions taken at the 
same point is (from (23)) a number and, thus, commutes with everything. 
Also the averaged current obeys the continuity equation 


ca 

5 ig HO ee ol, ng eee ee 

Ox, : (24) 
Due to (24) the problem of auxiliary condition is the same as in the 


usual quantum electrodynamics 


re) Fe 
ie, fu |=0. diy sa Ge 


We may develop the Lorentz condition in a Taylor series and show, 
as usual, that all terms of the expansion vanish (due to (11’) and (24’)) if 
two initial conditions 
oe nd) yee a a ey 26 
0X, ma ’ )| [t=tol 0x, Bie es MN tect - : ( ) 
are satisfied. 

By applying the method described by C. N. Yang (1950) we may 
construct the S-matrix for the extended source theory. Since the 
ingoing and outgoing waves satisfy the same commutation relations, there 
should exist a unitary matrix transforming the ingoing into the outgoing 
waves 


AM—S-1AMS, yor_s-ying = (a7) 
From (17), (19) and (27) we get 
[AD?S]=—}{S, fdte'DPB(x—a')j(e')}, 9. 9... (28") 
[p', S]=—(e/4t){S, fdtz'S(a—ar') fyb (x'), AR(x')}}, . . (28””) 
LW, SI=(€/4i){S, fate {AN(e'), Ha’ \Sla'—x)}. 28") 


From (17) and (28) the 8-matrix may be computed as a power series in 
the coupling constant e 


S=14g ERO ce ee, eee (29) 


n=1 
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The explicit computation of R® and R®) shows that the equations (17) 
and (28) are compatible in this approximation and yield the result 


Be 1) AeA pe tere ei) er aA 30") 
RAY=THRMYPALP+LQ, «2. . & .. 7(380"} 
where 
a ; an 
P= 5 Jfdte’dta” G(x"), j(e")JD®*a'—2"), . (30) 


Q= (ie?/8) ff d4x’dte” ([(x")y"S(a"—a'), yp (e')] 


HIM %e!), Sle’ —x")y"H'(2"y) {AOH"), AFH") 
(30) 


The S-matrix is obviously unitary in this approximation and differs. 
from the usual one merely by the fact that every A, has been replaced 
by the averaged A and every D has been replaced by the twice 
averaged D®"®. This is exactly the same replacement which converts 
the traditional equations of a local theory into our new fundamental 
equations (17) and (28). We may easily guess the form of the S-matrix 
elements in all higher approximations. They will differ from the 
traditional ones by the same replacements which convert the usual field 
equations into the new equations (17) and (28). However, in the 
S-matrix elements appear only the unperturbed potentials A‘ whence 
(taking account of (8)) the general rule for the transcription of the 
traditional S-matrix elements into the new form is very simple: we have 
to replace every D function by the twice averaged D®® while every A\ 
and every D® function (which will be obtained by taking a vacuum 
expectation value) remain unchanged. 

The very fact that it has been possible to construct a unitary S-matrix 
confirms the compatibility of the set of quantized equations (13) and (15). 

We have to examine finally whether the extended, source theory is 
able to improve the situation with respect to the convergence problem. 
The expressions of the vacuum polarization current and the photon 
self-energy remain unchanged (since the D function does not appear in 
the well-known tensor K,,,). On the other hand the expression for the 
self-mass of the electron derivable from the S-matrix is 


Sm=— 5 fatey"(S(a)D®* (x) +S(2)D(c))y" exp (—ip,,)- (31) 


If we assume a suitable function F(z) the first term will converge while 
the second term is left unchanged and diverges.* 


* A convergent expression for the electron self-energy could be obtained iff 
we used the integrations of the causal function De of Stueckelberg and Feynman 


in the complex plane instead of Schwinger’s D and D@) functions integrated as 
principal values. 
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§ 2, ANOTHER TyPE or A Non-Locat INTERACTION. | 


Let us investigate another type of non-local interaction derivable: from 
_ the following Lagrangian 


LL, a) 
where 
L’=s,A, oS Me le eee 
with i 
6, = (el2)Lb®, ye? 1,” 2 Ge a ey 


where b® and pk are the averaged spinor wave functions obtained with 
the aid of a general rule 


a (x)= { d’'G(e—z)a(x)) 5 1 ee oe 


G is again a delta-like function 


1 , 2 
G(x) = Qa [ d*k exp (th,z,)9(k. A), i-g 
where 
gig, 0) 1. 2S ee 
For a spinor field in vacuo we have 
yl OR) —9(—am3,' A) (O(n), 5 a 


where m is the rest mass of the electron so that the current s“)) in vacuo 
differs from the usual current merely by a charge renormalization 


e’=eg?(—m?, Xd). hie = si eee) 
Thus, it will be convenient to use g-functions with the property 

g(—m?, A)=1. OE ee ee ey 
For example we may assume | 


1 


ko )= ya E mR OF Mkin Y=exp {—AMKi+me)e}. (9) 


7) 
With such functions G the current s“? will be identical with the usual pen 
On the other hand the property § 1 (7) is lost : 


J d4xG(a)=g(0, A), : Riser aes a 


where g(0, A) is different from unity but tends to unity for A>0, 
The field equations derivable from (1) are 


LJA,=—s,, ey Fe (11) 
1 p(x) ee 
iar + mip(x)=te | dtx G(a—a' ybP(e')Av(w') . . (12) 


with a similar equation for ib. 
It should be noticed that, with this type of interaction, neither j,, nor 


: s 
satisfy the continuity oquations. 


7 


Aye This unusual feature of the field theory 
Seems acceptable since it is possible to guarantee the integral conservation 
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Jaws without continuity equations. In particular, the charge conserva. 
tion is:secured simply by the gauge invariance of the. first kind.* The. 
non-existence of the continuity equation for the current four-vector. 
should be interpreted as a fundamental uncertainty preventing a strict 
localizability of charge. 

The field equations (11) and (12) violate the postulate of gauge 
invariance of the second kind.* The resignation of this postulate seems 
also acceptable since we may replace it by a less stringent requirement 
that the experimental rest mass of the photon should be zero. This may 
be achieved either by formal renormalization or by the so-called realistic 
compensation (Jost and Rayski 1949). In this respect the new formalism 
by no means constitutes any drawback in comparison with the 
traditional quantum electrodynamics which was gauge invariant only 
formally so that we were also obliged to remove the (infinite !) self-mass 
of the photon by means of a renormalization or compensation. 

The equations (11) and (12) may be integrated by means of retarded 
potentials which yield 


A, (7)=Aj(a)+ (ie/2) | d4a’D**(w@—a')[P*(x'), y'P*(w')], .  (13) 
p(w) = (x) + (e/2t) Jdta'S® *\(a—ar') fy"pR(a’), A,(w')}, . (14) 
where the products have been written in a symmetrized form. Similar 


equations hold for the advanced potentials. By adding and subtracting 
the equations with the retarded and advanced potentials we get 


A, (2) =A (ax) + (te/2) fd4e’D(a—a')[P®(x"), -y"PR(a’)], (15) 
(x)= P'(a) + (e)24) [dt Sawa!) p(w’), A,(e')} (18) 


and 

Ai(az)— A(x) = (ie/2) [d4x’D(x—zx’) bP (a'), y'pB(a’), .. (17) 

(ar) — yp" (ae) = (e/24) [ d4a’S®(a—a’ ane Ale) } see (18) 

where 7‘ and A(” mean half the sum of the ingoing and outgoing waves. 
From (16) and ( 18) we get analogous formule for the averaged * 

PP (ar) = YOR (ar) + (6/21) fdta'S®*(w—a') fp (a'), A,(w')}, -  (16') 

po (ae) —p® M(x) = (2/21) Jax’ S™*(x—x') (ype (x'), A,(w')}. - (18') 

We quantize directly the field equations by assuming again that the 

ingoing and outgoing waves satisfy the well-known commutation relations 

for the separate fields in vacuo. The commutation relations for the 

perturbed wave functions will follow automatically from the field 


equations. 
' The S-matrix may be computed from the formule 


[A®, S]= — (ie/4)[S, [dt Diw—a')[PX(a’), ype’), «6 «+ (19) 
(HOR, S]——(e/4i)1S, fale’S"™w—a' YR’), Aw)» 20) 
- * For a definition of this, see Wentzel (1949), §§4, 11, 16. 
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with the aid of (15) and (16). There is little doubt in the compatibility” 
of our fundamental formule and we may immediately guess the S-matrix 
elements in any approximation. They will be quite similar to the usual 
ones with the only difference that every p', ° and S should be replaced 
by POR, POR, and S®®. Due to (6) and (8) the general rule for an 
immediate transcription of the usual S-matrix into the new form is simply :. 
every 8 should be replaced by S**. 

The new S-matrix elements yield non-vanishing probabilities only for 
processes in which energy, momentum and charge are conserved. Thus, 
integral conservation laws are secured in spite of the absence of the 
notion of densities satisfying continuity equations. 

As the perturbed current s, does not satisfy the continuity equation, 
the question arises as to the form of the auxiliary condition. In order 
to answer this question we introduce (15) into (16) and get a new equation. 


yb (a) = (w) +(€/2i) [dta’S®(a—2’) Gx 
x {ysp(x’), AO (a) + J d4a"D(x’—a")s,(a")} . . . (21) 


(and a similar equation for %) where only A®, but not A, appear 
explicitly. In the same way, by substituting (15) into (20) we get new 
formulae for the computation of the S-matrix in which only the 
unperturbed potentials A‘ appear explicitly. Now (15) may be treated 
as a definition and not as one of the equations. Thus, the perturbed 
potentials A, become secondary quantities without influence upon the 
fundamental equations. Now we may limit ourselves to the unperturbed 
A) and assume the Lorentz condition as a restriction upon the state vector 


0 
5 AD [= 0, 2g ee 


iu 
which reduces to two well-known initial conditions, as usual, due to the _ 
fact that the unperturbed potentials obey the Laplace equation 


AO a Lice oS ec 
The computation of the vacuum polarization current yields now 
«j,(X) vac —4e* i} d4r’K,,(w~—2)A,(x’), e . . be (24) 
where 
AAR GAY QA QAM /aA® AA y 
K = CO ES 5 SD en ee a R 
jw 0x, Ox, Ox, 02, wy ( Ox, Ox; mA 4°) < . (25) 
or 
(8,(%) )vao= —4e? fdte'L,,(a—a')A,(a'),  . . . (26) 


where L,, differs from K,, merely by the fact that every J is aver- 
aged twice, J**. The appearance of a 4® or 4®® changes essentially the 
situation in the problem of the induced current and secures convergence. 
The condition for the gauge invariance and the continuity equation for 
the induced current is, of course, violated as 
7] 0 
5— K,,=5*. ia, 7. st Bernal He Rae ONY LY 9 


* t 
0x, 
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is different from zero. However, the regularized delta function 
5®(x2)=G(a) ree a eee ee a (28) 
is ever and essentially different from zero only in the neighbourhood of 
zero while 0/dx,4 is odd, so that (27) vanishes if integrated over a small 
region about the origin. Hence, the continuity equation and the gauge 
invariance are still valid in the mean. The photon self-energy and the 
self-charge are finite so that the procedures of mass and charge 
renormalization are now mathematically correct and may be applied 
consistently. 
The electron self-energy derivable from the S-matrix is 
dm= —(e2/2) [d4xy"(S (x) D(x) +S®®(x)D (2) exp (—ip,a,), . (29) 
where the second term is convergent but the first term still diverges. 
By comparing (29) with §1 (31) we see that a full convergence may be 
achieved by combining both modifications discussed in this paper. The 
convergence will be generally secured if we assume the following 
Lagrangian 
= (te/4) {[p* (a), yp*(a)], AP (z)}, . 4 : “g (30) 
where % and ¢% are averaged by means of a delta-like function G(x) while 
Ae is averaged by means of another delta-like function F(x). The 
question of the definitive form of the delta-like function lies beyond the 
scope of the present investigation. We mention here only that other 
types of non-local field theories (e.g. the non-local field of Yukawa) 
may be quantized by the method of a direct quantization of field 
equations. An investigation of Yukawa’s non-local field will be 
published by the author elsewhere. 
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SUMMARY 


In continuation of a previous paper (Bishop and Hill 1951) it is 
~ econiectured that the work done in plastically deforming a polycrystal is 
approximately equal to that which would be done if the grains were free 
to deform equally. In conjunction with the principle of maximum plastic 
work, this enables the yield function of an aggregate to be calculated. This 
is done for an isotropic aggregate of face-centred cubic crystals, following a 
determination of the stresses needed to produce multi-slip. The theoretical 
yield criterion lies between those of Tresca and von Mises, in good agreement 
with observaton for copper and aluminum. It is shown further that the 
work-hardening of an aggregate would be a function only of the total 
plastic work if the grains hardened equally ; the departure from this 
functional relation is expressed explicitly in terms of the non-uniform 
hardening. 


§ 1. ResuME oF Previous Work. 


In a recent paper (Bishop and Hill 1951, henceforward referred to as BH) 
some general theorems were proved for an aggregate in which the crystals 
individually deform by slip according to the Schmid law. The theorems 
depend, in essence, on a principle of maximum plastic work for a homo- 
geneously deformed single crystal (BH, equation (9)). This principle 
states that if a crystal is caused to deform plastically through an incre- 
ment of strain de;,, the work done by the required stress o;; is not less than 
that done by any other stress o*, not violating the yield condition; thus 
(o4,—o7,) dep 2202 2 Se 
The yield condition is that the component shear-stress in any of the possible 
shp-directions, and over the associated slip-planes, cannot exceed the 
corresponding critical shear-stress in the current state of hardening. The 
critical shear-stress may vary from one slip-direction to another, due to 
differential hardening, or may be different in the two opposite senses of the 
same slip-direction, due to a microscopic Bauschinger effect. In a statis- 
tically homogeneous aggregate, in which cohesion between crystals is 
maintained by multi-slip, the same principle was proved for the 
macroscopic stress S,, and strain-increment dE,,. (BH, equation (16)) 
a 


t+ Communicated by the Authors. 
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From this it was shown to follow immediately that, if the yield criterion 
295 the aggregate is 

iC ip) =O yagtie deer eed a oa kcopnen aad 2) 
in a certain state of hardening and anisotropy, then the relation between 
the ratios of the components of the stress and strain-increment tensors is 


of 
dk, Bares te Sans aN ese et (3) 
‘The function f is, moreover, concave to the origin when regarded as 
surface in stress hyperspace, and does not involve the ace 
‘component of 8;;. A is a function of the stress and strain-history, con- 
‘trolling the magnitude of the strain-increment produced by a further 
Increment of stress. 


The plastic behaviour of an aggregate is thus completely determined 
by the functions f and h. The present paper is concerned with their 
calculation for an aggregate of face-centred cubic crystals. 


§2. THe Basic THEOREM. 


Since the actual stress distribution o,,; is in equilibrium, it follows by 
‘virtual work that 
S(%; dE ;;) dV=S,; dH;;, 


‘where the integral extends over a unit cube of aggregate (see BH, equation 
(14) and succeeding remarks). Let o% be a stress which, at each point 
-of the aggregate, would produce the strain dE;; in a free crystal having 
the local orientation and hardening. The distribution o* is, of course, 
not necessarily in equilibrium. Then, by (1), 

(o7;—o,;) dE, 50, 
‘the equality holding in general only when the deviatoric parts of o,;; and 
o* are equal. Hence é 
Si dH ,,=J(o,; dE A) dV <dk,,fo*, dv. . . . . (4) 


On the other hand, by integrating (1) through a unit cube of aggregate, 


we obtain 
S,, dE,,>J( (o% de;;) dV. 


Assuming that there can be no statistical correlation between any 
components of o%* and de;,, 


J(o% de,;) dV=Jo%, dV x J(de jj) dV= dH ,,\o%, dV, 
provided there is no rupture or relative sliding of the grains. Hence 
SAGA) SrA DPA Coeds Alpe SS Cus Sr ee Ara 29) 


On combining (4) and (5) we obtain the theorem on which the 
possibility of a calculation of the function f depends : 


adi, —dB,,jot drome tant. »(4) 
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In words : the actual work done is equal to the work that would be done: 
if the grains all underwent the same (macroscopic) strain. 

In particular, when the slip-directions in any one grain all have the 
same critical shear-stress + (in both senses), we may define an average 
value 7 for the aggregate, such that 


7=Jr dv. 


Then, if o% would produce a strain dE ;; in a grain with critical shear stress: 
zt, «t=7o0%/r would produce the same strain in a grain of the same 
orientation but with critical shear stress 7. Hence 

Tlo* dV=Jro%, dV=7la% dV, 
assuming that there is no correlation between o*, and 7. Therefore : 

foX% dV= J of dv, 
and so, from (6), 
Si dH ,,=dk,; i o* dV. . . . . . . . (7) 


We can thus evaluate the work using the average value of the critical 
shear stress in the aggregate. 

But equations (5) and (6) cannot be strictly true since (6) implies that 
_ the deviatoric part of o,; is equal to that of o*, all components of which 
are generally discontinuous across grain boundaries or wherever else the 
lattice orientation changes. (It is shown later that a general strain can 
only be produced by some one of a finite set of stresses.) However, we 
are inclined to think (though we have not found a rigorous proof) that 
equation (6) isnot much in error, and we assume it as an approximation. 


§3. MmtTHoD OF CALCULATION OF THE YIELD FUNCTION. 


The way in which f can be calculated is conveniently visualized geo- 
metrically. A strain increment dE;; can be considered as a free vector’ 
in stress hyperspace. According to (6) the scalar product §S,;dE,; is. 
calculable in terms of the slip properties of a single crystal. The extremity 
of the stress vector S;;, corresponding to a strain-increment dH,,, is thus. 
known to lie on a plane which is perpendicular to dE, and whose distance 
from the origin is 
dK ,,fo% dV 8 

(Although the point fo%dV lies on the plane, it does not, of course, 
necessarily coincide with §,;.) Now according to (3), the normal to the 
yield surface at the point S;; is parallel to dE iz, and se perpendicular to 
the plane. The plane is therefore tangential to the yield surface at 
S;,._ The surface is thereby found as the envelope of planes whose distances. 
from the origin are given by (8). Since possible strains have zero hydro- 
static part (being the result of simultaneous simple shears), the surface 
is cylindrical with generators parallel to the direction 8 Las 
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When the aggregate as a whole is isotropic the principal axes of the stress 
and strain-increment tensors are coincident, and we only need to consider 
the relation between their principal values represented in three dimensional 
Space. The aggregate will be statistically isotropic when the orientations 
of grains in the same state of hardening are randomly distributed ; this 
condition must be satisfied for each of the various states of hardening at 
any moment in the differentially-hardened aggregate (if not, the aggregate 
might be anisotropic despite the absence of a preferred orientation). 

The numerical calculations in the present paper are restricted to such 
an isotropic aggregate in which, moreover, the. slip-directions in any 
grain are equally hardened (so that equation (7) is applicable). The 
development of deformation textures and their effect on the yield surface 
is left to a later paper. 


§4. Muxti-Siie In A Face-CEnTRED CuBIC CRYSTAL. 

At ordinary temperatures and rates of strain, glide occurs in a face- 
centred cubic crystal on the octahedral planes in the directions of the 
octahedron edges. Each of the four distinct glide-planes contains three 
possible slip-directions, making twelve in all, each of which has two 
opposite senses. The positive senses of the slip-directions are arbitrarily 
chosen here according to the following Table; letters a, b, c, d refer to 
the four glide-planes, and with suffixes 1, 2, or 3, denote incremental 
shears in the respective positive senses. The components, referred to the 


Plane (111) (I11) (111) (111) 


Shear cme te OGD Pi Date 0g0 )- 625.69 Cg Oy, Og, vg 


Direction 011 101 .110 O11 101 110 011 101 TIO OIL. 101 110 


cubic axes, of the strain tensor de,,; due to simultaneous shears in the 
twelve directions are given by the equations (Taylor 1938 a) 


4/6 de ;=—4_+43—by+b3—Cy+C3—d2+ds, 
4/6 dégg= 4, —A, +b, —b3 +c, —¢3 +d, —ds, 
4/6 degg= —4,+4,—6,+b,—0 +Cg—4, + dp, 


ae) 
24/6 dégg= Ay—A3—b+b3+Cy—C3—d2 +s, 


24/6 des;= —4,+43+b,—bs +0,—C3—d, +s, 


. 24/6 deyo= A —g+b,—b,—C, + Cg—4, +g. | 
Any possible strain has five independent components (the hydrostatic 


part being zero), and therefore in general can only be produced by multi- 
slip over a group of directions containing an independent set of five. 
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Of the #2C,—792 sets of five shears, only 384 are independent. The 
remaining 408 include dependent combinations of type 


a@,+a,+a;=90 (144), 
a,+b,+d3=0 (228), 
a,—b,+¢,—d,=0 (36), 


or their equivalents ; these equations are to be interpreted as meaning 
that such combinations of wnit shears produce zero resultant strain. Of 
the 384 independent sets of five shears, 216 have two shears on each of 
two planes, and 168 have two shears on one plane only (the latter were 
apparently thought by Taylor to be dependent sets). ; 

If the components of a stress applied to the crystal are o;; referred to- 
the cubic axes, the shear stress components, multiplied by 1/6, are equal 
to 

A-G+H(q,), B+F—H(a),  C-F4+6(a), | 


A+G+H (b,), B—F-H(b,), C+F—G (), 


(10) 
A+G—H (c,), B+F+H(c,), C—F—G (cs), 


A—-G—-H (d,), © B—F+H(@,), C+F+G (4), 
where 
A=099—033, B=o33—04, C=01,;—%99, 
H=Ges d G==¢4; : Hea 


It will be observed that the 126 matrix of coefficients in the 
relations (10) between the 12 shear stress components and the 6 applied 
stress components is just the transpose of the 6 x 12 matrix of coefficients. 
in (9). This is a simple consequence of the virtual work equation 
g;; de,;=2X7 dy, and does not depend on any particular lattice geometry. 
It follows that we can always find a stress for which the component. 
shear stresses attain the critical values in prescribed senses in a given set- 
of five independent slip-directions (the respective minors being identical 
and non-zero). The critical value would usually be exceeded in one or 
more of the other seven directions, but, for any given strain, it is always 
possible to find at least one of the independent sets for which there 
exists a physically possible stress to operate the constituent shears. 
However, in evaluating the expression (8) there is no need to determine 
a physically possible combination of shears which are equivalent to the. 
strain, It is necessary merely to calculate the works done in the given 
strain by stresses not violating the yield condition, and to select from 
these works the greatest. In fact it is only necessary to make the choice 
from the works done by 56 particular stresses, which correspond ‘to the: 
‘* vertices’ of the polyhedral surface in stress space representing the 


yield criterion for the crystal. Proofs of these statements are left to a 
subsequent paper. 
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The 56 stresses may be classified in five groups, in which the typical 


values of (A, B, C, F, G, H)/4/67 are as follows when the critical shear 
stress 7 is the same in all the slip-directions : 


(i.) (1. —1, 0, 0, 0, 0). Tension or compression of amount 4/67 along 
a cubic axis. 


(ii.) (0, 0, 0, 1, 0, 0). A pure shear of amount 4/67 in a cubic plane 
and parallel to a Sil: axis. 

(ili.) (3, 3, —1, 0, 0, 44). A pure shear of amount 4/37 in a cubic 
plane and at 291° to the cubic axes. 

(iv.) ($, —%, 0, +3, +4, 0). Principal stresses + 4/67(1, —i, 0) with 
the second normal to an octahedral plane and the third along a slip- 
direction in that plane. 


(v.) (0, 0, 0, 3, +4, +4). Tension or compression of amount 34/67 
normal to an octahedral plane. 

The 56 stresses have 6 members in each of (i.) and (ii.), 12 in (iii.), 
24 in (iv.), and 8 in (v.). The critical shear-stress is attained in 8 slip- 
directions for (i.), (ii.), and (iii.), and in 6 slip-directions for (iv.) and (v.). 


§5. Tot Numericat MEeruHop. 


For a calculation of the yield function of an isotropic aggregate in which 
' there is no Bauschinger effect, it is only necessary to consider macroscopic 
strains dE,;, whose principal values are (1, —A, A—1) dE where $<A<1 
(Hill 1950, p. 17 et seq.). . The work done in such a strain of a free crystal 
has to be computed for each orientation of the strain axes to the cubic 
axes. If (0, 4, w) are the Eulerian angles of the strain axes with respect 
to the cubic axes, all essentially distinct orientations occur once only in 
the intervals 


0<$<n/4,  O<cot@<sing,  0<p<n, 


in view of the lattice symmetry and the assumption of equally hardened 
slip-directions. In other words, it is sufficient to restrict the major strain 
axis to one of the 48 identical spherical triangles in the stereographic 
projection, and allow the other axes to rotate through half a revolution. 
If dW is the work done on unit volume of crystal in a strain defined by 
given values of A, 0, ¢ and % (7 being the average critical shear stress), 
then, from (7), 


S,, dE ,(A)=J JJ dW sin 6 do dp dp/fff sin 6 dO dd dip. 


In the computation, net points were taken at 5° intervals of 6 and ¢, 
pug 18° intervals in . Five values of A were taken in steps of 3 between 

land 1. It was found convenient to express dW in terms of the principal 
PAlees of the 56 stresses and the angles between the stress and strain 
axes. The angles were read from a stereographic net on which were 
marked the axes of the 56 stresses. 
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If p denotes the perpendicular (8) from the origin to a tangent plane 
of the yield surface in principal stress space, the results of the calculations 


are 
1 


2 


5 te 
8 8 
V (8) pt 3:06 3-04 2-98 2-91 2-86 


The error involved in the integration is estimated to be not more than 
one unit in the second decimal place. 


§ 6. DISCUSSION OF RESULTS. 


Fig. 1 shows a typical 60° sector of the cross-section of the yield surface, 
obtained from the values of p as described in §3. The calculated curve 
cuts orthogonally the radii corresponding to pure shear and to pure 
tension, and lies between the Tresca hexagon and the von Mises circle 
when all are made to coincide for tension. In particular, the ratio of the 
yield stress in shear to that in tension is 2-86/(4/3 x 3-06)=0-540, 
compared with 0-500 and 0-577 for the Tresca and Mises criteria. 

The experimental data of Taylor and Quinney (1931) for copper and 
aluminium are also shown in the figure, and tend to lie between the 
Mises circle and the theoretical curve. The small discrepancy may be 
due to defects in the theory such as the approximation involved in 
equation (6), unequal hardening of slip-directions in individual crystals, or 
microscopic modes of distortion other than slip in the 12 directions. On_ 
- the other hand, the experimental data may not be completely reliable, 
due to anisotropy or to a slight uncertainty in the extrapolation used to 
circumvent the hysteresis loop. 

Also shown in the figure are the plastic potential curves computed 
from Taylor and Quinney’s measured relation between the directions of 
the stress and strain-increment vectors ; Taylor’s (1947) calculations for 
‘copper have been corrected in some instances. The agreement of the 
plastic potential curves with the directly measured yield loci suggests the 
substantial validity of (3), the theoretical derivation of which does not 
involve the further assumptions made in the present computation of the 
actual form of f. 

The comparison can also be made in terms of Lode’s diagram, but in 
our opinion, this over-magnifies small variations in the stress-strain 
relations, and is responsible for frequent statements that the Lévy—Mises 
equations are merely a moderate approximation. Thus, to the uncritical 
eye, a typical observed pair of Lode variables (0-5, 0-4), compared with 
the Lévy—Mises prediction (0-5, 0-5) might suggest an error of 20 per cent. 
However, it can be seen from the figure that the maximum difference 
in direction between the normals to the Mises circle and the measured 
plastic potential curves for copper and aluminium is only about 4° ; 


without any theory at all the strain-increment vector could be anywhere 
in the 360° range. 
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Fig. 1. 
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§ 7. INFLUENCE OF RESIDUAL STRESSES. 


When the external loads are removed from a plastically-deformed 
aggregate, a residual distribution of internal stress remains. According 
to the maximum work principle for the aggregate (BH, equation (16)), 
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the yield function for the aggregate does not depend in any way on these 
internal stresses, but only upon the intrinsic hardening of the grains (as 
specified by the current values of the critical shear stresses in the slip 
directions). This conclusion is directly due to the neglect of elasticity in 
the present theory. In effect, the theory assumes that the elastic moduli 
are indefinitely large, so that a single crystal remains rigid under 
increasing load until it yields plastically. Similarly, an aggregate remains 
rigid until all grains are stressed to their individual yield points, and there 
is therefore no hysteresis loop. During loading, the weakest grain becomes 
plastically stressed under a load which may be termed the elastic limit 
of the aggregate. The elastic limit is generally appreciably lower than the 
yield-point load of the aggregate and clearly depends directly on the 
residual stresses. With increasing load other grains in turn are plastically 
stressed, but are prevented from deforming (or hardening) by the remaining 
non-plastic, and therefore rigid, grains. Hence, apart from the restricted 
possibilities of a stress adjustment within the individual yield surfaces 
of such plastically stressed grains, the greater part of further increments of 
load is borne by the non-plastic grains. 

In an actual aggregate, elasticity permits plastic distortion and 
hardening of the weakest grains as soon as the elastic limit is reached. 
The immediate consequence is a hysteresis loop, greatly dependent on 
residual stresses (and correspondingly removable by a mild annealing), 
and with a breadth equivalent to a strain of elastic magnitude. When 
the hardening during the loop is small, it is naturally to be expected (in 
accord with observation) that the reloading branch of the loop will bend 
over fairly sharply to rejoin, virtually, the continuation of the previous 
stress-strain curve. Moreover, the local ordinate of the curve should be 
effectively identical with the yield-point stress calculated for a plastic- 
rigid aggregate. Many parallels might be instanced in stress analyses of 
structural parts on the basis of the macroscopic theory for a plastic-elastic 
homogeneous solid. 

Similar considerations apply to the hysteresis loop during reverse 
loading, and to the tendency of the residual stresses to act so as to lower 
the elastic limit (Bauschinger effect). However, a further distinction — 
must be made here. When the sign of dE,; is reversed, the value of S,, 
calculated from equation (6) would not be merely reversed when there is 
a microscopic Bauschinger effect ; that is, when the critical shear stress 
on a slip-plane depends on the sense of slip. Thus, even for the plastic- 
rigid model, one would expect in such circumstances a macroscopic 
Bauschinger effect on the yield-point. 


§8. Work HARDENING OF AN AGGREGATE. 

The calculated value of the tensile yield stress is 3-067. This agrees 
closely with the value obtained by Taylor (1938 b), who assumed that 
an aggregate actually deforms uniformly (whereas we have conjectured that 
the work done is as if it did). Taylor was able to determine the yield 
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stress directly from the work done in a pure extension since, for an 
isotropic aggregate, it follows from symmetry alone that the stress is 
uniaxial. In general, of course, it is also necessary to know the relations 
between the directions of the stress and strain-increment vectors. 

Taylor computed the work from the minimum shear principle. As we 
have proved (BH, equation (11)), this leads to the same result as the 
maximum work principle, but the latter is much more convenient. 
Taylor’s values for the works done in the individual orientations agree 
closely with ours (even though he selected the shears from only 216 of the 
384 independent sets). 

As Taylor showed, the formula S=mz(mE) for the tensile stress-strain 
curve of an aggregate, where 7(y) is the shear-hardening curve for a 
crystal, is in fairly good agreement with observation (when m=3-06, in 
the case of aluminium). It is assumed in the derivation that the virtual 
work equation S dE=7 dy can be written for an aggregate, where the 


averaged quantities + and 7 actually correspond on the shear-hardening 
curve. RK 


In general, however, 
dIW=S,,dE,=frZ|dy|dV and dt=f(dr) dV=fr(7)2 | dy | dV 


if hardening in multi-slip is a function only of the sum of the shears. 
Here r(7) denotes the rate of hardening dz/dy. Assuming that there is 
no correlation between 7 and 2'| dy|, 


dW=7{E|dy|dV and dz=r(7)fS| dy|dv. 
Hence 


dW=7 dzjr(r). . . ee ee (ID 


If all grains are equally hardened, r(7)=r(7) and dW=7dy.7 and 
the hardening of the aggregate (7. e. the scale of the yield surface) are 
then functions only of the total plastic work. The unequal hardening of 


the grains in such a way that r(r)Ar(t) is presumably responsible for the 
observed departures from this functional relation and from the formula 
S=mr(mE). 
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CXXIX. Note on the Boundary Layer on a Rotating Sphere *. 


By L. Howarrts, F.B.S., 
Department of Mathematics, University of Bristol +. 


[Received July 6, 1951.] 


SUMMARY. 


The solution of the boundary layer equations for the flow associated 
with a rotating sphere is considered and an approximate solution is 
obtained. The solution on either hemisphere is thought to be valid up 
to the immediate vicinity of the equator where it is invalidated by the 
interaction with the flow from the other hemisphere. 

A short appendix gives the flow due to a semi-infinite cylinder rotating 
in a uniform stream. 


§1. INTRODUCTION. 


Tue problem of the flow engendered by a sphere rotating about a 
diameter in otherwise undisturbed fluid has two features which are not 
without interest. In the first place it provides an example in which 
the curvature terms in the three-dimensional boundary layer equations 
are important. In the second place the boundary layers from the two 
hemispheres impinge on each other in the neighbourhood of the equator 
and their interaction is, in itself, of interest. 

The object of this note is to discuss in some detail the first of these 
aspects and to place the second on record as a problem still to be solved 
in detail. In explanation of this second problem, it may be said at once 
that the boundary layers on the two hemispheres may be regarded as 
originating at the poles (where the flow approximates to the rotating disk 
solution) and developing towards the equator where they impinge on 
each other. The boundary layer equations can deal adequately with 
most of this development but must fail to represent the region of 
interaction between the two impinging layers if for no other reason than 
the parabolic character of the equations. 


§ 2. THe Bounpary Layer Equations, 


We shail use spherical polar coordinates r, 0, 6 with r measured radially 
outwards from the centre of the sphere, 6 measured from the axis of 
rotation and ¢ the azimuth. In order to preserve w for the velocity 
normal to the surface (i.e. in the direction r increasing) we shall use 


* This paper forms a small part of one of the essays to which the Adams 
Prize was awarded in 1951. 


ft Communicated by the Author. 
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u, v for the velocities in the directions 0, 4 increasing respectively. The 
boundary layer equations can then be derived directly from the full 
equations of flow in spherical polar coordinates or from Howarth (1951) 
by placing K,=0, K,=—(1/a) cot 6, where a is the radius of the sphere. 
They are: 


u Ou Ou wv 02u 

a a6 me Oe ee aa Mi ae tr ee (1) 
U Ov Ov ow 02u 
Macnee a oh ey oree A re 2) 
TeGUa Ow Ka ) 
emer ue so eee orileetie (3) 


where use has been made of the fact that derivatives with respect to ¢ all 
’ vanish and that there is no imposed pressure gradient. 

Apart from the curvature terms (v?/a) cot 6 and (uv/a) cot 0, we should 
have the ordinary kind of two-dimensional boundary layer in w and w, 
with v determined by a convection equation. We may regard the (v?/a) cot 
term in (1) effectively as an imposed pressure gradient (whose intensity 
varies across the layer) essentially accelerating in character, since it is 
positive over the range 0 to 7/2 in 6 to which attention is restricted. At 
the surface, with 2 denoting the angular velocity of the sphere, 


v= Qa sin 0, 


so that 


2 . 
v : 
a cot 6= 22a sin 6 cos 6. 


The effect of this term is therefore a maximum when 0=7/4 and tends to 
zero as either the pole or the equator is approached. It is to be noted 
at once, in support of what has been said earlier, that it is most unlikely 
that forward integration of equations (1)-(3) from 6=0 to 7/2 will lead 
to a solution at 6=7/2 in which w=0 throughout the layer. In other 
words the boundary layer equations, which present us with the prospect 
of the layer from the two hemispheres meeting head-on at the equator, 
are inadequate to deal with the flow in this vicinity. Nevertheless, it 
is to be expected that the interaction will be limited to the immediate 
vicinity of the equator and that it will lead to an outflow in the equatorial 
plane of the same order as the mass flows through the boundary layers. 
Since the solution of the rotating disk problem is known to lead to inflow 
along the axis, we may reasonably expect inflow at the poles, and this 
combined with the outflow at the equator gives us a picture similar to 
the one envisaged by Stokes (1845). The arguments used here must 
however be regarded as essentially different from Stokes’s, his being 
based on pressure which, to the order of the approximation here, has 
been regarded as uniform. 
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§ 3. THE SoLuTION NEAR THE POLES. 


A solution of equations (1), (2), (3) in series in powers of @ is obtainable, 
though it is unlikely to converge sufficiently rapidly to make an approach — 
to the equatorial conditions practicable. To obtain it we put 


: 
BEN | 

é | 
u=—aQ[OF,+6F,+ ... |, = tee: ete Ee 


v=a Q[0G,+64G;+ ... ], 
w= (v2)"*[H,+0H;+ ...], 
and find the following equations for the first six functions 


F24+F/H,—G2=F*, | 


(5) . 


2F,G,+H,Gi=G;, 

2F+H/=0. 

4¥F,+F/H,+F5H,+4G6?7—2G,G,=F;, Mars 1. 
2F'.G,-+48,G,+-H,Gi+H,G,—4F,G,=Gs i ee ee 
AF 4 00 ee 


The boundary conditions are F,=F,=0, G,=1, Gz=—%4, H,=0, H,;=0 
at Z=0. 
F,>0, F,>0, G,>0, G;>0 as Z> oo. 

Equations (5) are, as mentioned earlier, the equations determining the flow 
due to a rotating disk Karman (1921). It is to be doubted whether, since 
the functions defined in equations (6), (7), (8) refer to this specific problem 
only, the labour involved makes their integration worth while. The 


significance of the expansion lies in giving the form of that solution near 
the poles. 


§ 4. THe SoLuTION BY THE KARMAN METHOD. 


We may look to the momentum integral equation to provide a more 
rapid though approximate solution. Integrating (1) through the boundary 
layer and making use of (3) we have, writing z=r—a, 


ue ‘aaa Ou 
She Se J 
| sputdeta | we | + | (u—v2) cot @ d= —va (5): 


that is 
0 


d /é 3 ou 
eal des | (we—vt) cot de=—0a (ZY , # Hi ee (9) 


since u=0 when z=0 and z=5, 
Similarly from (2) 


bee : + ov 
Re (uv) dz-a | mw | “fie jw cot 6dz=—va (=), 
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d ( dv 
a). uv dz-+ 2f uv cot 6dz=—va (=), wee ee (LO) 
since w=v=0 when z=6 and w=0 when z=0. 
We may now follow Karman’s approximate method of solution of the 
rotating-disk problem by putting 


2 
u=A[Z—323-+ 278] — 2 ain 9 cos 8 [$Z2—Z2-4 474), 


(11) 
v= Qa sin 6 [1—$Z+423], 
where Z=z/d ; these expressions satisfy the conditions 
07u 2°52 Ou 
0 aya sin @ cos 6 at Z=0 ; =O; 57 = at Z=1, | 
BO aio ow 0 at Z=0; =0 Ue att 
v= Qa sin 0, ir v=0, 57 =0 at Z= 
(12) 


the conditions on the second derivatives at z=0 being obtained direct 
from equations (1) and (2). The significance of A is that it is 3(0w/dz),. 
The quantities A and 6 are unknown and are to be determined by (9) 
and (10). Let us for convenience write A=(y 2762a/v) sin @ cos 0, d= (v/ 2)12A, 
so that x, A are the new unknowns. 

In the solution of the rotating-disk problem it is certainly true that 


[odes [wae era, Pulremerer eee) 
We may therefore make a tentative sagonueae to the present problem by 
neglecting the terms containing tp uw? dz on the left-hand side of (9). 
It then follows from (9) that 
x= | Sal7 oA 7i\3 dz eend Mes Wace ah (14) 
so that x= 0-2357. SE Oe tet e ke koe c.f (LO) 
It then follows from (10) that 
(5 +2 cot a) (0:00864 sin? 6 cos @ == =sin@, . . (16) 
so that 
3 sin 8 cos gue + (4 cos? é—sin? 0)\4*=173-6, . . . (17) 
and therefore 
& [sins 6 cost? 6 A4]= 231-5 sin? 6 cost? 6. . . . (18) 


Hence ! 


P21. 6 
er ae sin!’ 9 cos! 6 dé. mie (19) 


sin163 9 cos4!3 6 J 4 
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For small @ this solution should agree with the rotating-disk solution. 
From (19) we see that At=43-4 when @ is small, so that 
du Q3/2a 6 ov , 23208 
A=2:57, (=),-08 ae (5), 9 8a 
Cochran (1934), using this type of approximate solution for the rotating 
disk but retaining the terms in w?, gives the corresponding numerical 
coefficients as 2-79, 0:54 and —0-54 respectively. In the solution he 
obtained by numerical integration of the full equations the first of these 
coefficients is infinite, of course, and the other two are 0-51 and —0-62 
respectively. 
Returning now to the rotating sphere, the solution of (19) for A is given 
in Table I. The value y=0-2357 given in (13) leads to 


2 
DS Song = M10 23572 —0-52?-+ 0-292978— 0-0286284} . (20) 
and yg = 1-1 5240-528, Se ee ee 


where A, is the value 2-57 taken by A when 6=0. 


TaBLeE I, 


6 Umax Umax | i 
(degrees) aQ GOR (S2v)1? ook 


0 
10 
30 
50 
70 
80 
90 


0-000 0-000 
0-039 0-174 
0-110 0-500 
0-127 0-766 
0-180 0-940 
0-161 0-985 

0 1-000 


SSO9 


bo bo bo bo 
me bo bo bo 


TUR yw hy by by 
Sworn aan 
COON O~I 
eee 
om 
Serie] 


The right-hand sides of (20) and (21) are both functions of Z only and 
are given in Table IT. It will be seen that the approximation we have 
obtained is such that w and v retain the same form throughout the range 
0 to 7/2 for 0, the actual scale and amplitude alone varying with 0. 


TABLE II. 


z : WX v y urs v 
A2aQ sin 8 cos @ | aQ sin 6 AaQ sin 6 cos 6 aQ sin @ 
Ne 90-0000 1-0000 0-4 00-2133 0:4320 
5 0:1246 0-8505 0-6 0-1386 0-2080 
: 00-1944 0-7040 0-8 0-04.50 0:0056 


0-3 0-2205 0-5635 1-0 0-0000 0-0000 
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The maximum values of the right-hand sides of (20) and (21) considered 
as functions of Z are 0-22 and 1 respectively. The variation of the 
associated values w,,,, and v,,,, are included in Table I. The ratio 
Umax/Umax 18 also included in Table I. 

It follows from (20) that 


1 or) 
Steal udz=0-0184A3 sind cos@. =... . (22) 
0 


Furthermore, from (3) 


; IP Gh Wr oe) 
w= —a| gq] mdeteote | ud: | ; 


where wy is the value of w at the edge of the boundary layer. Hence, 
from (22), 


Ws=— 90-0184 (Qv)!? E (A? sin 6 cos 6)-+- A? cos? a| 
A (23) 


=0-0184 (Qy)12 E cos? 6—sin? @)A+ = sin 6 cos 6 oe ; | 


so that by equation (17) 
: 2 
Ws=—0-0184(Qv)H2 Jae aa j “| 
Values of w; are included in Table I. It will be noticed that the solution 
gives an inflow throughout the whole range. 
It follows from (19) that A (and hence 8) behaves like c cos-13 6 (where 
c is constant) in the vicinity of 0=7/2. In such circumstances, neglect 
of the first term in (9) is no longer justified and we must expect our 
solution to break down in this vicinity. However, in view of the nature 
of the singularity, we should expect the breakdown to be limited to the 
immediate vicinity of the equator where, anyway, the interaction with 
the flow from the other hemisphere will be important. Since this 
interaction is probably largely determined by the flux through the 
boundary layers, which our solution will give with little error, it is 
worth while examining our results in this neighbourhood further. It 
follows from (20) and (21) that, as 6>7/2, 


(24) 


v 
u— 0* and—~—> 1 


ai 


throughout the layer. The limiting values of both (du/dz)) and (dv/dz)y 
are zero. It appears therefore that the solution presents us with something 
akin to separation at the equator. 


) 
Although u— 0, | u dz remains finite and equal to 0-497 (Qv)!?, so that 
0 


* This result appears to contradict earlier remarks but it should be remembered 
that the value of u is unlikely to be correct in this region although the non-zero 


6 
value of | u dz (see below) may be more reliable. 
0 
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there is a net flux through the equatorial plane. This is to be expected 
from our solution, since as we have seen there is an inflow into the boundary 


layer and it can readily be verified that 
m2 
2a |. ud | = —2na? | Ws, sin 6 dé Bere Pk eT 
0=n/2 


0 0 
as continuity requires. 

This flux through the equatorial plane cannot of course take place 
from symmetrical considerations and must be annihilated by the flow 
from the other hemisphere. As indicated earlier on general grounds, a 
complete solution must take account of this interaction which one might 
reasonably expect to be confined to the vicinity of the equatorial plane 
and to give a mass outflow of the order of the total flux from the two 
hemispheres. 

On the basis of this solution it appears likely therefore that inflow 
will occur over a large part of the surface, the outflow necessary to 
maintain continuity being confined to the vicinity of the equatorial plane. 

Moreover, the general equations (1) and (2) confirm that inflow probably 
occurs at the edge of the boundary layer. At the outside of the layer these 
equations may be written approximately 


ou 07u 

Tee aD 2 
6 Op” Orr? (26) 

ov 0?v 
Ura ae 


where w, is the outflow at the edge of the layer. Now it is most unlikely 
that du/dr and dv/dr are other than negative and that 0?u/dr2 and 02v/dr? 
are other than positive in the outer parts of the layer and hence w, must 
be negative—that is there must be inflow. 


ALP PRN Dex, 


Another example of an essentially three-dimensional boundary layer 
is provided by the boundary layer along a thin hollow semi-infinite circular 
cylinder which is made to rotate with constant angular velocity w about 
its axis and placed in a uniform stream moving with velocity U parallel 
to the axis. 

It has been shown by Howarth (1951) amongst others that the curvature 
terms disappear from the equations of flow for the boundary layer along such 
a cylinder, so that we have the following equations : 


Ou ou 02u 
a UP Ming Arm Tels ch ky I eb) 
Ov os 02v 
Un ee a Sil ornate ees EL 


Ou ow | 
Fate Be Ae Oe eee ol SE ema) 
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where x IS measured parallel to the axis from the leading edge, y around 
a right section and z normal to the surface. 


Equation (28) is the boundary layer equation for flow along a flat 
plate, so that we have 


u=4U0F"’(Z), 
where Z=1 (=) Zz 
~ \va 
and BFR” —0, 


with F=F’=0 at Z=0 and F’=2 as Z-> 00 *. 
It follows from equation (30) that 


1 /fvU\12 


and then from (21) that 
v=twa(2—F’), 


since v=wa when z=0. 
The surface traction in this case has components 


3\ 1/2 
Tog Or Sd2 P (=) ? 


U\12 
T y= — 9:332 pwa (=) ° 


The drag on the cylinder is thus unaffected by its rotation and the 
frictional moment M which has to be overcome in order to maintain the 
rotation is 1-3282 mpwa*(vlU)!? for a length J of the cylinder ; the moment 
coefficient Cy, is given by 


Cy= ma2lU2 p U 


As a result of this nature can clearly be applied to determine the rate 
at which the spin would be lost if the rotation were not artificially 
maintained, provided the rate of loss were sufficiently slow for the 
aerodynamic problem to be considered steady at all stages. In such 
guise it is not without interest ballistically. 
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ABSTRACT. 


An attempt is made to derive theoretically the Simon formula for the 
variation of the melting point of a solid with pressure. The model used is 
similar to that of Lennard Jones and Devonshire, but the interpretation of 
the results is different. A formula of the right type results, but a detailed 
examination shows that the melting temperatures are much too high, and 
the reasons for the discrepancy are discussed. The theory predicts the 
absence of a solid-fluid critical point, and this is in disagreement with 
recent general work of Miinster. 


§ 1. InrRoDUCTION. 


Tue form of the solid-fluid transition curve at high pressures has been the 
subject of considerable discussion for some time. Particular attention has 
been fixed on the question whether a solid-fluid critical point exists above 
which there is continuity of state, as in the liquid-vapour transition. 
There were quite lively discussions on this matter between those who 
favoured the existence of such a critical point and those who were 
convinced that it could not exist. All these discussions were, however, 
based on very insufficient experimental data, and on the extrapolation of 
interpolation formule. The later and more complete experimental work 
of Bridgman showed no indication of the existence of such a critical point, 
and Bridgman considers that the evidence is fairly conclusively against it 
(see Bridgman 1949, p. 201, for a detailed discussion). On the other hand it 
has been pointed out by Simon (1937) that thermodynamic reasoning 
alone cannot suffice to decide either way, and also that there are no 
obvious reasons which exclude the one or the other possibility. 


Some years ago Simon proposed, on semi-empirical grounds, the following 
formula for the melting curve : 


Pa (mr) -1. Sn Cae ee 


(Here T, is the normal melting temperature, p, the internal pressure and c 
is a constant.) This formula is the first one to connect the melting curve 
with properties of the substance. It fits all the available experimental 


data with reasonable accuracy, the values of the internal pressure are 
a See aE LS en re ee ee 


* Communicated by Professor F. Simon, F.R.S. 
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found to be in good agreement with those calculated by van der Waals’ 
equation, and the value of ¢ turns out to be round about 2 for most sub- 
stances (about 4 for metals). Simon argued that if one wants to penetrate 
relatively to very high temperatures and pressures one has to carry out 
experiments on substances with low internal pressures, and he therefore 
used the substance with the lowest internal pressure, helium, as a model 
substance. Recent work on helium (Simon et al. 1950) has shown that 
the formula fits the experimental data with values of T/T, ranging up to 
50°.* It is of interest to examine the possibility of deriving theoretically 
a formula which fits experimentally over such a wide range. 

In 1939 Lennard Jones and Devonshire put forward a theory of melting 
analogous to the theory of order-disorder transitions in alloys. At 
sufficiently low temperatures they obtained isothermals of van der Waals 
shape with two stationary values, which they interpreted in the usual way 
as a first-order phase transition. At higher temperatures the two stationary 
values disappeared, thus indicating the existence of a critical point. 
However, Lennard Jones and Devonshire pointed out that they considered 
the model too crude to make any reliable predictions. This model was 
refined slightly by Blaisse (1946) without substantially altering the 
conclusions. On using an adaptation of the Born theory of melting, Blaisse 
was able to derive a formula of type (1), but the value of c obtained was 
less than 1. 

In the present paper we shall adopt a model similar to that of Lennard 
Jones and Devonshire. Our interpretation of the consequences of this 
model will differ, however, and we shall derive a formula similar to (1) for 
the melting curve. We shall also indicate the mathematical developments 
required to improve the theory. 


§ 2. FORMULATION OF THE PROBLEM. 


We shall throughout assume classical theory, interactions between 
molecules being given by central forces of limited range. Our results will 
therefore only be applicable to systems in which the binding forces can be 
regarded as central, such as the inert gases, and will not be valid for metals 
in which the bonding is of a different character. The partition function for 
an assembly of classical particles reduces to a product of two factors, a 
kinetic energy factor (27mkT)??/h3 per particle, and a configurational 
factor arising from the potential energy terms. The configurational factor 
can be conveniently specified in the following manner. Let us take an 
arbitrary lattice structure whose spacing is small compared with the 
intermolecular distance in the solid state. Suppose that molecules are 
restricted to occupy positions in the centre of a lattice cell. Between any 
two molecules in different lattice cells there is an interaction energy «,. 
The intermolecular forces can be ignored outside a certain radius, and for 


* The normal melting point of helium would correspond to about 1° abs. if 
quantum effects did not interfere. 
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any finite lattice spacing the «; form a finite set 7. In any possible 
configuration of the assembly let 7; be the number of pairs of molecules 
with interaction energy e;. The energy of this configuration is 


pr N ;€;- 
i=1 
Let W(n,,..- 2% .++,) be the number of possible configurations of the 
assembly in which there are n, pairs of energy «),.....- n, pairs of 
energy e,. The configurational partition function is then 
Lo Dy W(n, erareje nN; seis Ny) exp [— D4; n€,[kT]. : (2) 
Fights t=] 


The expression (2) is only exact in the limiting case of infinitely small 
lattice spacing. But any finite spacing which contains a sufficient number 
of points to give a fair description of the intermolecular energy curve 
should provide a reasonable approximation. 

The lattice used above is fixed and artificial and, provided that it is 
sufficiently fine, any given crystal structure can be fitted in to it. At lowest 
temperatures the system will choose the configuration of minimum energy, 
and if the law of force is of Lennard Jones type (A/r"—p/r"), this corre- 
sponds to the face-centred cubic structure (Lennard Jones and Ingham 
1925). An increase in pressure is represented by an increase in the density 
of particles, and the formula (2) also embraces the possibility of poly- 
morphic transitions. In fact the description’ is too general for our 
purposes, since we are concerned with the disordering of a particular 
crystal structure, and we shall therefore modifiy the model. 

We consider the number of particles as fixed, but the lattice spacing as 
variable, an increase in pressure corresponding to a decrease in lattice 
spacing ; we also confine our attention to a particular cyrstal structure 
(fig. 1). At constant volume our problem is now formulated as a 
generalized order—disorder transition in which many neighbour interactions 
must be taken into account ; the ratio of occupied to vacant cells is very 
much less than 4. It should be observed that thermal vibrations are 
automatically taken into account in this formulation, 


§ 3. APPROXIMATE MopDEL 


If we could solve the above problem exactly we should find that as the 
temperature was raised from zero a point would be reached at which the 
partition function had a singularity. This singularity would correspond 
to a first-order transition, and would be associated with the disappearance 
of long-range order; it could thus be interpreted as the melting point. 
Unfortunately the only reliable information we have concerning order— 
disorder transitions is in regard to systems containing equal ratios of 
constituents and in which only nearest neighbour interactions are taken 
mto account. We are therefore forced to use a drastic approximation in 
our model, and we must then decide whether any information provided by 
this approximation can be considered reliable. Following Lennard Jones. 
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and Devonshire, we take account of small displacements from equilibrium 
by assuming a harmonic vibration spectrum, the kinetic energy terms now 
being included in the partition function ; the gradual disordering of the 
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crystal structure is allowed for by assuming that the molecules can jump 
into a single interstitial position, as shown in fig. 2. An essential 
physical difference between the solid and fluid states is the presence or 
absence of long-range order, and the model can take this into account. 
We can thus divide the partition function into two components, 


F=F,,+F o ’ . - . . . . . . . (3) 


F,,,, being vibrational and F,, being given by order-—disorder theory. 

In their discussion of the problem Lennard Jones and Devonshire use a 
refined calculation for F,,,, averaging the inter-molecular energy over a 
cell to which a molecular is confined by its neighbours. We do not con- 
sider that the crudeness of the model warrants any such detailed calcula- 
tion, and we shall use a simple Debye model, so that 


F,,=—3N&T In KT/hv, ee eee 


v being a mean frequency. In fact we do not consider that the present 
model is capable of yielding any detailed information on the nature of the 
singularity, and hence on matters such as the magnitude of the latent heat. 
Nevertheless, we do believe that it is possible to derive a reasonable idea of 
the position of the singularity. This assumption is based on the experience 
of order—disorder transitions for which an exact solution is available in two 
dimensions (Onsager 1944, Wannier 1945). It is there evident that crude 
approximations are quite incorrect in regard to the nature of the singularity, 
but provide a fair approximation to its position. The exact solution shows 
that the singularity is associated with the disappearance of long-range 
order, so that we can speak with precision of an ordered and a disordered 
phase meeting at the singularity. This is essentially different from a 
liquid and vapour above the critical point, where there is complete 
continuity, and a single analytic function for the partition function. We 
shall interpret the above singularity as corresponding to the melting point. 

For F,. we shall use the Bethe approximation with a coordination 
number equal to 12 (Rushbrooke 1949, p. 306). Then if « is the repulsive 
energy between a molecule in equilibrium and an interstitial molecule, the 
temperature of the singularity is given by 


1 
kT .=ae | a= 2In 12/10 =287 | . woe 6 (5) 
To determine the corresponding pressure we must take account of the 
variation of « and » with volume, using the formula p=—dF/dv. For > 


we shall use Gruneisen’s assumption >=K/v’, so that for the vibrational 
term 
deg ere 3NkT op = 38NKTy 
Pub Ov aS re 7 On == >a tie ° oT et tows (6) 


For the intermolecular energy we shall assume a formula of Lennard Jones 
type, the attractive part of which can be neglected in the interstitial 


position. Thus 
7r\~—n v\—nl3 
e= Aey (=) =Ae€, (=) ’ . . ° . ° . (7) 


Melting Curve at High Pressures 1321 


where v is the equilibrium volume of the substance and A is a constant. 
(If the interstitial position is at a distance from equilibrium equal to 1/4/2 
of the spacing between atoms, the value of A is 2”/2.) Typical values of the 
constants for various substances are given in Fowler and Guggenheim 
(1939, p. 285). 

Hence we have 


co Ais de Nbn Acg /v \-(™/)-1 
Por aire Ov =—Nb 5 = (2) <iet's eu 's (8) 
127.10 
a —— 
(s3.229 203) 

Combining (6) and (8), and substituting from (7) in (5), we find 

De AF Uo (n/3)+1 nb Le Ne 

as Uo n/3 
and — =Aa (=2) Fen eae Oe Th A: Ae (10) 


The pressure w is determined by the lattice energy at the absolute zero 
and is related to the internal pressure, and 7 is related to the ordinary 
melting temperature. Therefore 


Dee ae sy nb/sa 
ae —x (=) B= (Aa)3/" . . . ° . . (11) 


At sufficiently high pressures this is a formula similar to (1) with c=1+3/n. 
ti ——9, C= 1-3a-and i n=12, c==1-25. 


oa 


§ 4. Discussion oF RESULT. 


It was pointed out by Simon (1937) that the existence of a universal 
formula of type (1) implied a law of corresponding states for melting. 
Lennard Jones and Devonshire (1939) observe that as long as the inter- 
molecular energy curves of various substances can be adequately 
represented by an expression with two free parameters (such as (7) with 
fixed n), we will beled to a law of corresponding states. However, Bridgman 
is opposed to the idea of corresponding states for melting (1931, p. 207), 
since it is clear experimentally that latent heats and volume changes on 
melting cannot be fitted into such a law. Nevertheless, a law of correspond- 
ing states should only be regarded as a rough approximation ; it may well 
be that a feature such as the melting point depends only in an average 
way on the shape of the intermolecular energy curve, whereas features such 
as latent heats and volume changes are much more sensitive to details. 
In this connection it is worth noting that the formula (11) is not very 
dependent on the coordination number used in the model. 

From the experimental point of view Simon (1951) has shown in a recent 
article that the data for the melting curve of a number of substances fit 
a law of corresponding states with considerable accuracy. 
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When we attempt a detailed comparison between (11) and experiment we 
find two serious discrepancies. The value of c is too low, and the values of Te 
required to fit formula (1) are about three times too large. If we re-examine 
our approximate model it is not difficult to trace the source of the dis- 
crepancy. The single interstitial position we have taken corresponds toa 
very high energy ; for n=12 the repulsive energy in this position is 64 
times its value in equilibrium. Clearly a violent degree of thermal 
agitation and disorder would be required before this position would 
become occupied with any frequency, and this is far more than is required 
for melting. Interstitial positions much closer to the equilibrium position 
are sufficient to provide the disorder required for melting, and these would 
correspond to much lower temperatures. They are also situated on a 
flatter region of the intermolecular energy curve, and should therefore give 
rise to a larger value of c, since (11) indicates that the value of ¢ varies 
inversely as the steepness of the curve. Confirmation of this interpre- 
tation is provided by X-ray diffraction experiments which show no very 
marked change in the shape of the first maximum between solid and 
fluid states. 

It may also be noted that the Lindemann melting formula, which is well 
satisfied experimentally, assumes that melting occurs when the average 
vibration amplitude is about one-tenth of the minimum distance between 
atoms in the solid lattice (Griineisen 1926). Our choice of the average 
interstitial position has been arbitrary, and if we were to take an interstitial 
position at this distance from equilibrium, our results would be in better 
agreement with experiment. 

Unfortunately, although one can account for the disagreement 
qualitatively, it is difficult to see how to take adequate account of it 
theoretically without tackling the more complicated problem of generalized 
order—disorder transitions. We could start with a problem of the type 
indicated in fig. 3 in which nearest and second neighbours both have a 
repulsive energy of interaction, so that the lowest energy state is as 
indicated. Then we should try to introduce more intervening repulsive 
interactions. Our resulting critical temperature would only be expected to 
tend to a non-zero value as the number of interactions increased 
indefinitely for a three-dimensional model, since Peierls (1936) has shown 
that thermal vibrations on their own are sufficient to destroy long-range 
order in two dimensions. 


§5. THE EXIsTencr or a Crrroan Porn. 


In a recent paper Miinster (1951) attempts to show very generally that 
all first-order phase transitions disappear at sufficiently high temperatures, 
and hence that these transitions end in a critical point. This conclusion 
contradicts the result of the present paper, and we believe it to be 
incorrect. Miinster’s argument is rather intricate, but Mr. R. O. Davies 
has pointed out to the writer that an error arises in the passage from a 
finite canonical ensemble of N particles to a grand canonical ensemble 
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which provides in principle for an indefinitely large N. The finite limit to 
which dE ,/0@ tends as @-> oo may well depend on N, and its behaviour for 
large N must be investigated. 
This point may be aptly illustrated by a particular physical example. 
Consider an order—disorder transition in a two-dimensional quadratic 
lattice with equal numbers of the two constituents. An exact solution is 
available (Onsager 1944), and although the transition is not first order, the 
specific heat is infinite at the Curie point, the fluctuations are large, and 
Miinster’s argument should apply. For any given number N the Curie 
point is given by 
SLA rel DMO Rta thie Sa) ghee, 12) 
where « depends on the intermolecular spacing, and hence on «/N (« being 


the area of the assembly). As N increases so does «, and hence T,. If 
there is no upper bound to « there is likewise no upper bound to T,. 


Big. 3, 


It seems clear, therefore, that a critical point could only arise in general if 
there were an upper limit to the intermolecular energy as the intermolecular 
distance decreases to zero. At the present moment there is no indication, 
theoretical or experimental, of the existence of such a limit. Theoretically 
the repulsive energy is of the form 

(7 extn ==71 0) Am oteyei te ieee pel 3) 
where p is a constant and R(r) is a polynominal containing small inverse 
powers of r (Slater 1928). Even if R(r) is replaced by a constant, P, the 
numerical values of P given in Fowler and Guggenheim, p. 298, would 
correspond, from equation (5), to melting temperatures 10° times the normal 


melting temperature. 
We may sum up by saying that we have indicated the possibility of 


deriving, for a classical model, a formula of type (1) from the melting curve. 
4U2 
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The parameters occurring in this formula can be related to the normal 
melting temperature and internal pressure. It should be emphasized, 
however, that the Simon formula (1) is of much greater general validity 
than the present theory since it applies even when quantum effects are 
present (Simon 1951). 
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CXXXI. CORRESPONDENCE. 


The Influence of Huchange Energy on the Specific Heat of 
Free Electrons in Metals. 


By A. B. Liprarp, 
King’s College, London *. 


[Received August 15, 1951.] 


In a recent paper Wohlfarth (1950) has published calculations of the 
effect of the exchange energy on the thermal properties of free electrons. 
His starting point is an expression (due to Koppe (1947)) for the free 
energy, F, in terms of the distribution function f(k). (His equation (2.1) 
is, in fact, only $F and should be multiplied by 2.) The first two terms 
in this expression are respectively the kinetic energy and the exchange 
energy of the electrons ; together they make up the internal energy, E. 
The third term is the entropy, 8S, multiplied by —T. As an approximation 
it is assumed that S is the same function of f(k) as it would be in the 
absence of exchange interactions. We should expect this to be a good 
approximation when the exchange energy is small compared with the 
kinetic energy, simply because it is correct when the exchange energy is 
actually zero. An integral equation for f(k) may be obtained by putting 


the variation 5F=0, namely 
dk,—C=kT In (Si”). mabe 


he ae 
f(hx) 


82m ki — nk, an CALE Ny cay 


¢ is a parameter introduced by the condition that the total number of 
electrons, N, is a constant. It is equal to the chemical potential 0F/0N. 
Owing to the difficulty of solving equation (1) its derivation by putting 
5F=0 is made the basis of a variation treatment. As Wohlfarth takes 
a rather crude variation function f(k) (equation 2.3), it seems desirable 
to obtain some confirmation that the results of its use are not seriously 
in error. 
Firstly we note the possibility, in principle, of solving (1) by a method 
of successive approximations, an nth approximation for f(£) being inserted 
on the L.H.S. to give an (n+1)th approximation on the right. We could 
therefore consider Wohlfarth’s variation function as a starting point in 
this iterative process. This we have done, and found that for low 
temperatures the next approximation was, in its variation with k, 
graphically indistinguishable from a simple Fermi-Dirac function, in 
which the temperature T is replaced by a parameter 7, say. 7 Was 
occ le tt anc ill all pete a eB aoe 2a a el TNE Sid, SEA: AOE 2 
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several times smaller than T, in general agreement with Wohlfarth’s 
result that the specific heat was greatly decreased by the exchange 


energy. ieee 
This means that probably a very good variational distribution function 


would be : & 
s0)= [expge (Sa, #—t) +1] kee 


The parameter ¢ may be found as usual by the condition that the total 
number of electrons is a constant, N. The second parameter 7 will be 
determined by putting dF/@7=0: when the exchange energy is not 
included, this method, of course, yields t=T, as is necessary. Now ~ 
Yokota (1949), in the course of a paper on the Thomas—Fermi—Dirac 
equation, has evaluated the exchange term in the expression for F, 
using the (incorrect) Fermi—Dirac distribution function. It is thus only 
necessary to substitute 7 for T in his result (6) to obtain F(z), since the 
kinetic energy and entropy expressions are well known. (There is a 
minor error in Yokota’s equation (6) which is corrected here.) In this 
way we obtain 


B38 (2) sat Sef 2 y )] 
aa AON Ws Tet ee Pe i he ait 
sen s(i+ =) pl 1+ 5 (met 3, oe 


Sur ARE GE). oe ee 
€9 = 
in which 
ake A h? /(3N\2/8 be 3N\¥8 
one “0o— 8m \nV > Sjeae a7V 


and y is a numerical constant. (In terms of y,, yg, and y3, defined in 
Yokota’s paper, 
y=3—4 In 24 12y,/7?—y.+2y, : 
the 3 appears incorrectly as a 2 in his paper.) We have evaluated y 
numerically and found it to be 0-648. The quantity «, is just the Fermi 
energy, and «;,/eg gives a convenient measure of the exchange energy. 
By minimizing F with respect to x (7. e. with respect to 7) we obtain as 
the equation determining a as a function of T 
akT e,[ 2 x y 1 
=—gw— 2] — pase pte a 
a =2 | Sin ( b+ Int 5) |, sents . (4) 
These two equations (3) and (4) are to be compared with Wohlfarth’s 
equations (2-6) and (2.7), namely 


F 3 5a 3 €; 2 13 
Ne, = (14+ 7) = aae Ee (in x—In 2— 3) | 
3RT x 
*9? “<i, Pt oos She 3 (5) 
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Apart from the substitution of 3 for 7 in the entropy terms, these 
two equations differ from (4) and (5) only in the replacement of 


i i 
(Inw+ 5—3) =1-32 by (in 2+ 3) =1-28, 1.e. a difference of 4 per 


cent. Finally, Wohlfarth’s equation (2-8) for the specific heat, C, is 
replaced by 


l ¢; e; {1 1 =1 
CO {i 5-gine— 2 (5 +2 Jina) > ee (7) 
where C® is the Sommerfeld specific heat }72NR(RT/e)). This again is 
equivalent’ to the replacement of In2+7/12 by Ina+4—y/2. This 
reaffirms the result that the specific heat is greatly decreased by the 
inclusion of the exchange energy. For it is clear from (4) and (6) that 
as T > 0, x + 0 and from (7) that as x > 0, C/C® > 0. 

This close agreement between the results obtained by using the first 
and second order variation functions is strong evidence for the validity 
of the variation method. This is confirmed by a comparison of the 
approximate and exact values of the thermodynamic functions for a 
perfect Fermi-Diracgas without exchange interactions (published elsewhere, 
Lidiard 1951). 


In conclusion I should like to thank Isaaki Yokota for supplying the 
derivation of his equation (6) which has been used here. 
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A Dislocation Reaction in the Face-Centred Cubic Lattice. 


By W. M. Lomer, 
Cavendish Laboratory, Cambridge*. 


[Received July 15, 1951.] 


Tuer marked difference in the initial rates of strain-hardening of hexagonal 
and cubic crystals has often been ascribed to some form of interference 
between slip on non-parallel planes in the cubic material. It is shown 
in this paper that dislocations on two intersecting slip planes can attract 
one another and combine to form a “ sessile ’”’ dislocation (Frank 1949), 
which may act as a considerable immobile nucleus of internal stress, and 


so cause the crystal to harden. 
Reece ee <a ee ee eee 
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Consider two dislocations on slip planes inclined at an angle @ as shown 
in fig. 1. (Only dislocation lines parallel to the line of intersection will 
interact strongly, so the problem is effectively two dimensional.) The 
force per unit length exerted on dislocation 2 by dislocation 1 is 697, 
where 6, is the Burgers vector of number 2, and 7 is the shear stress due 


Fig. 1. 


Coordinate system used for calculations. 


to number | resolved first into the slip plane of 2, and then into the slip 
direction. A full analysis shows that the force can be derived from the 
following energy functions :— 

pbb, log r 


1 and 2 both screw type E= 9 
7 


} » poids pq sin* 0 
1 and 2 both edge type E= Sale) {0s 6 log r+ er a 
1; edge and 2; screw, or vice versa, E=0. 
where .=modulus of rigidity, =Poisson’s ratio, 
b,, b,=Burgers vectors, 
p,q, 7 as diagram, 

In the close-packed cubic lattice the planes may be for instance (111) 
and (111), intersecting in [110] and making an angle cos! 1/3 with each 
other. The only dislocation lines considered are therefore those parallel 
to [110]. Possible displacement vectors (excluding “ half-dislocations ” 
(Heidenreich and Shockley 1948) for reasons explained later) are in the 
directions [110]; [101]; [011], in the first plane, and [110]; L101 
[011], in the second. The first one of each of these sets is a pure screw, 
and of no particular interest. The reaction to be studied in detail is 
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that between [101] and [011] which can add to give [110], a single 
dislocation of the same strength as each of the two combining ones. 
This resulting dislocation line. is. parallel to [110] and its displacement 
vector is [110]; these two directions are not contained in the same 
close-packed plane, and so the dislocation is sessile. 

The displacement vectors may be resolved into screw and edge 
components as follows, with a=length of unit cell edge :— 


a{lOl]=Ja ((1T0]+[113)) 
4a[011]=4a(—[110]+[112]) 


Fig. 2. 
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Interaction energy contour map for reaction [10I]+[011]—> [110]; energy 


2 


contour interval= i ergs per cm.—0-85 e.V. per atom plane in copper. 
TT 


giving the screw components a magnitude a/2,/2 and the edges / 3a/2+/2. 
With these values, and taking v=1/3 we get, adding all the energy 


contributions 
K= so 5 { 5 log Hee ae 


It is convenient to represent the energy as a function of the oblique 
coordinates p and q by the contour map of fig. 2, where the energy value 
for dislocation 1 at q and dislocation 2 at p is plotted at the point p, q. 


* gy, 6 
= 5 [101] i is a vector with components 5 5? 0, — 5 and so on. 
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It may be presumed that the representative point of the system will 
move down the line of greatest slope on the diagram, so that both 
dislocations will in practice move simultaneously. A typical “ trajectory ” 
for the system has been sketched in. 

There is also an attractive force between dislocations with Burgers 
vectors [101] and [101]. It is small and depends on the fact that the 
self-energy of a screw dislocation per unit length is less than that for an 
edge with the same displacement vector. The energy function for this 
pair is 


This is mapped in fig. 3; the trajectory shows that at first the dislocations 
adjust their angular relationship, and only when p/q is approximately —1 
do they both approach the centre. 


Fig. 3. 
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Interaction energy contour map for reaction [101]+[101] > [200]; contour 
interval as fig. 2. 


The details of this picture probably depend on the assumptions made 
such as the infinitely long, straight dislocations, and on neglecting the 
effect of the dissociation into half dislocations. Since it is believed that 
the halves will only separate by some 10 A. this cannot have a large 
effect on interactions at a considerable distance, but might play some 
part in the details of the actual process of coalescence of the dislocations 
The production of sessile dislocations by the interaction of mobile ones, 
and not only by stacking faults in the growth process, should almost 
certainly have some effect on the rate of hardening, but it is not yet 
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clear to what extent it will be important, or what the exact mechanism 
of the interaction of the sessile dislocations with those still mobile 
would be. 

Recent experimental work has shown that hardening is very closely 
associated with the occurrence of deformation bands (Honeycombe, 
in the press), and that the latter occur more frequently in crystals with 
several slip planes. This may imply that the effect dealt with in this 
paper is not important, although it is perhaps possible that sessile 
dislocations produced in this way help the formation and stabilization 
of the deformation bands. 
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Lattice Defects in Silver Bromide at Room Temperature. 


By H. D. Kerrx and J. W. MircHet, 
University of Bristol*. 


| [Received September 7, 1951.] 


As a result of the experiments of Tubandt and Eggert (1920), Wagner 
and Beyer (1936), Berry (1951) and others, it has now been established 
fairly conclusively that, at temperatures above 200°C., the lattice 
disorder in silver bromide is predominantly of the Frenkel type. Con- 
cerning the nature of the lattice disorder in this material at lower 
temperatures, however, little direct experimental evidence is available, 
and an attempt has therefore been made to determine the ‘“‘ soundness ”’ 
of pure silver bromide at room temperature from newly-measured values 
of its lattice-parameter and density. 

The experiment consisted of comparing the molecular weight of silver 
bromide as determined by X-rays, that is, as calculated by means of the 
formula (Straumanis 1949) 

3 

M,= Lee 

n 
(where k is a factor introduced to correct for an error in the value of the 
molecular weight of calcite used by Siegbahn, N, is the effective Avogadro's 
number, a is the lattice-parameter in kX, n is the number of molecules in 
the unit cell of the ideal lattice and d is the density in gm./c.c.), with 
the value, M, of the molecular weight which is known from accurate 


* Communicated by the Authors. 


1332 Correspondence 


chemical data. If D=(M,—M) is negative and greater than the 
experimental error, then it can be concluded that the silver bromide 
contains vacancies which may be located either at grain or mosaic 
boundaries or within the lattice (Schottky disorder and/or aggregates of 
vacant lattice sites); the experiment does not, however, enable the 
presence of Frenkel disorder to be detected since interstitial ions and 
vacant lattice sites in equal concentrations do not affect the macroscopic 
soundness of a crystal. 

The lattice-parameter of pure silver bromide was measured at 25° C. 
by means of the X-ray powder method, the results of two determinations 
carried out with a directly calibrated 19 cm. “ Unicam ” camera and with 
a 64mm. Straumanis camera agreeing to within 1 part in 4x 10* 
(a=5-7633 kX). The densities of large well-annealed polycrystalline 
specimens of silver bromide were determined by the displacement of 
doubly-distilled water using the method of Snoek (1950) and the results 
obtained with three such specimens agreed to within 1 part in 2x 10+ 
(d=6-4799 gm./c.c. at 18° C.). For the purpose of correcting the value 
of the lattice-parameter for temperature differences, the value of the 
thermal expansion coefficient given by Strelkow (1937) was used. 

The averaged result of six determinations of D gave the result 
(D/M)=+8-6x 10-> but, due to experimental errors and uncertainties 
in certain fundamental constants, this may be in error to the extent of 
+1:6 10-4. The maximum possible number of vacant lattice sites in 
pure silver bromide at room temperature represents, therefore, a fraction 
(1-6— 0-86) x 10-4=7-4 x 10-° of the total number of lattice sites in the 
ideal crystal. This is probably an over-estimate and, as allowance has to 
be made for the probably not inconsiderable number of vacancies 
associated with mosaic and grain boundaries, it is therefore unlikely that 
Schottky defects are present in sufficiently large numbers to contribute 
significantly either to ionic conductivity or to the mechanisms of latent 
image formation. The situation may well be different, however, in the 
sensitized silver bromide grains of photographic emulsions (Mitchell 
1951). 
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A Note on the Trace of the Product of Dirac’s Matrices. 


By L. M. Yana, 
The University, Edinburgh*. 


[Received September 7, 1951.] 


Ir is well known that the trace of the product of any odd number of 
Dirac’s matrix y, is zero. For a product of any even number of V8 
the trace can be written down simply in the following manner : 


Ns oe. Vn, Ped Snpne? Saat Tete Lad) 
UU 
2n factors 


where x{...} denotes taking the trace of the matrices in the bracket, 
and the summation sign indicates summing over all possible pairings 
({(2n) !]/[7 ! (2 !)"] in number), preceded by a positive or negative sign 
according to whether an even or odd number of transpositions are 
required to bring (n,n,)... (n,n) to (u,v)... (y, ¢), transposition within 
each pair not counted. 

The proof of the rule consists in the repeated application of the 
relation y,y,+Y,¥,=268,»5 


XV Va = + MnVSH=2V a = YnV Suv — XLV Vala + Mave 
=x (Va- + + YnVe}Sue— XLV 0 = MnV hua + - 


+x a RAQKO Yn} oucl—X{V Va coo YnV Yue 


Since any circular permutation of the factors inside the bracket does not 
alter the value of the trace, one has therefore 


> 


XY Vane VnV=KV a «=  VnV Suv — Xo «+ + VnVeb Ona 
Be ae Cee Yaseen ee (2) 


The sign preceding each term on the right-hand side of (2) is positive 
or negative according as the number of transpositions required to bring 
the array of indices to the initial order on the left-hand side of (2). The 
process can obviously be repeated and the argument holds for each 
subsequent reduction until finally each term contains n—1 Kronecker 0’s 
multiplied by a factor of the typical form x{y,y-}—46,-. Thus (1) is 
proved. 

For practical purposes, the possible pairings and their appropriate 
signs can most easily be obtained by labelling the 2n indices u,v, A... 7, ¢ 
on the circumference of a circle and connecting the pairs by lines drawn 
outside the circle. The sign is determined by counting the number of 
crossings of the lines in a given graph, positive for the even number of 
crossings and minus for odd ones. 
JIE ee ee ee ees 

* Communicated by Professor M. Born. 
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Three-Photon Decay of Positronium. 


By J. M. RADCLIFFE, 
University of Birmingham*. 


[Received August 31, 1951.] 


THE two spin states of positronium in its ground state have markedly 
different life-times against annihilation. Positronium in the singlet state 
has the shorter life-time (r=1-25 x 10-1 sec.) since it can decay into 
two photons. Such a decay from the triplet state is forbidden (see 
e.g. Yang 1950), the simplest decay allowed being one into three photons. 
Three different results of three calculations of the life-time against the 
latter process have been published (Lifshits 1948, Ivanenko and Sokolov 
1948, Ore and Powell 1949). 

As direct experimental determination of the life-time is now possible, 
it was thought worthwhile to repeat the calculation. Ore and Powell’s 
value for the mean life-time, 1-4 10-7 sec., was confirmed. 

A covariant approach, using the methods developed recently by 
Feynman, simplified the work considerably, especially in the sum over 
the polarizations of the photons. The matrix element was evaluated 
using zero-velocity free-particle wave functions, normalized so that the 
density was that of the electron at the positron in the ground state of 
positronium. (The effect of binding is to multiply this matrix element 
by a factor [1-++-0(«)], «=e?/hce). Also it has been shown that three-photon 
decay from the singlet state does not occur, at least in the limit of zero 
relative velocity of the two particles. This means that the matrix 
element which occurs can be conveniently summed over all spin states 
of the particles. 

The writer is indebted to Dr. G. E. Brown for helpful advice, and to 
the Department of Scientific and Industrial Research for the award of 
a maintenance grant, 
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CXXXIT. Notices of New Books and Periodicals received 


Josiah Willard Gibbs. By L. P. WuErnter. (Yale University Press. London : 
Geoffrey Cumberlege.) [Pp. 264.] Price 25s. 


“In the course of the thirty-two years of his teaching somewhat less than 

one hundred students attended Gibbs’ courses.” “Only a handful of these 

students still survive.” ‘‘ In the whole period of his teaching there were but a 

Hine dozen of dissertations which can be said to owe their inspiration to 
ibbs.”’ 

Lynde Phelps Wheeler, from whose book the above quotations are taken, 
is one of the few survivors of the lucky dozen. Having received three separate 
requests to write a biography of Gibbs, he was eventually persuaded 
reluctantly to undertake the difficult task. Mathematicians, physicists, chemists 
and all interested in the history of science will be grateful to him for carrying 
out this labour of love as thoroughly, honestly and conscientiously as humanly 
possible. 

The theme of this book is in the author’s own words “ Gibbs’ scientific work 
was Gibbs’. There is little else, that is live, to tell. In this biography, or any 
other biography of Gibbs, there is a tragedy for him who can read between 
the lines. Whereas Gibbs was honoured and recognized as at least an equal 
by such giants as Maxwell and Rayleigh, and at Yale he was undoubtedly 
well liked and respected by his immediate colleagues, when he died in 1903 
his genius was inadequately appreciated in his own country. No proper care 
was paid at the time of his death to preserving his correspondence. Not until 
a score of years had passed was there any active desire to know more about the 
man and then it was too late. Wheeler is able to tell us something of Gibbs 
as a lecturer, but scarcely an anecdote of his private life. 

An epitaph, similar to Wren’s, would seem apposite to Gibbs: “Si 
momentum requiris, opera lege mea ”’. 
HE. ASG. 


An Introduction to the Theory of Control in Mechanical Engineering. By R. H. 
MacMititan. (Cambridge: University Press, 1951.) 


Fundamentals of Automatic Control. By G. H. Farrinapon. (Chapman and 
Hall, London, 1951.) 


THE present applications of automatic control, or as it has recently been termed 
“* feed-back engineering ”’ are by far too numerous and wide-spread to mention. 
The study of the theory of control, which is now well formulated, is of major 
interest in almost all branches of engineering and of more than passing interest 
in the fields of economics and physiology, since the recent work of Wiener, 
Tustin, et al., have shown the dependence of many natural phenomena upon 
the same laws and equations applicable to mechanical systems. 

The foundations of the basic control theory can be traced to a classical paper — 
by Clerk Maxwell, “‘ On governors” (Proc. Roy. Soc., 16, 270, 1868), but the 
modern development dates essentially from the application of the Nyquist 
method of harmonic response analysis, originally developed for use with 
electronic circuits. There is now a very considerable literature on the subject, 
but no British text-book dealing exhaustively with the subject has previously 
appeared. It is thus gratifying to note the almost simultaneous appearance of 
the above two books. , 

The former, by MacMillan, claims only to be an introduction to the subject, 
but correlates excellently the very diverse fields of application which respond 
to the same basic theory and to identical methods of analysis, e.g. regulation, 
position control servo-mechanisms, and process control. The early chapters 
are mainly descriptive, leading to the concept of transfer functions and system 
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ions. The human operator and non-linear effects are discussed only 
Sera but a very full account is presented of the methods of transient 
and harmonic analysis by both analytical (Laplace) and. graphical (Nyquist 
diagram) methods. An extensive bibliography is provided and the book 
should be intelligible to any engineer or physicist and of particular value to 
the student in this subject. 

MacMillan uses throughout the nomenclature and notation of servo-mech- 
anisms, and it is to be regretted that there are significant differences between 
this and the corresponding nomenclature of process cntrol as used by Farringdon. 
This may confuse the general reader who is not aware of the independent lines 
of development and the consequently different nomenclature used in these 
allied fields. : 

Farringdon deals solely with the theory of process control: in general the 
treatment is basically similar to that of MacMillan, but greater emphasis is 
naturally placed on the types of control and discontinuities normally encountered - 
in process applications. There is an excellent discussion of analogous transfer 
stages, and the system equations utilising proportional, integral and derivative 
controls. Plant characteristics, and the control discontinuities and non- 
linearities such as are found in process regulating units are also discussed. 
There is a very interesting analysis of typical pneumatic controlling intruments. 
The mathematical treatment is generally adequate and not beyond the 
capabilities of the graduate in physics or engineering. 

This book will be of greatest interest to those workers concerned with process 
control—to the designer rather than the operator or user, but it will amply 
repay study by those interested in the other related fields, if only to correlate 
the extremely wide fields of application. A. P. 


ERRATA. 


Total Cross-Sections of the Elements for 156 MeV. Neutrons, by A. E. Taytor, 
T. G. Pickavancrz, J. M. Casses, and T. C. Ranpuz, 1951, Phil. Mag., 
42, 751. 


The experimental results reported in this paper were discussed briefly in the 
light of the optical model of nuclei. The discussion involved the mean free 
paths of neutrons in nuclear matter, and these were incorrectly calculated 
because wrong values had been adopted for the proton—proton cross-section. 
The error arose from the fact that the total cross-section for proton—proton 
“7 p27 

scattering is | 27o0(@) sin @ dé, and not | 27o(8) sin @d@, where o(@) is the 
0 0 

differential cross-section as conventionally defined. 

Correction of this point removes most of the difficulties which were thought 
to follow from the assumption that the neutron—neutron and proton—proton 
interactions are the same at high energies. Some discrepancy still remains, 
as the theoretical results of Jastrow (1951) have shown. 


R. JastRow, 1951, Phys. Rev., 82, 1951. 


A Theoretical Investigation of the Compression of a Ductile Material between 
Smooth Flat Dies, by A. P. Grenn, 1951, Phil. Mag., 42, 917. 


Equation below fig. 13 should read : ‘ 
e=—v; u me : Ist kind, 


m. : 
y=u; v= 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


